Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 222,958 Δ1 bp 100% coding (126/576 nt) gmhB → D,D‑heptose 1,7‑bisphosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913222,9580A.99.5% 976.2 / ‑2.2 198coding (126/576 nt)gmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (80/117);  total (80/118)
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.82e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 222 total aligned reads displayed.

CCTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTATT  >  NC_000913/222811‑223103
                                                                                                                                                   |                                                                                                                                                 
ccTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaa                                                                                                                                                 <  2:757460/150‑1 (MQ=255)
ccTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaa                                                                                                                                                 <  2:62696/150‑1 (MQ=255)
ccTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaa                                                                                                                                                 <  2:58788/150‑1 (MQ=255)
ccTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaa                                                                                                                                                 <  1:1998602/150‑1 (MQ=255)
ccTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaa                                                                                                                                                 <  1:1631973/150‑1 (MQ=255)
ccTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTACAGGAGATCGAAAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaa                                                                                                                                                 >  2:1245188/1‑150 (MQ=255)
    attaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaa                                                                                                                                             <  2:1033726/150‑1 (MQ=255)
    attaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaa                                                                                                                                             <  2:903926/150‑1 (MQ=255)
               ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcg                                                                                                                                 <  1:468796/150‑1 (MQ=255)
               ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTgcg                                                                                                                                 <  1:768297/150‑1 (MQ=255)
                  cccgggggCGCAGCGCGGCCCCGCCATTTTTCGTGACCGTGATGGCACCATTAATGTCGCTCACGGCTATGTCCATGCGCTCGACAACTTGGTATTTATCGTCGGGGTTATTGCCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtg                                                                                                                              <  1:233380/144‑1 (MQ=255)
                    aGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggt                                                                                                                            >  1:1051869/1‑150 (MQ=255)
                     ggTGGCGACGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggtg                                                                                                                           <  1:2083606/150‑1 (MQ=255)
                     ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggtg                                                                                                                           <  1:1070396/150‑1 (MQ=255)
                        ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGTTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTTGTGgta                                                                                                                        <  1:27372/150‑1 (MQ=255)
                        ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                        <  2:658381/150‑1 (MQ=255)
                           gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagta                                                                                                                     >  1:1244952/1‑150 (MQ=255)
                           gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagta                                                                                                                     >  2:526623/1‑150 (MQ=255)
                           gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagta                                                                                                                     >  2:1446882/1‑150 (MQ=255)
                           gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGGGTTATTGACGCCATGCGCGAGCT‑ACAAAAATGGGCTGTGCGCTGGTGgtcgta                                                                                                                     >  2:1547247/1‑150 (MQ=255)
                                gCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaa                                                                                                                <  1:589568/150‑1 (MQ=255)
                                gCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaa                                                                                                                <  1:2212949/150‑1 (MQ=255)
                                gCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaa                                                                                                                <  1:2516762/150‑1 (MQ=255)
                                   tACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggtgg                                                                                                                          >  2:1638476/1‑137 (MQ=255)
                                    aCCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAg                                                                                                            >  1:1452503/1‑150 (MQ=255)
                                     cccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGt                                                                                                           <  1:741198/150‑1 (MQ=255)
                                     cccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGt                                                                                                           <  1:796762/150‑1 (MQ=255)
                                     cccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGt                                                                                                           <  1:369293/150‑1 (MQ=255)
                                       cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                         >  2:1413925/1‑150 (MQ=255)
                                        gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTg                                                                                                        <  1:132995/150‑1 (MQ=255)
                                        gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTg                                                                                                        <  1:2579512/150‑1 (MQ=255)
                                         cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTgg                                                                                                       <  1:258890/150‑1 (MQ=255)
                                         cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGACTTTGCGCTGGTGGTAGTAACCAACCAGTCTgg                                                                                                       <  1:891422/150‑1 (MQ=255)
                                           aTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCa                                                                                                     <  2:2195975/150‑1 (MQ=255)
                                           aTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCa                                                                                                     <  2:2025688/150‑1 (MQ=255)
                                             ttttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                   >  2:2361821/1‑150 (MQ=255)
                                             ttttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                   >  1:2731683/1‑150 (MQ=255)
                                                  ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcg                                                                                              >  1:327406/1‑150 (MQ=255)
                                                  ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcg                                                                                              >  2:1091758/1‑150 (MQ=255)
                                                     aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCgg                                                                                           >  1:1001506/1‑150 (MQ=255)
                                                     aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCgg                                                                                           >  1:1880472/1‑150 (MQ=255)
                                                     aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCgg                                                                                           >  1:932737/1‑150 (MQ=255)
                                                     aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGGGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCACgg                                                                                           >  1:864901/1‑150 (MQ=255)
                                                        gTGATGGCACCATTAATGTCGATCTCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaa                                                                                        <  2:2431062/150‑1 (MQ=255)
                                                        gTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaa                                                                                        <  2:1136733/150‑1 (MQ=255)
                                                        gTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaa                                                                                        <  2:2481460/150‑1 (MQ=255)
                                                        gTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaa                                                                                        <  2:2126486/150‑1 (MQ=255)
                                                        gTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaa                                                                                        <  2:57027/150‑1 (MQ=255)
                                                        gTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaa                                                                                        <  2:454990/150‑1 (MQ=255)
                                                            tGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAAttt                                                                                    >  2:2529975/1‑150 (MQ=255)
                                                              gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                  >  1:223256/1‑150 (MQ=255)
                                                              gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                  <  2:2016907/150‑1 (MQ=255)
                                                                aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCg                                                                                >  1:884708/1‑150 (MQ=255)
                                                                aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCg                                                                                >  1:16548/1‑150 (MQ=255)
                                                                aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCg                                                                                >  1:2427962/1‑150 (MQ=255)
                                                                  cATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                              <  2:314922/150‑1 (MQ=255)
                                                                  cATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                              <  2:171599/150‑1 (MQ=255)
                                                                         gTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGt                                                                       >  2:107819/1‑150 (MQ=255)
                                                                               cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                                 <  1:1343604/150‑1 (MQ=255)
                                                                               cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                                 <  1:1913262/150‑1 (MQ=255)
                                                                               cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                                 <  1:293386/150‑1 (MQ=255)
                                                                                aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑CAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAAc                                                                >  2:2664239/1‑150 (MQ=255)
                                                                                aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAAc                                                                >  2:2567828/1‑150 (MQ=255)
                                                                                      aTGTCCATGAGATCGACAACTTTTAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAc                                                          <  2:1828815/150‑1 (MQ=255)
                                                                                      aTGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAc                                                          <  1:1539885/150‑1 (MQ=255)
                                                                                             tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑‑AAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgga                                                  <  2:1979968/150‑1 (MQ=255)
                                                                                                gaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatg                                                <  2:282005/150‑1 (MQ=255)
                                                                                                 aTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgg                                               >  2:1688741/1‑150 (MQ=255)
                                                                                                 aTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgg                                               >  1:1742061/1‑150 (MQ=255)
                                                                                                 aTCGACAACTTTGAATTTATCGACGGTGTTAATGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACAGAGtggatgg                                               >  2:1754552/1‑150 (MQ=255)
                                                                                                  tCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgga                                              <  2:2306808/150‑1 (MQ=255)
                                                                                                    gACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACt                                            >  2:250772/1‑150 (MQ=255)
                                                                                                     aCAAGTTTGAATTGCTCGACGGTGTTCTTGACGCCATGCGCGAGCT‑AAAAAATGGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACCGTTTGACACGCTGACCGAGTGGAGGGACTg                                           <  2:1465829/150‑1 (MQ=255)
                                                                                                     aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTg                                           <  1:2549886/150‑1 (MQ=255)
                                                                                                     aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTATGAAACGCTGACCGAGTGGATGGACTg                                           >  2:2105099/1‑150 (MQ=255)
                                                                                                        aCTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTc                                        >  1:1664391/1‑150 (MQ=255)
                                                                                                          tttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGATCGc                                      >  1:1789183/1‑150 (MQ=255)
                                                                                                                tttATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCgg                                <  2:2198909/150‑1 (MQ=255)
                                                                                                                tttATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCgg                                <  2:2252546/150‑1 (MQ=255)
                                                                                                                   aTCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAcc                             >  1:2281673/1‑150 (MQ=255)
                                                                                                                   aTCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAcc                             >  1:846032/1‑150 (MQ=255)
                                                                                                                   aTCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAcc                             >  1:577203/1‑150 (MQ=255)
                                                                                                                     cGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCga                           <  1:180837/150‑1 (MQ=255)
                                                                                                                     cGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCga                           <  1:1300010/150‑1 (MQ=255)
                                                                                                                     cGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCga                           <  1:1157679/150‑1 (MQ=255)
                                                                                                                     cGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCga                           <  2:2359299/150‑1 (MQ=255)
                                                                                                                     cGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCga                           <  2:2176099/150‑1 (MQ=255)
                                                                                                                          gtgtTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGt                      <  2:659782/150‑1 (MQ=255)
                                                                                                                           tgtTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTc                     <  1:2322168/150‑1 (MQ=255)
                                                                                                                           tgtTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTc                     <  1:767190/150‑1 (MQ=255)
                                                                                                                           tgtTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTc                     <  1:2322195/150‑1 (MQ=255)
                                                                                                                              tATTGACGCGATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTATATTTTCCGAAGCTCAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt                  <  2:2331511/150‑1 (MQ=255)
                                                                                                                              tATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt                  <  2:520672/150‑1 (MQ=255)
                                                                                                                              tATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt                  <  1:506223/150‑1 (MQ=255)
                                                                                                                              tATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt                  <  2:825902/150‑1 (MQ=255)
                                                                                                                              tATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCACGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt                  <  2:864901/150‑1 (MQ=255)
                                                                                                                                    cGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt            >  2:1368884/1‑150 (MQ=255)
                                                                                                                                      ccATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATgg          >  2:368146/1‑150 (MQ=255)
                                                                                                                                        aTGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTa        <  1:1413925/150‑1 (MQ=255)
                                                                                                                                              gAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCtatt  >  2:436381/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTATT  >  NC_000913/222811‑223103

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