Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,750,268 C→T 6.4% P27P (CCG→CCA ydhU ← putative cytochrome b subunit of YdhYVWXUT oxidoreductase complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,750,2680CT 6.4% 189.5 / 2.3 59P27P (CCG→CCAydhUputative cytochrome b subunit of YdhYVWXUT oxidoreductase complex
Reads supporting (aligned to +/- strand):  ref base C (31/24);  new base T (2/2);  total (33/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 = 7.72e-01

AGCCAGCATGCCAGTAACAAAAATCCGCAAACTTCATGCACCGCAACCAGACTTTTAACCGCAGTTGCGCCCACCATCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGA  >  NC_000913/1749981‑1750556
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                
aGCCAGCATGCCAGTAACAAAAATCCGCAAACTTCATGCACCGCAACCAGACTTTTAACCGCAGTTGCGCCCACCATCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCg                                                                                                                                                                                                                                                                                                 >  2:247321/1‑289 (MQ=255)
                aCAAAAATCCGCAAACTTCATGCACCGCAACCAGACTTTTAACCGCAGTTGCGCCCACCATCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGt                                                                                                                                                                                                                                                                                 >  2:289368/1‑289 (MQ=255)
                        ccGCAAACTTCATGCACCGCAACCAGACTTTTAACCGCAGTTGCGCCCACCATCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACtgtg                                                                                                                                                                                                                                                                         >  2:136468/1‑289 (MQ=255)
                           cAAACTTCATGCACCGCAACCAGACTTTTAACCGCAGTTGCGCCCACCATCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGACCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGc                                                                                                                                                                                                                                                                       >  2:230314/1‑288 (MQ=255)
                                          gCAACCAGACTTTTAACCGCAGTTGCGCCCACCATCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTa                                                                                                                                                                                                                                                          >  1:160571/1‑286 (MQ=255)
                                                   cgtttAACCCAAGTTGCGCCACCCACCGAAAAGGGATTAACCAGCCCACGGGCAAGCAAAAAACCAAAGGGTAACCCATCCAACCATTGCCTAAACTTGCCCGTTTGGTTGTAAAGTAAAACATTCTAACAAAGACGGGGGGCGGTTTTCTCCCCGCCACGGCCACGAACAAAGCCATCCGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGaaa                                                                                                                                                                                                                                             <  2:366207/287‑1 (MQ=255)
                                                     tttAACCGCAGTTGCGCCCACCATCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGCTCGCTAACTGGCTCTGaa                                                                                                                                                                                                                                              >  1:706573/1‑287 (MQ=255)
                                                      ttAACCGCAGTTGCGCCCACCATCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACTGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGaa                                                                                                                                                                                                                                              >  1:546109/1‑286 (MQ=255)
                                                            gCAGTTGCGCCCACCATCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAAATGTGGCCAATAGTTCGCTAACTGGCCTTGAAACTGTTCa                                                                                                                                                                                                                                      >  2:20183/1‑288 (MQ=255)
                                                                      ccaccaTCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGa                                                                                                                                                                                                                           <  2:4547/289‑1 (MQ=255)
                                                                      ccaccaTCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGa                                                                                                                                                                                                                           <  1:289368/289‑1 (MQ=255)
                                                                      ccaccaTCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGa                                                                                                                                                                                                                           <  1:247321/289‑1 (MQ=255)
                                                                       caccaTCTCAAAGTGATGTATCACCCCACGTGCACGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCCTAATCTGACCGCTTTGCTGTAAAGATAAATCTGATCACCATGATCGGTGGCTGATTTCTTCACGCCACGATCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGAc                                                                                                                                                                                                                          <  2:424372/289‑1 (MQ=255)
                                                                                      aTTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACTGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGAc                                                                                                                                                                                                                          >  1:735276/1‑274 (MQ=255)
                                                                                      aTTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACTGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGAc                                                                                                                                                                                                                          <  2:735276/274‑1 (MQ=255)
                                                                                                    tGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATAc                                                                                                                                                                                                                >  1:604159/1‑270 (MQ=255)
                                                                                                         aGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCCCACCCGCCAAAATTCGGGAGTGAACTGTCGCACATAGTTCGCTAACTGGCCCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTCCCTTTTTAATCAAATGttt                                                                                                                                                                                        >  2:708730/1‑288 (MQ=255)
                                                                                                                aTACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTCCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCACAACATCGGTAGTGAACTGTGGCACAAAGTTCGCTAACTGGCTCTGAAACTGTCCAGCATGGAGCGACGGGTTCCTAGCTTCTCCTTTTTAATACGATGTTTACCaaa                                                                                                                                                                                  >  2:562263/1‑288 (MQ=255)
                                                                                                                        agTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATa                                                                                                                                                                         >  1:470475/1‑289 (MQ=255)
                                                                                                                                  tCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACAt                                                                                                                                                                                                                                                                   <  2:700644/189‑1 (MQ=255)
                                                                                                                                  tCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACAt                                                                                                                                                                                                                                                                   >  1:700644/1‑189 (MQ=255)
                                                                                                                                  tCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGt                                                                                                                                                               >  2:15220/1‑289 (MQ=255)
                                                                                                                                           gCCATAAACTGACCGATTAGCGGTACAGATAAATCTTCTCACCATGGTCGGTGGCTGATTCCGTCACGCCACGAGCCGGAAGCAAGGCATGCAGCCCCAGCACCAACCGCATCCCACCCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAGATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCa                                                                                                                                                      <  2:706573/289‑1 (MQ=255)
                                                                                                                                                 aTCTGACCGCTTAGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGt                                                                                                                                                <  1:729962/289‑1 (MQ=255)
                                                                                                                                                    tGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTa                                                                                                                                               >  1:170170/1‑287 (MQ=255)
                                                                                                                                                    tGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCAACAGGTa                                                                                                                                               >  1:226939/1‑287 (MQ=255)
                                                                                                                                                     gACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAg                                                                                                                                              <  2:170170/287‑1 (MQ=255)
                                                                                                                                                      aCCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAACATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTg                                                                                                                                           >  2:147347/1‑289 (MQ=255)
                                                                                                                                                      aCCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCa                                                                                                                                                                                                                   >  1:358999/1‑217 (MQ=255)
                                                                                                                                                      aCCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCa                                                                                                                                                                                                                   <  2:358999/217‑1 (MQ=255)
                                                                                                                                                          cTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGa                                                                                                                                                                                                                           >  1:241019/1‑205 (MQ=255)
                                                                                                                                                          cTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGa                                                                                                                                                                                                                           <  2:241019/205‑1 (MQ=255)
                                                                                                                                                             tgatttAAAGTTAAACCTGAGCACCATGATCGCTGGCCTATATCTGCACGCCGCGAGACCGAAGCAAGGCATGCAGCCCCAGCACCAACCACAGCCCAACCAGCAATACTCCGGCAAGGATCAACCACAGCGGCCAAAATTCCGGAGTGAACTGTGGCACATCTCTCGCTAAGTGGCTCTGAAACAGAAGAGCCAGTTGCGACGTGTTCATACATTTTCCTTTTTAATCAAATGTTGATCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAAGTAACTGATATTGa                                                                                                                                    <  2:280789/284‑1 (MQ=255)
                                                                                                                                                              gCTGTAAAGATAAATCTGCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGAt                                                                                                                                   <  1:610666/289‑1 (MQ=255)
                                                                                                                                                               cTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCACCGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGAt                                                                                                                                   <  2:457069/288‑1 (MQ=255)
                                                                                                                                                                 gTAAAGATAAATCTTGTCACCATGATCGGTGGCTGATTTCTGCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCATACACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGt                                                                                                                                <  1:630501/289‑1 (MQ=255)
                                                                                                                                                                 gTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGt                                                                                                                                <  1:136468/289‑1 (MQ=255)
                                                                                                                                                                         aaaTCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTatttattt                                                                                                                        >  1:66646/1‑289 (MQ=255)
                                                                                                                                                                                             gTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCgg                                                                                                                                                        <  1:36553/237‑1 (MQ=255)
                                                                                                                                                                                             gTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCgg                                                                                                                                                        >  2:36553/1‑237 (MQ=255)
                                                                                                                                                                                               ggCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAAtt                                                                                                   >  1:286286/1‑288 (MQ=255)
                                                                                                                                                                                                     tttCGTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTggg                                                                                             <  1:20183/288‑1 (MQ=255)
                                                                                                                                                                                                              gCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACTGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTcgc                                                                                   <  2:546109/289‑1 (MQ=255)
                                                                                                                                                                                                                         gAAGCAAGGCATTCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATc                                                                        <  1:667550/289‑1 (MQ=255)
                                                                                                                                                                                                                              aaGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGc                                                                   <  2:66646/289‑1 (MQ=255)
                                                                                                                                                                                                                                          accAGAACCTACCACATCCCACCCAGCAATACTCCGTCAATGATCAACCACATCAGCCAGCATTCGGGAGTGCACTGTGGCCCATAGTTCGCTCACTGGCTCTGACACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATATTGATCGAAGCGGCGATACACACGAGGTTCCGCGGCACCAGGTAGATGATATCGCTAGTATTCATTTTGCTGCACCCACGCCTGATTTTCTGGGCTATCTTGGCGACCACAGATCAATTCGTGTCCCGGGCa                                                       >  2:95182/2‑289 (MQ=255)
                                                                                                                                                                                                                                                accaaccaCATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGgcg                                                   <  1:15220/287‑1 (MQ=255)
                                                                                                                                                                                                                                                         aTCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAATTTCGGGAGTGAACTGTGCCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACAGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGAAACAAATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCCACCACGCCTGAATTTCTAGGCTATCTTCGCGACGAATGCTCAATGCGGGTTCCGGGCAGgctacgaacacaat                                        >  2:284336/1‑277 (MQ=255)
                                                                                                                                                                                                                                                         aTCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGaaa                                                                                                                                                                                                                                             <  1:711025/92‑1 (MQ=255)
                                                                                                                                                                                                                                                         aTCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGaaa                                                                                                                                                                                                                                             >  2:711025/1‑92 (MQ=255)
                                                                                                                                                                                                                                                               aCCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCgg                                  >  2:197542/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                    cAATACTCCCGCCATGATCAACCACAACGGCCAAAATTCGGGAGTGCACTGTGGCCAATAGTTCCCAAACTGGCTCTGAAACTGTTCAGCAGCTTGCGACGGATACGTACATTTTCCTTTTTAATCACATGTTCCTCGAACTGGCGATCCAGATGCAGTTTCTCGGCCCCAGGTAGCCGATACTGATAGTATTTATTTTGCTGCAACCACGCCTGAATATCTGGTCTATGTTCGCGACCAGAGATCAATGACTTTCCCGGTCAGGCCCTGACGCAAATgtgtcgaacg                              >  2:245438/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                          tCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTg                       <  1:197542/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATTCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAg                   >  2:701391/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                              gCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAg                   >  2:370525/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                              gCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAg                   >  2:223642/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                            cacCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACt           >  2:455021/1‑283 (MQ=255)
                                                                                                                                                                                                                                                                                            cacCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACt           <  1:455021/283‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            cacCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACTCAgcg      >  2:67269/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                               cGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGa  >  2:439014/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                
AGCCAGCATGCCAGTAACAAAAATCCGCAAACTTCATGCACCGCAACCAGACTTTTAACCGCAGTTGCGCCCACCATCGCAAAGTGATTTATCAGCCCACTGGCCAGCAACAATACAAAGAGTAACGCATTCGACCAGTGCCATAATCTGACCGCTTTGCTGTAAAGATAAATCTTCTCACCATGATCGGTGGCTGATTTCTTCACGCCACGAGCACGAAGCAAGGCATGCAGCCCCAGCACCAACCACATCCCAACCAGCAATACTCCGGCAATGATCAACCACACCGGCCAAAATTCGGGAGTGAACTGTGGCACATAGTTCGCTAACTGGCTCTGAAACTGTTCAGCATGTTGCGACGGGTTCATACATTTTCCTTTTTAATCAAATGTTGACCGAACCGGCGATACAGATGCGGTTTCCCGGCACCAGGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGA  >  NC_000913/1749981‑1750556

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