Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,563,690 C→T 5.7% intergenic (+172/+97) glpD → / ← yzgL sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding/pseudogene, periplasmic solute binding protein homology

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,563,6900CT 5.7% 399.9 / 9.4 125intergenic (+172/+97)glpD/yzgLsn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding/pseudogene, periplasmic solute binding protein homology
Reads supporting (aligned to +/- strand):  ref base C (61/55);  new base T (3/4);  total (65/60)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

Only 100 of 126 total aligned reads displayed.

GATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGACACCCAGCATGGCAATAACCACCGCCAAAAATTGCGCCAGTATGGGGATGCCGAAAA  >  NC_000913/3563450‑3563962
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                
gATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAAtt                                                                                                                                                                                                                                                >  2:72816/1‑275 (MQ=255)
              tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGNTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACNTTCACGTCGAATCCGGCATGAACAAAGAGCACTTTGTTAAAAATCNANACTACTTTTTAATTTAAGGAAATaacac                                                                                                                                                                                                                                  >  3:33493/1‑271 (MQ=255)
              tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCg                                                                                                                                                                                                                                   >  2:43472/1‑274 (MQ=255)
               gtgGGTAAGTGGCGGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGc                                                                                                                                                                                                                                  <  4:467159/274‑1 (MQ=255)
                  aGTCAGTGGCTGCGGGAGTATAAGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATAAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAg                                                                                                                                                                                                                              <  2:13590/275‑1 (MQ=255)
                             ggtggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTCTCAGGCCTACGTGGTTNATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATTAACAAATCGNACTTTGTCCAACAACCAACCTACTTTTTAANNTATGGAATTTCCCCACACCCCaaccac                                                                                                                                                                                                                     >  3:69779/1‑268 (MQ=255)
                                   gTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATAcc                                                                                                                                                                                                              <  1:58355/274‑1 (MQ=255)
                                                  gTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACCAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACt                                                                                                                                                                                               >  2:21420/1‑274 (MQ=255)
                                                    tATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAg                                                                                                                                                                                            <  4:43471/275‑1 (MQ=255)
                                                     aTCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAACCACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAg                                                                                                                                                                                            >  2:105122/1‑274 (MQ=255)
                                                                                   aTCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGcc                                                                                                                                                                                                                            >  1:135080/1‑212 (MQ=255)
                                                                                       ggTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACg                                                                                                                                                                                                                        <  3:135080/212‑1 (MQ=255)
                                                                                                                        aaaCACAAAACTGCCTGATGCGCTACTCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGc                                                                                                                        <  1:69779/275‑1 (MQ=255)
                                                                                                                                     ccTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCTCCAGCATAAAAATTTTGAGCAGAAGGTAAAttg                                                                                                           >  2:332703/1‑275 (MQ=255)
                                                                                                                                     ccTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAAttg                                                                                                           >  4:91337/1‑275 (MQ=255)
                                                                                                                                         aTGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCa                                                                                                                                                                                                                                 <  4:20723/153‑1 (MQ=255)
                                                                                                                                              cTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGTGCGCTTTGttca                                                                                                                                                                                                                                                                      <  7:928961/111‑4 (MQ=255)
                                                                                                                                               tACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGt                                                                                                                                                                                                                                                                         <  7:381371/107‑1 (MQ=255)
                                                                                                                                               tACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGt                                                                                                                                                                                                                                                                         >  8:381371/1‑107 (MQ=255)
                                                                                                                                               tACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGAt                                                                                                 <  1:33493/275‑1 (MQ=255)
                                                                                                                                               cactcgtgtCAGGTCAACGTGATTTCTGCAATATATTGAATTTGCATGGTCATGGAGGCCCGCTAAGACCTTCACGTCGCATCCGGCATGAACAAGTTGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCACACTCCATCATACCCAACGTAACAAAACTCAGGCAAGTTCCCCACGTGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGAt                                                                                                 <  4:265685/267‑1 (MQ=255)
                                                                                                                                                     tATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTa                                                                                                                                                                                                                                                               <  6:1041614/111‑1 (MQ=255)
                                                                                                                                                     tATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTa                                                                                                                                                                                                                                                               >  8:486829/1‑111 (MQ=255)
                                                                                                                                                      aTCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGata                                                                                          >  2:155755/1‑275 (MQ=255)
                                                                                                                                                       tCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACt                                                                                                                                                                                                                                                                             <  6:357223/95‑1 (MQ=255)
                                                                                                                                                       tCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACt                                                                                                                                                                                                                                                                             >  5:357223/1‑95 (MQ=255)
                                                                                                                                                        cAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACttt                                                                                                                                                                                                                                                      <  6:737589/117‑1 (MQ=255)
                                                                                                                                                        cAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACttt                                                                                                                                                                                                                                                      >  5:737589/1‑117 (MQ=255)
                                                                                                                                                        cAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGc                                                                                                                                                                                                                       >  7:703574/1‑148 (MQ=255)
                                                                                                                                                           gCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTAc                                                                                      >  2:334374/1‑274 (MQ=255)
                                                                                                                                                              tACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGt                                                                                                                                                                                                                                                                         >  5:222269/1‑92 (MQ=255)
                                                                                                                                                              tACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTc                                                                                                                                                                                                                                               >  5:68856/1‑118 (MQ=255)
                                                                                                                                                              tACGTGGTGTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAAc                                                                                  <  4:334373/275‑1 (MQ=255)
                                                                                                                                                                  tGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAAt                                                                                                                                                                                                                                                 <  6:939494/112‑1 (MQ=255)
                                                                                                                                                                  tGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAAt                                                                                                                                                                                                                                                 >  5:939494/1‑112 (MQ=255)
                                                                                                                                                                  gggttttataaaTATATTGAAGTTTCACGGTTGTGTAGGCCAGTTAAGACGTTCCCGTCGCATAAGGCATGAACAAAGCGCACTTTGTCAAAACTCTACCCTACTTTTTAATACAGTGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACGAAACTCAGGCGCGTTCACCACTCGCCCCGAGAAAAAGAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCATAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCa                                                                              <  4:526224/265‑1 (MQ=255)
                                                                                                                                                                          gCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCaataat                                                                         <  2:91338/272‑1 (MQ=255)
                                                                                                                                                                             atatatTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAAt                                                                                                                                                               >  2:71455/1‑183 (MQ=255)
                                                                                                                                                                              tatatTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTTCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAAAAGGTAAATTGTTTGCGTGCTAAGGATATTACAAACTTCAATAATACCGGCa                                                                  >  2:295922/1‑275 (MQ=255)
                                                                                                                                                                              tatatTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGa                                                                                                                          >  2:285770/1‑219 (MQ=255)
                                                                                                                                                                               atatTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTc                                                                                                                                                                                                                                               <  7:763617/101‑1 (MQ=255)
                                                                                                                                                                               atatTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTc                                                                                                                                                                                                                                               >  8:763617/1‑101 (MQ=255)
                                                                                                                                                                                  ttGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCa                                                                                                                                                                                             <  8:703574/148‑1 (MQ=255)
                                                                                                                                                                                  ttGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAg                                                                                                                      <  4:285769/219‑1 (MQ=255)
                                                                                                                                                                                   tGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCt                                                                                                                                                                                                                                                           <  7:263638/85‑1 (MQ=255)
                                                                                                                                                                                   tGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCt                                                                                                                                                                                                                                                           >  8:263638/1‑85 (MQ=255)
                                                                                                                                                                                   tGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTAc                                                                                                                                                                                                                                     <  5:705327/107‑1 (MQ=255)
                                                                                                                                                                                   tGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTAc                                                                                                                                                                                                                                     >  6:705327/1‑107 (MQ=255)
                                                                                                                                                                                       tttGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCa                                                                                                                                                                                                                                              <  5:389152/94‑1 (MQ=255)
                                                                                                                                                                                          gCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAACTCTAACCTACTTTTTAATTCAGGGAATTACCGCCAAGCCCACTCCTTCATACCCACTGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATTCTCCCAGCGCCAGCATANAAATTTTGATACGAAGGTAAATTGTGCCCGTCAGAGGGATATAAAAAACTTCAATAATACCAGCCAATCCCAAAcc                                                       >  3:117219/1‑274 (MQ=255)
                                                                                                                                                                                          gCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACt                                                                                                                                                                              >  4:324161/1‑155 (MQ=255)
                                                                                                                                                                                          gCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTcc                                                             >  2:244603/1‑268 (MQ=255)
                                                                                                                                                                                          gCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCGAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCAAGCATAAAAATTTCGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGACAccc                                                      >  2:353044/1‑275 (MQ=255)
                                                                                                                                                                                          gCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTCCGTCAAAAACCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCTCACGCCATCATACTCAACGACACAAAACTCAGGCACGTTCCCGACTCGCGCCGAGAAAAATAGCATTAATGCGCCAACCGCCAGCATAAAAATTTTGAGAAGAAGGTACATTTTTGGCATTATATGGAATTTACAAACTTCAATAATACAGGCAAGTCCGAAAcc                                                       >  4:412042/1‑274 (MQ=255)
                                                                                                                                                                                             tGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATAcc                                                                                                                                                                                                              >  8:425976/1‑120 (MQ=255)
                                                                                                                                                                                             tGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATAcc                                                                                                                                                                                                              <  7:425976/120‑1 (MQ=255)
                                                                                                                                                                                              ggTCTTTTCGGACAGATAAAACGTTAACCGCGCATCAGGCATGAACAAAGCTCACTTTTCCAAAAATCTAACCTACCTTTTAATACAGGGAATTACCGCAAAGCCCACGCACTCATACCCAACGTAACAAAACTCAGGCACGGTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGACACCCAGCa                                                  <  4:21420/275‑1 (MQ=255)
                                                                                                                                                                                              ggTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGcc                                                                                                                                                                          <  2:324162/155‑1 (MQ=255)
                                                                                                                                                                                              ggTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGaca                                                         <  4:244602/268‑1 (MQ=255)
                                                                                                                                                                                                 cTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCa                                                                                                                                                                                                                                              <  6:798028/84‑1 (MQ=255)
                                                                                                                                                                                                        gCCAGATAAGAAGTTAACGGCGCATCGGGCATGAACAAAGCGCACTTTGTCAAAAATCTATCCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACGCAGGCACGTTCCCCACGCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGACACCCAGCATGGCAATAac                                        <  4:72815/275‑1 (MQ=255)
                                                                                                                                                                                                         ccAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACtt                                                                                                                                                                                                                                                                            <  8:284921/46‑1 (MQ=255)
                                                                                                                                                                                                          cAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGa                                                                                                                          >  4:333231/1‑191 (MQ=255)
                                                                                                                                                                                                          cAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAAtt                                                                                                            >  4:97074/1‑205 (MQ=255)
                                                                                                                                                                                                          cAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGACAccc                                                      >  2:340401/1‑259 (MQ=255)
                                                                                                                                                                                                          cAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGACACCCAGCATGGCAATAacc                                       >  2:218335/1‑274 (MQ=255)
                                                                                                                                                                                                           aGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAg                                                                                                                                                                                                                                             <  6:244989/75‑1 (MQ=255)
                                                                                                                                                                                                           aGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAg                                                                                                                                                                                                                                             >  5:244989/1‑75 (MQ=255)
                                                                                                                                                                                                            gATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGAAACCCAGCATGGCAATAaccac                                     >  1:49063/1‑274 (MQ=255)
                                                                                                                                                                                                              tAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAg                                                                                                                      <  2:333232/191‑1 (MQ=255)
                                                                                                                                                                                                              tAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGACACCCAGCa                                                  <  4:340400/259‑1 (MQ=255)
                                                                                                                                                                                                               aaGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAAtt                                                                                                                                                                                                                                                >  7:438698/1‑68 (MQ=255)
                                                                                                                                                                                                                   cGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGt                                                                                                                                                                                                        >  6:39026/1‑104 (MQ=255)
                                                                                                                                                                                                                   cGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGt                                                                                                                                                                                                        <  5:39026/104‑1 (MQ=255)
                                                                                                                                                                                                                     ttCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGtt                                                                                                                                                                                     >  6:352479/1‑121 (MQ=255)
                                                                                                                                                                                                                        aCGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTccc                                                                                                                                                                                  <  8:585311/121‑1 (MQ=255)
                                                                                                                                                                                                                        aCGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTccc                                                                                                                                                                                  >  7:585311/1‑121 (MQ=255)
                                                                                                                                                                                                                         cGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGc                                                                                                                                                                                                                             >  6:1006031/1‑77 (MQ=255)
                                                                                                                                                                                                                         cGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGc                                                                                                                                                                                                                             <  5:1006031/77‑1 (MQ=255)
                                                                                                                                                                                                                         cGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGtt                                                                                                                                                                                     >  7:1013909/1‑117 (MQ=255)
                                                                                                                                                                                                                         cGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTcc                                                                                                                                                                                   <  8:685458/119‑1 (MQ=255)
                                                                                                                                                                                                                         cGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTcc                                                                                                                                                                                   >  7:685458/1‑119 (MQ=255)
                                                                                                                                                                                                                           tCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCgaga                                                                                                                                                                    >  6:692875/1‑132 (MQ=255)
                                                                                                                                                                                                                           tCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCgaga                                                                                                                                                                    <  5:692875/132‑1 (MQ=255)
                                                                                                                                                                                                                            cccaccctggctaaaCAAAGCGCACTTTGTCAAAAATCTAAACACCTCTTTAAATCATGCAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAACACTCCGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACCTCAATAATACCGGCAAGTCCGACACCCAGCATGGCAATAACCACCGCCAAAAATTgcgc                      <  4:218334/260‑1 (MQ=255)
                                                                                                                                                                                                                            cGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGcc                                                                                                                                                                                                                      <  7:897135/81‑1 (MQ=255)
                                                                                                                                                                                                                            cGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGcc                                                                                                                                                                                                                      >  8:897135/1‑81 (MQ=255)
                                                                                                                                                                                                                                tCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTc                                                                                                                                                                                                                                               >  5:825967/1‑52 (MQ=255)
                                                                                                                                                                                                                                tCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCaa                                                                                                                                                                                                           <  6:566982/88‑1 (MQ=255)
                                                                                                                                                                                                                                 ccgacctGATCAACGGGCACTTTGTCAATAATCTAACCTACTTTTTCATCCCTGGAAGTACCGCAAAGGCCACGCCATCATATCCAACGTGACAAACTTCATGCACTTTCCCCGCACGCCCCGAGAAAAATAGCATTAATGCGCCCGGCGCCAGCATAAAAAGTTTGCGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGACACCCAGCATGGCAATAACCACCGCCAAAAATTGCGCCAGTa                 <  4:295921/267‑1 (MQ=255)
                                                                                                                                                                                                                                     cATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACtttt                                                                                                                                                                                                                                                     >  6:25372/1‑41 (MQ=255)
                                                                                                                                                                                                                                     cATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACtttt                                                                                                                                                                                                                                                     <  5:25372/41‑1 (MQ=255)
                                                                                                                                                                                                                                        tgaCAAAGTGCGCTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACaaa                                                                                                                                                                                                  <  6:415167/87‑1 (MQ=255)
                                                                                                                                                                                                                                        tgaCAAAGTGCGCTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACaaa                                                                                                                                                                                                  >  5:415167/3‑89 (MQ=255)
                                                                                                                                                                                                                                          aCAAAGTGCGCTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATa                                                                                                                                                                                                                >  7:549548/1‑73 (MQ=255)
                                                                                                                                                                                                                                          aCAAAGTGCGCTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATa                                                                                                                                                                                                                <  8:549548/73‑1 (MQ=255)
                                                                                                                                                                                                                                           cAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATa                                                                                                                                    >  6:651099/1‑148 (MQ=255)
                                                                                                                                                                                                                                            aaaGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATaa                                                                                                                                   >  7:623242/1‑148 (MQ=255)
                                                                                                                                                                                                                                              aGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGACACCCAGCATGGCAATAACCACCGCCAAAAATTGCGCCAGTATGGGGATGCCGaaaa  >  2:217834/1‑275 (MQ=255)
                                                                                                                                                                                                                                                cgcACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAAt                                                                                                                               >  7:339365/1‑148 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                
GATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAATTTTGAGCAGAAGGTAAATTGTTGGCGTGATAAGGATATTACAAACTTCAATAATACCGGCAAGTCCGACACCCAGCATGGCAATAACCACCGCCAAAAATTGCGCCAGTATGGGGATGCCGAAAA  >  NC_000913/3563450‑3563962

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: