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

TTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTATTATT  >  NC_000913/222813‑223106
                                                                                                                                                 |                                                                                                                                                    
ttattaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaa                                                                                                                                                   >  2:382768/1‑149 (MQ=255)
 tattaAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCAATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaa                                                                                                                                                  >  1:796615/1‑149 (MQ=255)
     aaGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatg                                                                                                                                              <  1:887403/149‑3 (MQ=255)
      aGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatgg                                                                                                                                             >  1:78692/1‑146 (MQ=255)
       ggAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatggg                                                                                                                                            <  1:457254/149‑5 (MQ=255)
       ggAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTaaaaaaatggg                                                                                                                                            <  1:378599/149‑5 (MQ=255)
        gAGCTATAAACGGTGGCGAAGAGCGTCCCCGCAATTTTTCTTGACCGTGAGGGCCCCATTAATGTCGATCCCGGCTATGTCCAGGAGATCGCCAACTTTGAATTTTTCGACGGGTTTTTTGCGGCCATGCGCGAGCTaaaaaaaaggcc                                                                                                                                           >  2:1124027/1‑145 (MQ=255)
            aataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGAGGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTg                                                                                                                                      <  2:887893/148‑1 (MQ=255)
                  aGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtgg                                                                                                                                <  1:177614/149‑1 (MQ=255)
                   ggTTGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggt                                                                                                                               <  2:206609/149‑1 (MQ=255)
                   ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGGACCATTAATGTCGATCACGGCTATGGCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggt                                                                                                                               <  2:642802/149‑1 (MQ=255)
                    gTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCTACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggtg                                                                                                                              <  1:236970/149‑1 (MQ=255)
                    gTGGCGAAGAGCGGACCCGCAATTTTTCTTGAGAGTGATGGCACCATTAATGTCAATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTTTTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCtggtg                                                                                                                              <  1:123487/149‑1 (MQ=255)
                       gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACAATTAATGTCGATCACGGATATTTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                           <  2:96407/149‑1 (MQ=255)
                            gagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaac                                                                                                                      <  1:402166/149‑1 (MQ=255)
                             agCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaacc                                                                                                                     >  1:923826/1‑149 (MQ=255)
                               cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaa                                                                                                                   <  2:1125495/149‑1 (MQ=255)
                               cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaa                                                                                                                   <  2:532763/149‑1 (MQ=255)
                                gTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                                  <  1:204678/149‑1 (MQ=255)
                                 tACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGAGCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaacc                                                                                                                 <  1:1161324/149‑1 (MQ=255)
                                  aCCCGCAATTTTTCTTGACCGTGATGGCACCATTACTGTCGATCACGGCTATGTCCATGAGCTCGACAACTATGAATTTATCGACGCTGTTATTGACGCCATCCGCGAGCT‑AAAAAACTGGGCTTTGCGCTGGTGGTAGTAaccaacca                                                                                                                <  2:133141/149‑1 (MQ=255)
                                      gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                            >  1:535801/1‑149 (MQ=255)
                                       cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGTTCGACAACTTTGATTTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTg                                                                                                           >  2:932396/1‑149 (MQ=255)
                                          tttttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCa                                                                                                        >  1:1031002/1‑149 (MQ=255)
                                            tttCTTGACCGTGATGGCACCTTTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                      <  2:800210/149‑1 (MQ=255)
                                                   aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGATTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgcg                                                                                               >  2:777317/1‑149 (MQ=255)
                                                    ccGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCgg                                                                                              >  1:635673/1‑149 (MQ=255)
                                                                   tAATGTCGATGACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGTGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCTAAGc                                                                               <  1:454519/149‑1 (MQ=255)
                                                                     aTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcac                                                                             >  1:32905/1‑149 (MQ=255)
                                                                      tGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcaca                                                                            <  2:950472/149‑1 (MQ=255)
                                                                             cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGATGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaa                                                                     <  1:66335/149‑1 (MQ=255)
                                                                                ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACg                                                                  <  2:1062509/149‑1 (MQ=255)
                                                                                   tATGTCCATGAGATCGACAACTTTGATTTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTg                                                               <  1:482334/149‑1 (MQ=255)
                                                                                      gTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCTAAGCACAATTTGAAAAGCTGAcc                                                            >  2:228528/1‑149 (MQ=255)
                                                                                      gTCCATGAGATCGACAACTTTGAATTTATCGACGGGGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTTCGCTGGTGGTATTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAACCGCTGAcc                                                            >  2:667490/1‑149 (MQ=255)
                                                                                       tCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCg                                                           >  2:278481/1‑149 (MQ=255)
                                                                                        ccATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTATATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                                          >  2:1175388/1‑149 (MQ=255)
                                                                                        ccATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                                          >  1:180744/1‑149 (MQ=255)
                                                                                         cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCCACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAg                                                         >  1:1070634/1‑149 (MQ=255)
                                                                                          aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGt                                                        <  1:164033/149‑1 (MQ=255)
                                                                                           tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAACTGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtg                                                       <  2:163758/149‑1 (MQ=255)
                                                                                              taTCGACAACTTTGAATTTATCGACGGTGTTATTGACGACATGCGCGAGCT‑AAAAAAATGGGCTTTGCGATGTTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTATATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggat                                                    <  1:757006/148‑1 (MQ=255)
                                                                                                tCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCAGACCGAGtggatgg                                                  <  1:753986/149‑1 (MQ=255)
                                                                                                 cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgga                                                 <  2:520893/149‑1 (MQ=255)
                                                                                                  gACATCTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                                <  1:204444/149‑1 (MQ=255)
                                                                                                  gACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                                >  1:345146/1‑149 (MQ=255)
                                                                                                   aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTAGAAACGCTGACCGAGACGATGGACt                                               >  1:1118037/1‑149 (MQ=255)
                                                                                                    aaaCTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCGCAGTTTGAAACGCTGACCGAGTGGATGGACTg                                              <  1:1018093/148‑1 (MQ=255)
                                                                                                         ttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGc                                         >  1:19202/1‑149 (MQ=255)
                                                                                                          tGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGGAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCt                                        >  2:1043477/1‑149 (MQ=255)
                                                                                                           gAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTg                                       >  1:845586/1‑149 (MQ=255)
                                                                                                            aaTTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGAAGTACCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTgg                                      <  1:694802/149‑1 (MQ=255)
                                                                                                             aTTTATCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGc                                     >  1:162906/1‑149 (MQ=255)
                                                                                                                  tCGACGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAcc                                <  2:923826/149‑1 (MQ=255)
                                                                                                                     aCGGTGTTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgag                             <  2:527306/149‑1 (MQ=255)
                                                                                                                        gtgtTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTAGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATATACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATCGACTGGTCGCTGGCGGACCGAGATg                          <  1:1052822/149‑1 (MQ=255)
                                                                                                                         tgtTATTGACGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGt                         <  1:932396/149‑1 (MQ=255)
                                                                                                                            tATTGACGCCATGCGCGAGCT‑AAAAACATGGGCTGTGCGCTCGTGGTTGTAACCAACCAGTCTGGCATTGCGCGCGGTAAATTTACCGAAGCACAGATTGAAACGCTGACCGAGTGGATGGACTGGGCGCTGGCGGACCGAGATGGCGa                      <  2:476312/149‑1 (MQ=255)
                                                                                                                               tGACGCCAGGCGCGAGCT‑AAAAAACGGGGCGTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAGGGCTCGCGGTACAATTACCGAAGCACAGCTTGAAACGCTGACCGAGGGGATGGACGGGTCGCTGGCGGACCGAGATGTCGATCt                   <  1:1059449/149‑1 (MQ=255)
                                                                                                                                 aCGCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTgg                 <  1:749413/149‑1 (MQ=255)
                                                                                                                                   gCCATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTAACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt               <  1:282897/149‑1 (MQ=255)
                                                                                                                                   cccATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt               <  2:938815/148‑1 (MQ=255)
                                                                                                                                    ccATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGACATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg              <  2:507901/149‑1 (MQ=255)
                                                                                                                                     cATGCGCGAGCT‑AATAAAATTGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATgg             <  1:163200/149‑1 (MQ=255)
                                                                                                                                     cATGCGCGAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATgg             >  1:1126323/1‑149 (MQ=255)
                                                                                                                                      aTGCGCGAGCT‑AAAAAAATGTGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGt            <  2:535801/149‑1 (MQ=255)
                                                                                                                                         cgcgAGCT‑AAAAACATGGGCTTTGCGCTGGTGCTAGAATCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCTTTCTGGATGGTATc         <  2:361655/149‑1 (MQ=255)
                                                                                                                                          gcgAGCT‑AAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCt        >  1:930943/1‑149 (MQ=255)
                                                                                                                                                 aaaaaaaaTGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGAACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTattatt  >  2:1000826/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTATTAAGGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCTATTATT  >  NC_000913/222813‑223106

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