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.97.6% 801.8 / 8.1 169coding (126/576 nt)gmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base . (80/85);  total (83/86)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
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 185 total aligned reads displayed.

CCTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTATT  >  NC_000913/222811‑223103
                                                                                                                                                   |                                                                                                                                                 
  ttattaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaa                                                                                                                                               >  2:1679548/1‑150 (MQ=255)
    attaAGGAGCTATACAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaa                                                                                                                                             <  2:3183807/150‑1 (MQ=255)
     ttaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaat                                                                                                                                            <  2:1994055/150‑2 (MQ=255)
     aaaaGCAGCAATAAAAGCTGCCGAAGAGCGTACCCCCAAATTTACTTGACCGTGATGGCACCAATAATGTCGATCACGGCTATGTCCATGAGATCGACAACATTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaat                                                                                                                                            <  2:1618499/148‑2 (MQ=255)
      taAGGAGCTCTACCAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatg                                                                                                                                           <  1:2289420/150‑3 (MQ=255)
      taAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatg                                                                                                                                           <  2:1372764/150‑3 (MQ=255)
      taAGGAGCTATAAAAGGGGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatg                                                                                                                                           <  2:2116668/150‑3 (MQ=255)
            gCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACACCTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑CAAAAAATGGGCttt                                                                                                                                    >  2:1987074/1‑150 (MQ=255)
            gCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCttt                                                                                                                                    <  1:444912/150‑1 (MQ=255)
            gCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCttt                                                                                                                                    >  1:2259036/1‑150 (MQ=255)
             cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTg                                                                                                                                   >  2:100336/1‑150 (MQ=255)
                    aGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggt                                                                                                                            <  2:919312/150‑1 (MQ=255)
                    aGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggt                                                                                                                            <  2:1022082/150‑1 (MQ=255)
                     ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggtg                                                                                                                           <  1:1027284/150‑1 (MQ=255)
                       tGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTggtggt                                                                                                                         >  1:1536272/1‑150 (MQ=255)
                       tGGCGAAGAGCGTACCCGCAATTGTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACGTTGAATTTATCGCCGGTGTTATTGACGCCATGCGGGAGCT‑GTATACATGGGCTTTGCGCTggtggt                                                                                                                         >  1:1471597/1‑150 (MQ=255)
                        ggCGACGCGCGTACCCGCAATTGTTCGTGACCGTGATGGCACCATTACTGTCGCTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑ACAACAATGGGCTTTGCGCTGGTGgta                                                                                                                        <  2:3082151/150‑1 (MQ=255)
                        ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                        <  2:870991/150‑1 (MQ=255)
                        ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                        <  1:3129252/150‑1 (MQ=255)
                         gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                                       >  1:2744375/1‑150 (MQ=255)
                          cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                                      >  2:2839412/1‑150 (MQ=255)
                          cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                                      <  2:2504014/150‑1 (MQ=255)
                           gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgtagta                                                                                                                     >  2:773503/1‑150 (MQ=255)
                            aaGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaa                                                                                                                    <  1:2263581/150‑1 (MQ=255)
                            aaGAGCGTACCCGCAATTTTTCTGGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaa                                                                                                                    <  1:319726/150‑1 (MQ=255)
                                 cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                               >  1:670602/1‑150 (MQ=255)
                                 cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                               >  2:1087254/1‑150 (MQ=255)
                                   tACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAaccaacca                                                                                                             <  1:2695947/150‑1 (MQ=255)
                                      ccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                                          >  2:2651821/1‑150 (MQ=255)
                                       cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                         <  1:2370927/150‑1 (MQ=255)
                                       cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                         <  2:484679/150‑1 (MQ=255)
                                       cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                         <  1:2049921/150‑1 (MQ=255)
                                         cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTgg                                                                                                       >  1:2777917/1‑150 (MQ=255)
                                         cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTgg                                                                                                       <  2:2491621/150‑1 (MQ=255)
                                          aaTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGc                                                                                                      <  1:100336/150‑1 (MQ=255)
                                          aaTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGc                                                                                                      <  1:1081599/150‑1 (MQ=255)
                                             ttttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                   <  1:2159467/150‑1 (MQ=255)
                                                  ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcg                                                                                              >  2:2419541/1‑150 (MQ=255)
                                                  ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcg                                                                                              >  2:1349188/1‑150 (MQ=255)
                                                    gACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgcg                                                                                            >  2:1786930/1‑150 (MQ=255)
                                                    gACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCAAGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgcg                                                                                            >  2:136462/1‑150 (MQ=255)
                                                     aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCgg                                                                                           >  1:1811654/1‑150 (MQ=255)
                                                         tGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaaa                                                                                       <  1:2003651/150‑1 (MQ=255)
                                                              gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                  <  1:1450378/150‑1 (MQ=255)
                                                              gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                  <  1:1208767/150‑1 (MQ=255)
                                                                aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCg                                                                                >  1:463683/1‑150 (MQ=255)
                                                                aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCg                                                                                >  1:2995932/1‑150 (MQ=255)
                                                                   aTTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAg                                                                             >  2:2269127/1‑150 (MQ=255)
                                                                      aaTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTAGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcac                                                                          >  1:794967/1‑150 (MQ=255)
                                                                           cGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGttt                                                                     <  1:3163740/150‑1 (MQ=255)
                                                                             aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                                   >  1:292615/1‑150 (MQ=255)
                                                                             aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                                   >  1:293722/1‑150 (MQ=255)
                                                                                aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAAc                                                                >  2:1449075/1‑150 (MQ=255)
                                                                                aCGGCTATGTCCATAAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAAc                                                                >  1:586283/1‑150 (MQ=255)
                                                                                  ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                              <  2:1783112/150‑1 (MQ=255)
                                                                                  ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                              >  2:3116288/1‑150 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         >  2:866535/1‑150 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         >  2:1926793/1‑150 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         >  2:1926720/1‑150 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         <  2:620228/150‑1 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         <  2:707295/150‑1 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGGGTTATTGAGGCCATGCGCGTGCT‑AACAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         <  2:1081521/150‑1 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTGTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         <  2:305770/150‑1 (MQ=255)
                                                                                           cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                              <  2:2283544/125‑1 (MQ=255)
                                                                                            aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtg                                                    <  2:2744375/150‑1 (MQ=255)
                                                                                             tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgg                                                   >  2:153092/1‑150 (MQ=255)
                                                                                             tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgg                                                   >  2:1813920/1‑150 (MQ=255)
                                                                                              gagaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgga                                                  <  2:2910958/150‑1 (MQ=255)
                                                                                               agaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggat                                                 >  2:674015/1‑150 (MQ=255)
                                                                                                   cGACAACATTGAATTTATCGACGGTGTTATTGACGACATGCTAGAGAT‑AAAAAAATGGGCTTTGAGCTGGTGGTATTAACCAACCAGTCTGGAATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAACGAGTGGATGGAc                                             <  2:2198805/150‑1 (MQ=255)
                                                                                                       aaCTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGt                                         <  1:2763623/150‑1 (MQ=255)
                                                                                                          tttGAATTTATCGACTGTGTTATTGACGCCATGAGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGc                                      <  1:2447500/150‑1 (MQ=255)
                                                                                                            tGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTg                                    >  1:443743/1‑150 (MQ=255)
                                                                                                            tGAATTTATCGACGGAGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGAAGTAACCCACCAGTCTGGCATTGCTCGCGGTAAATTTACAGAAGCACAGTTTTAAACGCTGACCGAGTGGATGGACTGGGCGCTg                                    >  2:3177212/1‑150 (MQ=255)
                                                                                                              aaTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTgg                                   >  2:476431/1‑150 (MQ=255)
                                                                                                              aaTTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGc                                  >  1:2527562/1‑150 (MQ=255)
                                                                                                               aTTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCg                                 <  2:1787708/150‑1 (MQ=255)
                                                                                                                tttATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCgg                                >  1:913861/1‑150 (MQ=255)
                                                                                                                tttATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCgg                                >  2:1161969/1‑150 (MQ=255)
                                                                                                                  tATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAc                              <  1:1423374/150‑1 (MQ=255)
                                                                                                                      gACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgag                          >  2:346595/1‑150 (MQ=255)
                                                                                                                      gACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgag                          >  2:449567/1‑150 (MQ=255)
                                                                                                                         ggTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATg                       <  1:1580601/150‑1 (MQ=255)
                                                                                                                             ttATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGa                   <  1:1934587/150‑1 (MQ=255)
                                                                                                                             ttATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGa                   >  2:1890061/1‑150 (MQ=255)
                                                                                                                              tATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt                  <  1:878808/150‑1 (MQ=255)
                                                                                                                               aTTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATc                 >  1:330694/1‑150 (MQ=255)
                                                                                                                                   aCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa             <  1:1578341/150‑1 (MQ=255)
                                                                                                                                   aCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa             <  1:1589800/150‑1 (MQ=255)
                                                                                                                                   aCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa             <  1:717302/150‑1 (MQ=255)
                                                                                                                                   aCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa             <  2:638065/150‑1 (MQ=255)
                                                                                                                                    cGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt            <  1:57852/150‑1 (MQ=255)
                                                                                                                                     gCCCTGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCGTTGCTCGCGGTAAATTTCCCGAAGCACAGTGTGAAACGCTGACCGCGAGGCTGGACTGGTCGCTGGCGGACCGCGATGGCGAGCTGGATg           >  2:2058249/1‑150 (MQ=255)
                                                                                                                                     gCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg           >  2:2855541/1‑150 (MQ=255)
                                                                                                                                     gCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg           >  1:161753/1‑150 (MQ=255)
                                                                                                                                     gCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg           >  2:1586463/1‑150 (MQ=255)
                                                                                                                                        aTGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGATAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTa        <  1:389109/150‑1 (MQ=255)
                                                                                                                                           cgcgAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCt     >  2:1776195/1‑150 (MQ=255)
                                                                                                                                              gAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCtatt  >  1:2471724/1‑150 (MQ=255)
                                                                                                                                              gAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCtatt  >  1:2690926/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTATT  >  NC_000913/222811‑223103

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