Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,906,826 T→C 7.8% intergenic (‑259/+27) bglG ← / ← phoU transcriptional antiterminator of the bgl operon/negative regulator of PhoR/PhoB two‑component regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,906,8260TC 7.8% 181.6 / 5.5 64intergenic (‑259/+27)bglG/phoUtranscriptional antiterminator of the bgl operon/negative regulator of PhoR/PhoB two‑component regulator
Reads supporting (aligned to +/- strand):  ref base T (35/24);  new base C (3/2);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

AGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGC  >  NC_000913/3906542‑3907119
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                     
agaATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAACACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTACTTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATCAAGTTATATATAAAAAATCCAAATAATTAAGTTATTGGGAttt                                                                                                                                                                                                                                                                                       >  1:450225/1‑301 (MQ=255)
  aaTTTTGGTTATTGGGATGTTCATATTAAGGACC‑TTTTTATAAATAAAAAAACCCGACTTCACTAGTATTCTCTGTTTATGTCAGGGTTTGGCTGCGAATGCAGTAACAATCCAGTCAATTATTAATGGGTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGATTTGtt                                                                                                                                                                                                                                                                                    <  1:245769/301‑2 (MQ=255)
      ttGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGATATTGGGATTTGTCTgg                                                                                                                                                                                                                                                                                 >  1:219621/1‑301 (MQ=255)
          tGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGaa                                                                                                                                                                                                                                                                             >  2:142610/1‑301 (MQ=255)
              tttcATGTTCATAGCAAGGTCCTTTTTTATTAACAAAAAAACCCGAGTTCACCATTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACTATCCAGTCATTTATTAATGGGTTTTATAACGAACATCCAGGTTCGGAAATTAATTTACTTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAattat                                                                                                                                                                                                                                                                         <  1:570066/298‑1 (MQ=255)
               tGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTAATATATAACAAAACCAAATAATAAAGTTATTGGGATTTGTATGGTGAattatt                                                                                                                                                                                                                                                                        >  1:32848/1‑301 (MQ=255)
                   tGTTCATAGCAAGGACCTTTTTTTTAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTc                                                                                                                                                                                                                                                                    <  1:242570/301‑1 (MQ=255)
                    tttCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCg                                                                                                                                                                                                                                                                   <  2:450789/300‑1 (MQ=255)
                     ttCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGc                                                                                                                                                                                                                                                                  <  2:375324/301‑1 (MQ=255)
                           gCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCtt                                                                                                                                                                                                                                                            >  2:239768/1‑301 (MQ=255)
                                  ccTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGc                                                                                                                                                                                                                                                     >  1:474127/1‑301 (MQ=255)
                                   tttttttATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGcc                                                                                                                                                                                                                                                    <  2:32899/300‑1 (MQ=255)
                                           taAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTCTTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGCAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTACTGTTAGGCATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGttt                                                                                                                                                                                                                                            >  1:355843/1‑301 (MQ=255)
                                            aaaCAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGtttt                                                                                                                                                                                                                                           >  1:248162/1‑300 (MQ=255)
                                               cAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATcc                                                                                                                                                                                                                                        >  1:39685/1‑301 (MQ=255)
                                                 aataatctcGACTTTACTCGTGTGCTCTTGTTATGTCAGGTTTTGCCGTCGAATGCAGTAACATTCCAGTCATTTATTAATGGTGTTTATAACGAACATCCATGTTCGTAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAg                                                                                                                                                                                                                                      <  1:304625/293‑1 (MQ=255)
                                                 aaaaaaCCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAg                                                                                                                                                                                                                                      >  1:37221/1‑301 (MQ=255)
                                                      acctgacataaccagagaatactGGTGAAGTCGGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCaa                                                                                                                                                                                                                                 >  2:102450/22‑300 (MQ=255)
                                                       cccGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATc                                                                                                                                                                                                                                <  2:474728/301‑1 (MQ=255)
                                                               cACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATACTTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCgc                                                                                                                                                                                                                        >  1:389630/1‑300 (MQ=255)
                                                                                 ttGTCAGGTTTTGTCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTGTTTATAACGAACATCCAGGGTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGcc                                                                                                                                                                                                      <  2:356327/300‑1 (MQ=255)
                                                                                   gTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGcccc                                                                                                                                                                                                    >  2:94098/1‑301 (MQ=255)
                                                                                      gggTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTc                                                                                                                                                                                                 <  1:102326/300‑1 (MQ=255)
                                                                                            tGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACgtaa                                                                                                                                                                                           >  2:121960/1‑300 (MQ=255)
                                                                                               acgctcttccgatctcAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGtagtag                                                                                                                                                                                        <  1:431686/286‑1 (MQ=255)
                                                                                                 ctttccctactcgactctcttccgatctttaATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                                                                                                                                                                      <  1:555478/273‑1 (MQ=255)
                                                                                                              cAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagagatcggaagagcaa                                                                                                                                                                         >  2:432250/1‑287 (MQ=255)
                                                                                                                aTCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACaaaa                                                                                                                                                                       >  2:483384/1‑300 (MQ=255)
                                                                                                                aTCCAGTCATTTATTAATGGTTTTTATAACCAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGACAATAGCAATGCGCTATTGATAAAAATCTGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTCGGCATTTATAAAGTTATATATAACAAATCCCCATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCCGTTTCTCCAGCTACTCCCCCCCGACGTGACCGAAACCCCGCCCCTTCACGTAGTATACAAAAAACCCACaaaa                                                                                                                                                                       >  1:551697/1‑300 (MQ=255)
                                                                                                                        aTTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGc                                                                                                                                                               >  2:87718/1‑301 (MQ=255)
                                                                                                                             ttaATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaaagatcggaagagaacacgtctgaactcc                                                                                                                                                          >  2:556192/1‑273 (MQ=255)
                                                                                                                                  tgTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGcc                                                                                                                                                     <  2:37272/300‑1 (MQ=255)
                                                                                                                                   gTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATACTTAAGTTATTGGGATTTGTCTGGTGACTTATTTGTCGCTATCTTTCCCCGCCAGCCGTTTATCCACCTCATCGCCACCGACGTGACGGACATCCTGCCCCTTCACGTAGTAGACGATAAACCCACAAAAATCCtacccacggccaccc                                                                                                                                                    >  1:46593/1‑287 (MQ=255)
                                                                                                                                    tttttATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCaa                                                                                                                                                   <  2:39737/301‑1 (MQ=255)
                                                                                                                                      ttttTAACGAACATCCAGTTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATa                                                                                                                                                 <  2:502921/301‑1 (MQ=255)
                                                                                                                                           aaCGAACATCCAGGTTTGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTc                                                                                                                                            <  1:87594/301‑1 (MQ=255)
                                                                                                                                                   ttcatgTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCTATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAAc                                                                                                                                    <  1:121836/296‑1 (MQ=255)
                                                                                                                                                         ttCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAAAAcacacac                                                                                                                              >  1:92214/1‑295 (MQ=255)
                                                                                                                                                                 ttaatttaatTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGACCGATAGAACGCGCGCAGAACacc                                                                                                                      >  1:398539/1‑299 (MQ=255)
                                                                                                                                                                   tatttaattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAgcgc                                                                                                                    <  2:390144/300‑1 (MQ=255)
                                                                                                                                                                   aatttaattGCGTGCAGAGAAAATAGCAATGAGTTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATGTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCGAATACGTTCGATAGAACGCGCGCAGAACAgcgc                                                                                                                    <  1:482783/301‑1 (MQ=255)
                                                                                                                                                                           tGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCagtaagta                                                                                                            >  2:569714/1‑301 (MQ=255)
                                                                                                                                                                                  agagAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCgg                                                                                                     >  1:293948/1‑301 (MQ=255)
                                                                                                                                                                                  agagAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCgg                                                                                                     >  2:574836/1‑301 (MQ=255)
                                                                                                                                                                                                     tATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTcc                                                                                  <  2:219864/301‑1 (MQ=255)
                                                                                                                                                                                                                   tGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGtcat                                                                    <  2:248429/301‑1 (MQ=255)
                                                                                                                                                                                                                         tGCTCGCTGTTATTAACTTTTTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTg                                                              <  2:399053/301‑1 (MQ=255)
                                                                                                                                                                                                                           ttcggatttttttattttGTGTAATTTTAGGAATTTATTAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGGCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGAc                                                            <  1:239501/286‑1 (MQ=255)
                                                                                                                                                                                                                            tCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGAca                                                           >  2:393308/1‑300 (MQ=255)
                                                                                                                                                                                                                                  ttatgtaCTTTGTGTAATTTTAGGTATTTATAAAGTTTTATATAACAAATCCCAATATTTAATTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAAt                                                     <  1:142485/295‑1 (MQ=255)
                                                                                                                                                                                                                                   tattaACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATa                                                    >  1:457620/1‑301 (MQ=255)
                                                                                                                                                                                                                                     ttaACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCAACTCTTCCATCATATAGGTCATCAGCAGACGAACAACAcc                                                  >  1:97342/1‑301 (MQ=255)
                                                                                                                                                                                                                                       aaCTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAAAAATACCtt                                                >  1:599550/1‑301 (MQ=255)
                                                                                                                                                                                                                                              tgtAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTc                                         >  2:78158/1‑301 (MQ=255)
                                                                                                                                                                                                                                                    tttAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATc                                   <  2:600311/301‑1 (MQ=255)
                                                                                                                                                                                                                                                              tcgtgtgctcttccgatctaaTCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTa                         <  2:88474/282‑1 (MQ=255)
                                                                                                                                                                                                                                                               ataAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTAt                        >  1:84948/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                              aCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCGCGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATaa         >  2:422876/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                              aCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCCTCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATAc         >  1:105684/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                aaaTCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGAAAAATACGt       >  1:492839/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                 aaTCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTaagatcggaagagcgtcgtt      >  1:88350/1‑282 (MQ=255)
                                                                                                                                                                                                                                                                                   tCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTAcc    >  2:403372/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                   tCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGAAAAATACGTAcc    >  2:290487/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                     ccAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGACGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACCATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGc  >  1:190335/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                     
AGAATTTTGGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAGTAACAATCCAGTCATTTATTAATGGTTTTTATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGC  >  NC_000913/3906542‑3907119

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑