Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 223,342 T→G D170E (GAT→GAG 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_000913223,3420TG100.0% 232.3 / NA 65D170E (GAT→GAGgmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (35/30);  total (35/30)

GATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCG  >  NC_000913/223196‑223488
                                                                                                                                                  |                                                                                                                                                  
gATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATACATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGAAACAAAAGTGATGGTGCGTACGGGGAAACCTATTACGCCTGCAGCAGAAAACGCGGCGGAGTggg                                                                                                                                                >  2:91866/1‑151 (MQ=255)
     tttGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGtt                                                                                                                                            >  1:604730/1‑150 (MQ=255)
      ttGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTa                                                                                                                                           >  2:100308/1‑150 (MQ=255)
      ttGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTa                                                                                                                                           >  1:90716/1‑150 (MQ=255)
          atatTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATa                                                                                                                                       <  2:113544/150‑1 (MQ=255)
          atatTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAg                                                                                                                                      <  1:362233/151‑1 (MQ=255)
                atatGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGc                                                                                                                                >  2:87195/1‑151 (MQ=255)
                atatGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGc                                                                                                                                >  2:681676/1‑151 (MQ=255)
                 tatGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGc                                                                                                                                >  1:841577/1‑150 (MQ=255)
                  atGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGc                                                                                                                                >  2:386632/1‑149 (MQ=255)
                        gCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGc                                                                                                                        >  1:700023/1‑151 (MQ=255)
                        gCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGc                                                                                                                        >  1:81499/1‑151 (MQ=255)
                          ttcttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCg                                                                                                                      >  1:274877/1‑151 (MQ=255)
                            cttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGc                                                                                                                     >  2:486233/1‑150 (MQ=255)
                               atatGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGc                                                                                                                 <  1:31682/151‑1 (MQ=255)
                                    gtggGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATa                                                                                                             >  2:810385/1‑150 (MQ=255)
                                      gggCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATaaa                                                                                                           >  2:4625/1‑150 (MQ=255)
                                        gCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATaaaaaa                                                                                                        <  1:440529/151‑1 (MQ=255)
                                        gCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATaaaaa                                                                                                         >  1:448305/1‑150 (MQ=255)
                                            tAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAagc                                                                                                      >  2:217965/1‑149 (MQ=255)
                                             aaaTTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAagcagc                                                                                                   >  2:609625/1‑151 (MQ=255)
                                             aaaTTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAagcagc                                                                                                   <  1:684102/151‑1 (MQ=255)
                                                 tAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCaaa                                                                                                <  2:215920/150‑1 (MQ=255)
                                                 tAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGAAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCaa                                                                                                 >  1:793585/1‑149 (MQ=255)
                                                           gcaggcagCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGcaca                                                                                      >  2:386158/1‑150 (MQ=255)
                                                                cagCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGAtt                                                                                >  2:884387/1‑151 (MQ=255)
                                                                 agCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGAtt                                                                                <  1:810403/150‑1 (MQ=255)
                                                                    ggTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTaaaa                                                                            <  2:527340/151‑1 (MQ=255)
                                                                       tGCGGCGAACGTGGGAACAAAAGCGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAAAGATTaaa                                                                             >  1:362620/1‑147 (MQ=255)
                                                                          ggcgAACGTGGGAACAAAAGTGTTGGTGCGTACGGTTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAg                                                                      <  1:208384/151‑1 (MQ=255)
                                                                          ggcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATg                                                                        >  1:472338/1‑149 (MQ=255)
                                                                          ggcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGa                                                                       <  1:805558/150‑1 (MQ=255)
                                                                           gcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGc                                                                     <  2:212/151‑1 (MQ=255)
                                                                             gAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCg                                                                    <  1:535306/150‑1 (MQ=255)
                                                                                 gTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGa                                                               <  1:71045/151‑1 (MQ=255)
                                                                                  tGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGa                                                               >  2:153643/1‑150 (MQ=255)
                                                                                  tGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGa                                                               >  1:724402/1‑150 (MQ=255)
                                                                                   gggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGaa                                                              <  2:152041/150‑1 (MQ=255)
                                                                                   gggAACAAAAGTGCTGGGGCGTACGGGTAAACCTACTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGaa                                                              <  2:458415/150‑1 (MQ=255)
                                                                                    ggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAAt                                                            >  1:811883/1‑151 (MQ=255)
                                                                                     gAACAAAAGTGCTGGTGCGTACGGGTAAACCGATTACGCCTGAAGCAGAAAACGCGGGGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAAt                                                            <  1:706119/150‑1 (MQ=255)
                                                                                         aaaagggagggggCGTAGGGGTAAACCTATAAGGCGTGAAGAAGAAAACGCGGGGAAGTGGGTGTTAAGTAGCCAGGCAAACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATaaaaa                                                       <  2:22390/139‑1 (MQ=255)
                                                                                                 tGGTGCGTGCGGTTAAACCTATTACACGTTAAACAAAAAACGCGGGCGAGTGGTTGTTAAATGACCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCAt                                                  <  2:693268/148‑1 (MQ=255)
                                                                                                  ggTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCAtttt                                               <  1:645576/150‑1 (MQ=255)
                                                                                                  ggTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCAtttt                                               <  2:699029/150‑1 (MQ=255)
                                                                                                     gCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTcc                                            <  1:533978/150‑1 (MQ=255)
                                                                                                     gCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTcc                                            <  2:644164/150‑1 (MQ=255)
                                                                                                         aCGGGTAACCCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCtt                                        <  2:654600/150‑1 (MQ=255)
                                                                                                            ggTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCt                                    >  2:294500/1‑151 (MQ=255)
                                                                                                               aaaCCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTc                                  <  2:317513/150‑1 (MQ=255)
                                                                                                               aaaCCTATTACGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTc                                  <  2:90715/150‑1 (MQ=255)
                                                                                                                        aCGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGAc                        >  1:646014/1‑151 (MQ=255)
                                                                                                                         cGCCTGAAGCAGAAACCGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACt                       <  1:885930/151‑1 (MQ=255)
                                                                                                                          gCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTc                      >  2:150925/1‑151 (MQ=255)
                                                                                                                             tGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCt                   >  2:818362/1‑151 (MQ=255)
                                                                                                                               aaGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCtat                 >  2:73565/1‑151 (MQ=255)
                                                                                                                                  cAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATaa               >  1:596452/1‑150 (MQ=255)
                                                                                                                                     aaaaCGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATgcg           >  2:756532/1‑151 (MQ=255)
                                                                                                                                     aaaaCGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATgc            <  1:908495/150‑1 (MQ=255)
                                                                                                                                         cgcgGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCt        <  2:763007/150‑1 (MQ=255)
                                                                                                                                         cgcgGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCt        <  1:863078/150‑1 (MQ=255)
                                                                                                                                         cgcgGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTc       >  2:414084/1‑151 (MQ=255)
                                                                                                                                          gcggcgGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTcc      >  1:917607/1‑151 (MQ=255)
                                                                                                                                            ggcggAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCa     <  1:578491/150‑1 (MQ=255)
                                                                                                                                               ggAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCg  >  2:52361/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCCGACTCCCTATAATGCGCCTCCATCG  >  NC_000913/223196‑223488

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