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)
            gCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACACCTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTCAAAAAATGGGCttt                                                                                                                                    >  2:1987074/1‑150 (MQ=255)
            gCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCttt                                                                                                                                    <  1:444912/150‑1 (MQ=255)
            gCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCttt                                                                                                                                    >  1:2259036/1‑150 (MQ=255)
             cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTg                                                                                                                                   >  2:100336/1‑150 (MQ=255)
                    aGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggt                                                                                                                            <  2:919312/150‑1 (MQ=255)
                    aGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggt                                                                                                                            <  2:1022082/150‑1 (MQ=255)
                     ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggtg                                                                                                                           <  1:1027284/150‑1 (MQ=255)
                       tGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTggtggt                                                                                                                         >  1:1536272/1‑150 (MQ=255)
                       tGGCGAAGAGCGTACCCGCAATTGTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACGTTGAATTTATCGCCGGTGTTATTGACGCCATGCGGGAGCTGTATACATGGGCTTTGCGCTggtggt                                                                                                                         >  1:1471597/1‑150 (MQ=255)
                        ggCGACGCGCGTACCCGCAATTGTTCGTGACCGTGATGGCACCATTACTGTCGCTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTACAACAATGGGCTTTGCGCTGGTGgta                                                                                                                        <  2:3082151/150‑1 (MQ=255)
                        ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                        <  2:870991/150‑1 (MQ=255)
                        ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                        <  1:3129252/150‑1 (MQ=255)
                         gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                                       >  1:2744375/1‑150 (MQ=255)
                          cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                                      >  2:2839412/1‑150 (MQ=255)
                          cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                                      <  2:2504014/150‑1 (MQ=255)
                           gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtagta                                                                                                                     >  2:773503/1‑150 (MQ=255)
                            aaGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaa                                                                                                                    <  1:2263581/150‑1 (MQ=255)
                            aaGAGCGTACCCGCAATTTTTCTGGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaa                                                                                                                    <  1:319726/150‑1 (MQ=255)
                                 cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                               >  1:670602/1‑150 (MQ=255)
                                 cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                               >  2:1087254/1‑150 (MQ=255)
                                   tACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAaccaacca                                                                                                             <  1:2695947/150‑1 (MQ=255)
                                      ccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                                          >  2:2651821/1‑150 (MQ=255)
                                       cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                         <  1:2370927/150‑1 (MQ=255)
                                       cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                         <  2:484679/150‑1 (MQ=255)
                                       cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                         <  1:2049921/150‑1 (MQ=255)
                                         cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTgg                                                                                                       >  1:2777917/1‑150 (MQ=255)
                                         cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTgg                                                                                                       <  2:2491621/150‑1 (MQ=255)
                                          aaTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGc                                                                                                      <  1:100336/150‑1 (MQ=255)
                                          aaTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGc                                                                                                      <  1:1081599/150‑1 (MQ=255)
                                             ttttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                   <  1:2159467/150‑1 (MQ=255)
                                                  ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcg                                                                                              >  2:2419541/1‑150 (MQ=255)
                                                  ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcg                                                                                              >  2:1349188/1‑150 (MQ=255)
                                                    gACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgcg                                                                                            >  2:1786930/1‑150 (MQ=255)
                                                    gACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCAAGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgcg                                                                                            >  2:136462/1‑150 (MQ=255)
                                                     aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCgg                                                                                           >  1:1811654/1‑150 (MQ=255)
                                                         tGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaaa                                                                                       <  1:2003651/150‑1 (MQ=255)
                                                              gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                  <  1:1450378/150‑1 (MQ=255)
                                                              gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                  <  1:1208767/150‑1 (MQ=255)
                                                                aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCg                                                                                >  1:463683/1‑150 (MQ=255)
                                                                aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCg                                                                                >  1:2995932/1‑150 (MQ=255)
                                                                   aTTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAg                                                                             >  2:2269127/1‑150 (MQ=255)
                                                                      aaTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTAGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcac                                                                          >  1:794967/1‑150 (MQ=255)
                                                                           cGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGttt                                                                     <  1:3163740/150‑1 (MQ=255)
                                                                             aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                                   >  1:292615/1‑150 (MQ=255)
                                                                             aTCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGa                                                                   >  1:293722/1‑150 (MQ=255)
                                                                                aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAAc                                                                >  2:1449075/1‑150 (MQ=255)
                                                                                aCGGCTATGTCCATAAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAAc                                                                >  1:586283/1‑150 (MQ=255)
                                                                                  ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                              <  2:1783112/150‑1 (MQ=255)
                                                                                  ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGc                                                              >  2:3116288/1‑150 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         >  2:866535/1‑150 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         >  2:1926793/1‑150 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         >  2:1926720/1‑150 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         <  2:620228/150‑1 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         <  2:707295/150‑1 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTTTGAATTTATCGACGGGGTTATTGAGGCCATGCGCGTGCTAACAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         <  2:1081521/150‑1 (MQ=255)
                                                                                       tGTCCATGAGATCGACAACTGTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                         <  2:305770/150‑1 (MQ=255)
                                                                                           cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGaa                                                                              <  2:2283544/125‑1 (MQ=255)
                                                                                            aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtg                                                    <  2:2744375/150‑1 (MQ=255)
                                                                                             tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgg                                                   >  2:153092/1‑150 (MQ=255)
                                                                                             tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgg                                                   >  2:1813920/1‑150 (MQ=255)
                                                                                              gagaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgga                                                  <  2:2910958/150‑1 (MQ=255)
                                                                                               agaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggat                                                 >  2:674015/1‑150 (MQ=255)
                                                                                                   cGACAACATTGAATTTATCGACGGTGTTATTGACGACATGCTAGAGATAAAAAAATGGGCTTTGAGCTGGTGGTATTAACCAACCAGTCTGGAATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAACGAGTGGATGGAc                                             <  2:2198805/150‑1 (MQ=255)
                                                                                                       aaCTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGt                                         <  1:2763623/150‑1 (MQ=255)
                                                                                                          tttGAATTTATCGACTGTGTTATTGACGCCATGAGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGc                                      <  1:2447500/150‑1 (MQ=255)
                                                                                                            tGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTg                                    >  1:443743/1‑150 (MQ=255)
                                                                                                            tGAATTTATCGACGGAGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGAAGTAACCCACCAGTCTGGCATTGCTCGCGGTAAATTTACAGAAGCACAGTTTTAAACGCTGACCGAGTGGATGGACTGGGCGCTg                                    >  2:3177212/1‑150 (MQ=255)
                                                                                                              aaTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTgg                                   >  2:476431/1‑150 (MQ=255)
                                                                                                              aaTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGc                                  >  1:2527562/1‑150 (MQ=255)
                                                                                                               aTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCg                                 <  2:1787708/150‑1 (MQ=255)
                                                                                                                tttATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCgg                                >  1:913861/1‑150 (MQ=255)
                                                                                                                tttATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCgg                                >  2:1161969/1‑150 (MQ=255)
                                                                                                                  tATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAc                              <  1:1423374/150‑1 (MQ=255)
                                                                                                                      gACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgag                          >  2:346595/1‑150 (MQ=255)
                                                                                                                      gACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgag                          >  2:449567/1‑150 (MQ=255)
                                                                                                                         ggTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATg                       <  1:1580601/150‑1 (MQ=255)
                                                                                                                             ttATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGa                   <  1:1934587/150‑1 (MQ=255)
                                                                                                                             ttATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGa                   >  2:1890061/1‑150 (MQ=255)
                                                                                                                              tATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGAt                  <  1:878808/150‑1 (MQ=255)
                                                                                                                               aTTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATc                 >  1:330694/1‑150 (MQ=255)
                                                                                                                                   aCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa             <  1:1578341/150‑1 (MQ=255)
                                                                                                                                   aCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa             <  1:1589800/150‑1 (MQ=255)
                                                                                                                                   aCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa             <  1:717302/150‑1 (MQ=255)
                                                                                                                                   aCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGa             <  2:638065/150‑1 (MQ=255)
                                                                                                                                    cGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt            <  1:57852/150‑1 (MQ=255)
                                                                                                                                     gCCCTGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCGTTGCTCGCGGTAAATTTCCCGAAGCACAGTGTGAAACGCTGACCGCGAGGCTGGACTGGTCGCTGGCGGACCGCGATGGCGAGCTGGATg           >  2:2058249/1‑150 (MQ=255)
                                                                                                                                     gCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg           >  2:2855541/1‑150 (MQ=255)
                                                                                                                                     gCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg           >  1:161753/1‑150 (MQ=255)
                                                                                                                                     gCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg           >  2:1586463/1‑150 (MQ=255)
                                                                                                                                        aTGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGATAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTa        <  1:389109/150‑1 (MQ=255)
                                                                                                                                           cgcgAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCt     >  2:1776195/1‑150 (MQ=255)
                                                                                                                                              gAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCtatt  >  1:2471724/1‑150 (MQ=255)
                                                                                                                                              gAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCtatt  >  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: