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)
            aataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGAGGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTg                                                                                                                                      <  2:887893/148‑1 (MQ=255)
                  aGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtgg                                                                                                                                <  1:177614/149‑1 (MQ=255)
                   ggTTGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggt                                                                                                                               <  2:206609/149‑1 (MQ=255)
                   ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGGACCATTAATGTCGATCACGGCTATGGCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggt                                                                                                                               <  2:642802/149‑1 (MQ=255)
                    gTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCTACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggtg                                                                                                                              <  1:236970/149‑1 (MQ=255)
                    gTGGCGAAGAGCGGACCCGCAATTTTTCTTGAGAGTGATGGCACCATTAATGTCAATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTTTTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCtggtg                                                                                                                              <  1:123487/149‑1 (MQ=255)
                       gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACAATTAATGTCGATCACGGATATTTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGgta                                                                                                                           <  2:96407/149‑1 (MQ=255)
                            gagCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaac                                                                                                                      <  1:402166/149‑1 (MQ=255)
                             agCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaacc                                                                                                                     >  1:923826/1‑149 (MQ=255)
                               cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaa                                                                                                                   <  2:1125495/149‑1 (MQ=255)
                               cGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaa                                                                                                                   <  2:532763/149‑1 (MQ=255)
                                gTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaac                                                                                                                  <  1:204678/149‑1 (MQ=255)
                                 tACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGAGCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTaaccaacc                                                                                                                 <  1:1161324/149‑1 (MQ=255)
                                  aCCCGCAATTTTTCTTGACCGTGATGGCACCATTACTGTCGATCACGGCTATGTCCATGAGCTCGACAACTATGAATTTATCGACGCTGTTATTGACGCCATCCGCGAGCTAAAAAACTGGGCTTTGCGCTGGTGGTAGTAaccaacca                                                                                                                <  2:133141/149‑1 (MQ=255)
                                      gCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCt                                                                                                            >  1:535801/1‑149 (MQ=255)
                                       cAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGTTCGACAACTTTGATTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTg                                                                                                           >  2:932396/1‑149 (MQ=255)
                                          tttttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCa                                                                                                        >  1:1031002/1‑149 (MQ=255)
                                            tttCTTGACCGTGATGGCACCTTTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAtt                                                                                                      <  2:800210/149‑1 (MQ=255)
                                                   aCCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGATTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTcgcg                                                                                               >  2:777317/1‑149 (MQ=255)
                                                    ccGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCgg                                                                                              >  1:635673/1‑149 (MQ=255)
                                                                   tAATGTCGATGACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGTGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCTAAGc                                                                               <  1:454519/149‑1 (MQ=255)
                                                                     aTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcac                                                                             >  1:32905/1‑149 (MQ=255)
                                                                      tGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGcaca                                                                            <  2:950472/149‑1 (MQ=255)
                                                                             cACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGATGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGaa                                                                     <  1:66335/149‑1 (MQ=255)
                                                                                ggCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACg                                                                  <  2:1062509/149‑1 (MQ=255)
                                                                                   tATGTCCATGAGATCGACAACTTTGATTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTg                                                               <  1:482334/149‑1 (MQ=255)
                                                                                      gTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCTAAGCACAATTTGAAAAGCTGAcc                                                            >  2:228528/1‑149 (MQ=255)
                                                                                      gTCCATGAGATCGACAACTTTGAATTTATCGACGGGGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTTCGCTGGTGGTATTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAACCGCTGAcc                                                            >  2:667490/1‑149 (MQ=255)
                                                                                       tCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCg                                                           >  2:278481/1‑149 (MQ=255)
                                                                                        ccATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTATATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                                          >  2:1175388/1‑149 (MQ=255)
                                                                                        ccATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGa                                                          >  1:180744/1‑149 (MQ=255)
                                                                                         cATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCCACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAg                                                         >  1:1070634/1‑149 (MQ=255)
                                                                                          aTGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGt                                                        <  1:164033/149‑1 (MQ=255)
                                                                                           tGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAACTGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtg                                                       <  2:163758/149‑1 (MQ=255)
                                                                                              taTCGACAACTTTGAATTTATCGACGGTGTTATTGACGACATGCGCGAGCTAAAAAAATGGGCTTTGCGATGTTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTATATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggat                                                    <  1:757006/148‑1 (MQ=255)
                                                                                                tCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCAGACCGAGtggatgg                                                  <  1:753986/149‑1 (MQ=255)
                                                                                                 cGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGtggatgga                                                 <  2:520893/149‑1 (MQ=255)
                                                                                                  gACATCTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                                <  1:204444/149‑1 (MQ=255)
                                                                                                  gACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGAc                                                >  1:345146/1‑149 (MQ=255)
                                                                                                   aCAACTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTAGAAACGCTGACCGAGACGATGGACt                                               >  1:1118037/1‑149 (MQ=255)
                                                                                                    aaaCTTTGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCGCAGTTTGAAACGCTGACCGAGTGGATGGACTg                                              <  1:1018093/148‑1 (MQ=255)
                                                                                                         ttGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGc                                         >  1:19202/1‑149 (MQ=255)
                                                                                                          tGAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGGAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCt                                        >  2:1043477/1‑149 (MQ=255)
                                                                                                           gAATTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTg                                       >  1:845586/1‑149 (MQ=255)
                                                                                                            aaTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGAAGTACCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTgg                                      <  1:694802/149‑1 (MQ=255)
                                                                                                             aTTTATCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGc                                     >  1:162906/1‑149 (MQ=255)
                                                                                                                  tCGACGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGAcc                                <  2:923826/149‑1 (MQ=255)
                                                                                                                     aCGGTGTTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCgag                             <  2:527306/149‑1 (MQ=255)
                                                                                                                        gtgtTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTAGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATATACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATCGACTGGTCGCTGGCGGACCGAGATg                          <  1:1052822/149‑1 (MQ=255)
                                                                                                                         tgtTATTGACGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGt                         <  1:932396/149‑1 (MQ=255)
                                                                                                                            tATTGACGCCATGCGCGAGCTAAAAACATGGGCTGTGCGCTCGTGGTTGTAACCAACCAGTCTGGCATTGCGCGCGGTAAATTTACCGAAGCACAGATTGAAACGCTGACCGAGTGGATGGACTGGGCGCTGGCGGACCGAGATGGCGa                      <  2:476312/149‑1 (MQ=255)
                                                                                                                               tGACGCCAGGCGCGAGCTAAAAAACGGGGCGTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCAGGGCTCGCGGTACAATTACCGAAGCACAGCTTGAAACGCTGACCGAGGGGATGGACGGGTCGCTGGCGGACCGAGATGTCGATCt                   <  1:1059449/149‑1 (MQ=255)
                                                                                                                                 aCGCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTgg                 <  1:749413/149‑1 (MQ=255)
                                                                                                                                   gCCATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTAACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt               <  1:282897/149‑1 (MQ=255)
                                                                                                                                   cccATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGAt               <  2:938815/148‑1 (MQ=255)
                                                                                                                                    ccATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGACATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATg              <  2:507901/149‑1 (MQ=255)
                                                                                                                                     cATGCGCGAGCTAATAAAATTGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATgg             <  1:163200/149‑1 (MQ=255)
                                                                                                                                     cATGCGCGAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATgg             >  1:1126323/1‑149 (MQ=255)
                                                                                                                                      aTGCGCGAGCTAAAAAAATGTGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGt            <  2:535801/149‑1 (MQ=255)
                                                                                                                                         cgcgAGCTAAAAACATGGGCTTTGCGCTGGTGCTAGAATCCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCTTTCTGGATGGTATc         <  2:361655/149‑1 (MQ=255)
                                                                                                                                          gcgAGCTAAAAAAATGGGCTTTGCGCTGGTGGTAGTAACCAACCAGTCTGGCATTGCTCGCGGTAAATTTACCGAAGCACAGTTTGAAACGCTGACCGAGTGGATGGACTGGTCGCTGGCGGACCGAGATGTCGATCTGGATGGTATCt        >  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: