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.100.0% 1469.8 / NA 294coding (126/576 nt)gmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (130/164);  total (130/164)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 323 total aligned reads displayed.

TGCCTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTAT  >  NC_000913/222809‑223102
                                                                                                                                                     |                                                                                                                                                
 gCCTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaa                                                                                                                                                 >  2:5090366/1‑150 (MQ=255)
  ccTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaa                                                                                                                                                >  1:444273/1‑150 (MQ=255)
  ccTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaa                                                                                                                                                >  2:4290855/1‑150 (MQ=255)
    ttattaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaa                                                                                                                                              >  2:4583385/1‑150 (MQ=255)
    ttattaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaa                                                                                                                                              >  1:5739915/1‑150 (MQ=255)
       ttaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaat                                                                                                                                           <  2:3740072/150‑2 (MQ=255)
          aGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatggg                                                                                                                                        <  1:3353058/150‑5 (MQ=255)
              gCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCttt                                                                                                                                   >  2:3451505/1‑150 (MQ=255)
              gCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCttt                                                                                                                                   <  1:2319560/150‑1 (MQ=255)
               cTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTg                                                                                                                                  <  1:3472643/150‑1 (MQ=255)
                        gTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggtgg                                                                                                                         <  1:2518060/150‑1 (MQ=255)
                          ggCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                       <  2:2099051/150‑1 (MQ=255)
                           gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                                      >  2:43453/1‑150 (MQ=255)
                           gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtag                                                                                                                      >  1:1540857/1‑150 (MQ=255)
                            cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                                     >  1:2217129/1‑150 (MQ=255)
                            cGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtagt                                                                                                                     >  2:806438/1‑150 (MQ=255)
                             gAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgtagta                                                                                                                    <  1:5425189/150‑1 (MQ=255)
                              aaGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaa                                                                                                                   >  2:219866/1‑150 (MQ=255)
                               agagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaac                                                                                                                  <  1:4823212/150‑1 (MQ=255)
                                    gTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaacc                                                                                                             >  2:57717/1‑150 (MQ=255)
                                     tACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAaccaacca                                                                                                            <  2:3834104/150‑1 (MQ=255)
                                        ccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTc                                                                                                         >  1:1864634/1‑150 (MQ=255)
                                         cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                        <  2:3910142/150‑1 (MQ=255)
                                         cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                        >  1:2664368/1‑150 (MQ=255)
                                          gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTg                                                                                                       >  2:2555283/1‑150 (MQ=255)
                                          gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGGCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTg                                                                                                       >  2:5716507/1‑150 (MQ=255)
                                           cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTgg                                                                                                      <  1:2003347/150‑1 (MQ=255)
                                            aaTTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGc                                                                                                     <  2:3393976/150‑1 (MQ=255)
                                               ttttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                  <  1:3559730/150‑1 (MQ=255)
                                                tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTg                                                                                                 >  2:75958/1‑150 (MQ=255)
                                                     tGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgc                                                                                            >  2:3276216/1‑150 (MQ=255)
                                                       aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCgg                                                                                          <  1:2947061/150‑1 (MQ=255)
                                                       aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGAGGTAGTAACCAACCAGTCTGGCATTGCTCGCgg                                                                                          >  2:3606411/1‑150 (MQ=255)
                                                        ccGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGt                                                                                         >  2:1111738/1‑150 (MQ=255)
                                                          gTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTaa                                                                                       <  1:1015066/150‑1 (MQ=255)
                                                                gCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTAc                                                                                 <  2:4971600/150‑1 (MQ=255)
                                                                  aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCg                                                                               >  2:552893/1‑150 (MQ=255)
                                                                  aCCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCg                                                                               >  1:3400904/1‑150 (MQ=255)
                                                                        aaTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcac                                                                         >  2:3208580/1‑150 (MQ=255)
                                                                            tCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGtt                                                                     <  2:1931854/150‑1 (MQ=255)
                                                                            tCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGtt                                                                     >  2:1304068/1‑150 (MQ=255)
                                                                             cGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGttt                                                                    <  1:612671/150‑1 (MQ=255)
                                                                              gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTg                                                                   >  1:4602033/1‑150 (MQ=255)
                                                                              gATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTg                                                                   >  1:1700191/1‑150 (MQ=255)
                                                                                tCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaa                                                                 <  1:3021315/150‑1 (MQ=255)
                                                                                tCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaa                                                                 >  2:5211397/1‑150 (MQ=255)
                                                                                 cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAATAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaaa                                                                <  2:2863769/150‑1 (MQ=255)
                                                                                  aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAAc                                                               >  2:1476597/1‑150 (MQ=255)
                                                                                  aCGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAAc                                                               <  2:4523376/150‑1 (MQ=255)
                                                                                     gCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCt                                                            <  2:2322733/150‑1 (MQ=255)
                                                                                       tATGTCCATGATATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGa                                                          >  2:660677/1‑150 (MQ=255)
                                                                                       tATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGa                                                          <  1:2200304/150‑1 (MQ=255)
                                                                                       tATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGa                                                          <  2:1386030/150‑1 (MQ=255)
                                                                                       tATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGa                                                          <  2:2783035/150‑1 (MQ=255)
                                                                                         tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                        <  2:4049986/150‑1 (MQ=255)
                                                                                         tGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGAcc                                                        <  2:3770296/150‑1 (MQ=255)
                                                                                           tCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                                      <  1:4927863/150‑1 (MQ=255)
                                                                                            ccATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAg                                                     <  1:476012/150‑1 (MQ=255)
                                                                                             cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGt                                                    <  2:1918986/150‑1 (MQ=255)
                                                                                             cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGt                                                    >  1:1183360/1‑150 (MQ=255)
                                                                                              aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtg                                                   >  2:2458425/1‑150 (MQ=255)
                                                                                               tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTCAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgg                                                  >  2:3138455/1‑150 (MQ=255)
                                                                                               tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtgg                                                  >  2:3869897/1‑150 (MQ=255)
                                                                                                 agaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggat                                                <  1:2465557/150‑1 (MQ=255)
                                                                                                 agaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggat                                                <  2:1989532/150‑1 (MQ=255)
                                                                                                  gaTCGACACCTTTGAATTTCTCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatg                                               <  2:5294301/150‑1 (MQ=255)
                                                                                                  gaTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatg                                               <  2:2487756/150‑1 (MQ=255)
                                                                                                   aTCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgg                                              <  2:2133184/150‑1 (MQ=255)
                                                                                                    tCGACAACTTTGTATTTATCGACGGTGTTCTTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTCACCAACCAGTCTGGCAGTGCTCGCGGTAAATTTACCGACGCACAGTTTGAAACGCTGACCGAGtggatgga                                             <  1:4068225/150‑1 (MQ=255)
                                                                                                     cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                            <  1:3848270/150‑1 (MQ=255)
                                                                                                       aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACCGTTTGAAACGCTGACCGAGTGGATGGACTg                                          >  2:3694762/1‑150 (MQ=255)
                                                                                                            tttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGc                                     >  2:902058/1‑150 (MQ=255)
                                                                                                            tttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGc                                     >  2:4761408/1‑150 (MQ=255)
                                                                                                            tttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGc                                     >  1:3952991/1‑150 (MQ=255)
                                                                                                            tttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGc                                     >  2:1559047/1‑150 (MQ=255)
                                                                                                             ttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCt                                    <  2:5413388/150‑1 (MQ=255)
                                                                                                               gAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTgg                                  <  2:3033672/150‑1 (MQ=255)
                                                                                                               gAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTgg                                  <  1:4266362/150‑1 (MQ=255)
                                                                                                               gAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTgg                                  <  1:125944/150‑1 (MQ=255)
                                                                                                                 aTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCg                                <  1:1932444/150‑1 (MQ=255)
                                                                                                                 aTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCg                                <  1:5536044/150‑1 (MQ=255)
                                                                                                                  tttATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCgg                               >  1:2836485/1‑150 (MQ=255)
                                                                                                                   ttATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGa                              <  1:5136683/150‑1 (MQ=255)
                                                                                                                   ttATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGa                              >  2:3189414/1‑150 (MQ=255)
                                                                                                                     aTCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAATATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAcc                            <  1:2641314/150‑1 (MQ=255)
                                                                                                                     aTCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATTGACTGGTCGCTGGCGGAcc                            <  1:5041613/150‑1 (MQ=255)
                                                                                                                     aTCGACGGTGTGATTGACGCCATGAGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCTAATCACCGTTGGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAcc                            <  2:3988559/150‑1 (MQ=255)
                                                                                                                      tCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCg                           >  1:800551/1‑150 (MQ=255)
                                                                                                                       cGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCga                          <  2:1049262/150‑1 (MQ=255)
                                                                                                                        gACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgag                         >  1:105280/1‑150 (MQ=255)
                                                                                                                           ggTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATg                      >  1:718256/1‑150 (MQ=255)
                                                                                                                               ttATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGa                  >  2:2679594/1‑150 (MQ=255)
                                                                                                                                 aTTGACTCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGTTAGTAACCAACCAGTCTGGCATTGCTCGCGGTTAATTTACCGAAGCTCAGTTTGAAACGCTGACCGAGTGGATTGACTGGTCGCTGGCGGACCGAGATGTCGATc                <  1:3656757/150‑1 (MQ=255)
                                                                                                                                  ttGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCt               >  1:2487271/1‑150 (MQ=255)
                                                                                                                                  ttGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCt               >  2:1789588/1‑150 (MQ=255)
                                                                                                                                   tGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTg              >  1:1877169/1‑150 (MQ=255)
                                                                                                                                      cGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt           <  2:468533/150‑1 (MQ=255)
                                                                                                                                        ccATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATgg         <  1:3147939/150‑1 (MQ=255)
                                                                                                                                             cgcgAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCt    >  1:3729301/1‑150 (MQ=255)
                                                                                                                                               cgAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCtat  <  1:4348490/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TGCCTTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTAT  >  NC_000913/222809‑223102

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: