Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,511,669 G→A 5.9% A6T (GCC→ACC)  ydcS → putative ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,511,6690GA 5.9% 230.2 / 4.0 68A6T (GCC→ACC) ydcSputative ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (34/30);  new base A (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TTCAAGAACGGCTGTAATAGCGTTTAATTTAATTCCTCTTAGATTGGGTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGA  >  NC_000913/1511393‑1511957
                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                
ttCAAGAACGGCTGTAATAGCGTTTAATTTAATTCCTCTTAGATTGGGTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGcagcag                                                                                                                                                                                                                                                                                        >  2:332733/1‑287 (MQ=255)
        cGGCTGTAATAGCGTTTAATTTAATTCCTCTTAGATTTGGTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTgcg                                                                                                                                                                                                                                                                              <  1:360661/289‑1 (MQ=255)
                     gTTTAATTTAATTCCTCTTAGATTGGGTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTATACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGccc                                                                                                                                                                                                                                                                                               <  1:771944/259‑1 (MQ=255)
                     gTTTAATTTAATTCCTCTTAGATTGGGTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGccc                                                                                                                                                                                                                                                                                               >  2:771944/1‑259 (MQ=255)
                     gTTTAATTTAATTCCTCTTAGATTGGGTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCgcag                                                                                                                                                                                                                                                                                           >  2:135380/1‑263 (MQ=255)
                     gTTTAATTTAATTCCTCTTAGATTGGGTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCgcag                                                                                                                                                                                                                                                                                           <  1:135380/263‑1 (MQ=255)
                        taatttaattCCTCTTAGATTGGGTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTACCCGCAGCAGCCTGTGTGCGCTCAGCATGACaataa                                                                                                                                                                                                                                                              >  2:579030/1‑289 (MQ=255)
                                    tctTAGATTGGGTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCAcg                                                                                                                                                                                                                                                  <  2:479269/289‑1 (MQ=255)
                                              ggTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAAcc                                                                                                                                                                                                                                        >  2:194380/1‑289 (MQ=255)
                                              ggTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAAc                                                                                                                                                                                                                                         >  1:485795/1‑288 (MQ=255)
                                                        aTTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGCCATTTGCCCGCAACATCCTGTGTGCGGCCATCATGACAATAATGACCGCTCACGCCGCCGAcccccgcacaaat                                                                                                                                                                                                                              >  1:644242/1‑279 (MQ=255)
                                                          ttCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCgcag                                                                                                                                                                                                                                                                                           <  1:11412/226‑1 (MQ=255)
                                                          ttCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCgcag                                                                                                                                                                                                                                                                                           >  2:11412/1‑226 (MQ=255)
                                                                             tcgtcaCTCCTTTACTATACTCCAGAAGATAGCCTTACAGACGGCATAATGCGCGATACCTCACAACCTGAATAAAGTTTCTCAGGTGCAAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGAAAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGa                                                                                                                                                                                                         <  2:644242/284‑1 (MQ=255)
                                                                             tCCTCACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGa                                                                                                                                                                                                         <  2:618787/289‑1 (MQ=255)
                                                                             tCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGa                                                                                                                                                                                                         <  2:95619/289‑1 (MQ=255)
                                                                              ccTAAATCCTTGACTATACTCCAGAAGATAACCTTACTAGCGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTATGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGa                                                                                                                                                                                                         <  1:293010/288‑1 (MQ=255)
                                                                                tAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATa                                                                                                                                                                                                                       <  1:126121/272‑1 (MQ=255)
                                                                                tAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATa                                                                                                                                                                                                                       >  2:126121/1‑272 (MQ=255)
                                                                                      cTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATa                                                                                                                                                                                                                       >  2:594685/1‑266 (MQ=255)
                                                                                      cTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATa                                                                                                                                                                                                                       <  1:594685/266‑1 (MQ=255)
                                                                                        tGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGATCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCCCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGAtatta                                                                                                                                                                                               >  2:293010/1‑288 (MQ=255)
                                                                                        tGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTa                                                                                                                                                                                                                                   <  2:36609/252‑1 (MQ=255)
                                                                                        tGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTa                                                                                                                                                                                                                                   >  1:36609/1‑252 (MQ=255)
                                                                                        tGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGAtatta                                                                                                                                                                                               >  1:39043/1‑288 (MQ=255)
                                                                                              taCTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAg                                                                                                                                                                                                                                                                         >  1:150369/1‑208 (MQ=255)
                                                                                              taCTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAg                                                                                                                                                                                                                                                                         <  2:150369/208‑1 (MQ=255)
                                                                                                              cTTACAGACGGCAGAATGCGCTGTAGCTAACAACCTGAATAAATTTTCTCAGGGGCGAAGTGGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGTGATACATCGaa                                                                                                                                                                        <  2:418257/289‑1 (MQ=255)
                                                                                                                tACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAAc                                                                                                                                                                       >  1:772449/1‑288 (MQ=255)
                                                                                                                      cGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACa                                                                                                                                                                <  2:39043/289‑1 (MQ=255)
                                                                                                                      cGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACa                                                                                                                                                                <  2:772449/289‑1 (MQ=255)
                                                                                                                      cGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACa                                                                                                                                                                <  1:258170/289‑1 (MQ=255)
                                                                                                                               gcgcGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATaa                                                                                                                                                       >  1:731734/1‑289 (MQ=255)
                                                                                                                                   ggTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACaa                                                                                                                                                   >  1:26141/1‑289 (MQ=255)
                                                                                                                                   ggTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACaa                                                                                                                                                   >  1:620521/1‑289 (MQ=255)
                                                                                                                                        cTCACAACCTGAATAAAGTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                                                                                              <  1:728223/289‑1 (MQ=255)
                                                                                                                                         tCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                                                                                              <  1:634728/288‑1 (MQ=255)
                                                                                                                                             aaCCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTggg                                                                                                                                         <  2:128244/289‑1 (MQ=255)
                                                                                                                                             aaCCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTggg                                                                                                                                         >  1:128244/1‑289 (MQ=255)
                                                                                                                                                                            gtgCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTACCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTgcg                                                                                                           <  1:579030/288‑1 (MQ=255)
                                                                                                                                                                             tgCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATCAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTgcgag                                                                                                         >  2:199294/1‑287 (MQ=255)
                                                                                                                                                                             tgCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGAc                                                                                                                   <  2:7842/279‑1 (MQ=255)
                                                                                                                                                                             tgCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGAc                                                                                                                   >  1:7842/1‑279 (MQ=255)
                                                                                                                                                                              gCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                                                                                              >  2:652412/1‑251 (MQ=255)
                                                                                                                                                                              gCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                                                                                              <  1:652412/251‑1 (MQ=255)
                                                                                                                                                                                 tGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTACCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGa                                                                                                     >  1:64656/1‑289 (MQ=255)
                                                                                                                                                                                  gCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGa                                                                                                                                                                                                         <  2:296207/188‑1 (MQ=255)
                                                                                                                                                                                  gCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGa                                                                                                                                                                                                         >  1:296207/1‑188 (MQ=255)
                                                                                                                                                                                          ccgTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGATCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGaaa                                                                                             >  2:386222/1‑288 (MQ=255)
                                                                                                                                                                                                       aaCAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAg                                                                                                                 <  2:207873/255‑1 (MQ=255)
                                                                                                                                                                                                       aaCAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAg                                                                                                                 >  1:207873/1‑255 (MQ=255)
                                                                                                                                                                                                       aaCAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACtt                                                                                >  2:164765/1‑288 (MQ=255)
                                                                                                                                                                                                       aaCAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTc                                                                               <  2:731734/289‑1 (MQ=255)
                                                                                                                                                                                                         cAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCg                                                                             <  2:485795/289‑1 (MQ=255)
                                                                                                                                                                                                                                  ggcTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTACCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCa                                                     <  2:64656/288‑1 (MQ=255)
                                                                                                                                                                                                                                              aattaCGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGa                                        >  2:99687/1‑289 (MQ=255)
                                                                                                                                                                                                                                               attaCGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGAt                                       <  1:332733/289‑1 (MQ=255)
                                                                                                                                                                                                                                                        ggAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACg                              >  2:306931/1‑289 (MQ=255)
                                                                                                                                                                                                                                                          aGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTTCGCTCAGCATGACAAAAATGACCGCTCAGGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACg                              <  1:756952/287‑1 (MQ=255)
                                                                                                                                                                                                                                                              aGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACGGCCGCGACTTCCGATGAAATGGGCAGTCTGATGCCCAAAGGGGGTTACGCTCTGGTTACGGAAtac                        >  1:108611/1‑287 (MQ=255)
                                                                                                                                                                                                                                                              aGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGAACTGGTTACGGCATcc                        >  2:428815/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                    tATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTTCGCTCAGCATGACAATAATGACAGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAAGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGAAAAACAATACGACTGGGTAACGCAAGTCAAAAAAGAGACAGTCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGACATGGGCAGTCTGATGACGAAAGGGGGTTGGGAGCTGGTTACggatacgggcggc                  >  1:500637/1‑278 (MQ=255)
                                                                                                                                                                                                                                                                    tATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACg                              >  1:398650/1‑277 (MQ=255)
                                                                                                                                                                                                                                                                    tATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACg                              <  2:398650/277‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGAt                  >  2:418770/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                    tATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGAt                  >  2:406067/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                    tATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGAt                  <  2:16232/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACAGGACGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCt          >  2:740840/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                    gCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGa  <  1:99687/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                
TTCAAGAACGGCTGTAATAGCGTTTAATTTAATTCCTCTTAGATTGGGTAATATGAATTTCGAATAGCAGTCATATTTCCTAACTCCTTGACTATACTCCAGAAGATAACCTTACAGACGGCATAATGCGCGGTAGCTCACAACCTGAATAAATTTTCTCAGGGGCGAAGGTGTGCCTGCAAGCCGCCGTCTATGGTTAAACAAGGAGATATTTTTACGGCACGGCGGCTGAACAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGA  >  NC_000913/1511393‑1511957

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