Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 4,296,061 C→T 16.6% intergenic (+266/+376) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,296,0610CT16.6% 142.1 / 21.5 88intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (33/38);  new base T (8/7);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 7.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCC  >  NC_000913_3_bme_pgi/4295913‑4296208
                                                                                                                                                    |                                                                                                                                                   
gtgtgttAAGGCGTTTATGCCGCATCAGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAc                                                                                                                                                     <  1:381814/143‑1 (MQ=12)
    gATAAGTCGTTTATGCCGTGTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                                 <  1:179333/149‑1 (MQ=255)
    gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                                 <  1:214699/149‑1 (MQ=17)
      tAAGGCGTTTATGCCGCGTTCGATATCAATGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                               <  1:154691/149‑1 (MQ=255)
        aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAta                                                                                                                                             >  1:586409/1‑148 (MQ=17)
          gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCg                                                                                                                                           >  2:404299/1‑149 (MQ=17)
            gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATACCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCCCATACGACATCAAACCCtgctgcacc                                                                                                                                         >  2:423347/1‑146 (MQ=11)
             tttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTACGTCGGATAAGGCGTTTACGCCGCATCCGACATCCATACCTGAGGCGACg                                                                                                                                        >  2:462795/1‑149 (MQ=11)
                    cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCATACGCAGGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGCCGGATAAGCAGTTTACGCCGCATCCGACATGCATGCCTGATGCGACGCTTGcct                                                                                                                                 >  2:146587/1‑148 (MQ=11)
                       aTCCGACATCAACGCCCGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                              <  2:586595/149‑1 (MQ=38)
                        tacGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAACAGAAAAAAACAAAACAGCGCAATAGCCTGATTTCGCGAGATTTTGTAGGTCAAATAAGGCATTTAAGCCGAATTGCACATCAATGCCTGATGCGTCGCTTTCcgcgac                                                                                                                             >  2:545446/3‑147 (MQ=2)
                        tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                             >  1:756081/1‑149 (MQ=21)
                        tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                             >  2:841919/1‑149 (MQ=38)
                        tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                             >  1:93358/1‑149 (MQ=37)
                        tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCCCATCCTACATCAACCCCTGATCCGACGCTTGCcgccac                                                                                                                             >  1:324431/1‑146 (MQ=11)
                          cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCGCCGACATCAATGCCTGATCAGAGCCTTGCCGCGtttt                                                                                                                           >  2:441629/1‑149 (MQ=11)
                          cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACg                                                                                                                                                    >  1:513795/1‑124 (MQ=32)
                           gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGGCAGCGCAATAGCCTGATTTAGCGTGACTTTGTAGGTCGGATAAGGCGTTTACGCAGCATCCGACATAAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                          >  1:834506/1‑149 (MQ=11)
                            aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                         >  2:628668/1‑149 (MQ=25)
                            aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGACTTACCGTGATTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAACGCCTGAGACGCAGCTTGCCGCGTCTTAt                                                                                                                         >  2:297533/1‑149 (MQ=11)
                            aCATCAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                         <  1:524404/149‑1 (MQ=255)
                               tgaaCGCCTGATGCGACGCTTAACGCGTTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                      <  1:422021/147‑1 (MQ=255)
                               tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCTTACGCTAGACAGCGTAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                      <  1:242970/149‑1 (MQ=255)
                               tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                      <  1:733926/149‑1 (MQ=39)
                                 aatgCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTTGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                    <  1:123335/146‑1 (MQ=0)
                                 aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                    >  2:846817/1‑149 (MQ=37)
                                   cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATg                                                                                                                                             >  1:471655/1‑122 (MQ=32)
                                      cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                       >  2:268386/1‑141 (MQ=37)
                                      cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                       <  1:268386/141‑1 (MQ=39)
                                      cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGGGAATTTGTAGGTTGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCCACGCTTGCCGCGTCTTATCAGCCCTAcc                                                                                                               >  1:886896/1‑148 (MQ=11)
                                       tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTTCCGCGTCTTTTCAGGCCTAccc                                                                                                              >  2:280345/1‑147 (MQ=19)
                                        gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGAGTCCGAGAGCAAAGCCTGATGCGGCGCTTGCCACGTATGATCAGGCATAcccc                                                                                                             >  1:568678/1‑149 (MQ=11)
                                         aTGCGCCGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                            <  1:317234/149‑1 (MQ=255)
                                         aTGCGACGTTTAACGCGTCTTATCAGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctta                                                                                                            <  1:121873/149‑6 (MQ=11)
                                           gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                              <  2:749354/129‑1 (MQ=25)
                                           gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                              >  1:749354/1‑129 (MQ=38)
                                           gCGAAGCTTAACGCGTCTTATGAGTCTTAAGCCAGACGGCGCAATGGTATGATTAAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                          <  1:392466/149‑1 (MQ=11)
                                              aCGCTTAACGCGTCTTATCAGTCCTACGCCAGACAGCGTAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAg                                                                                                       <  1:158117/149‑1 (MQ=255)
                                                gCTTACGGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                                                                                                     <  1:267277/149‑1 (MQ=255)
                                                gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATGCTAGACAgtg                                                                                                     >  1:474605/1‑147 (MQ=255)
                                                  ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAgcgga                                                                                                   >  1:754449/1‑147 (MQ=255)
                                                    accgcgTCTTATCAGGCCTAGGCCAGAGGGCGCAATAGCCTTATTGAGCGTCATTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGTCAGCGCAAt                                                                                                 <  2:733400/147‑1 (MQ=11)
                                                    aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACGGCGCAAt                                                                                                 >  1:39675/1‑149 (MQ=37)
                                                    aaCGCGGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAAt                                                                                                 <  2:807207/149‑1 (MQ=38)
                                                         gTCTTATCGGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGTGTTTTTTTGTAGGTCCGACAACCCCCTAACGCCACACCCGACCTTACCCCCTGATCCTACGCGTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCt                                                                                            <  2:127551/149‑1 (MQ=2)
                                                         gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCt                                                                                            <  1:770845/149‑1 (MQ=38)
                                                         atCTTATCATCCCTACGTGAGACACGACATTGGCCTTATTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCt                                                                                            <  2:802653/148‑1 (MQ=11)
                                                               tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                                                                                                  >  1:323777/1‑105 (MQ=35)
                                                               tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                                                                                                  <  2:323777/105‑1 (MQ=35)
                                                                cAGGCTGACCCCAGACAGCGCAATCGCCTTATTTCGCCTGTTTTTGTAGGTCGGATAGGGCATATACCCCGCATACTACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAg                                                                                     <  2:847544/149‑1 (MQ=11)
                                                                 aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGGTTTCGc                                                                                    >  1:600095/1‑149 (MQ=255)
                                                                 aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGACAGCGCAATAGCCTGATTTAGc                                                                                    >  1:600118/1‑149 (MQ=38)
                                                                     cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGGCGCTATCTGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGa                                                                                >  1:658410/1‑149 (MQ=11)
                                                                      tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTCAt                                                                               >  1:221408/1‑149 (MQ=38)
                                                                         gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGTCCTACGCCAGACAGCGCAATAGCCTGATTTATCTTGAtttt                                                                            >  1:697886/1‑149 (MQ=255)
                                                                           cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGt                                                                          >  1:284033/1‑149 (MQ=37)
                                                                            tgcatgtgcAATAGGCTGATTTAGCTTGATTTGTTAGGTCCGATAAGGCGTTTACGCCGCATCCGATATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                         <  2:761309/142‑1 (MQ=11)
                                                                            agatggcggAATAGCCTGATTTAGCGTGATTTTTTAGGGCGTATAAGGCCTTTCCGCCGTATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                         <  2:750855/144‑1 (MQ=11)
                                                                            aGACAGCGGAATGCGCTGATTTGTCGTGATTTTTGAGTGCTGTGGAGGCGTTTACGCCGCACCAGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTa                                                                         <  2:384528/149‑1 (MQ=11)
                                                                            aGACAGCGATATGGGCTGATTTAGGGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                         <  2:512253/149‑1 (MQ=11)
                                                                              aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGAGCAATAGCCTTATCTCGCGTGATTTTATAga                                                                       >  1:279989/1‑148 (MQ=11)
                                                                                cgttcagtAGCCTTATGTAGCTTGATTTGGTAGTTCGTATCATGCGTTTACGCCTCTGCGTACATGATCGCCTGATGCACCGCTTGCCGCGTTGTATCAGGCTTACGCCAGACAGCGCAATATGCTGATTTAGCGTGATTTTGTAGGTc                                                                     <  2:671124/142‑1 (MQ=2)
                                                                                       tAGCCTGATTTAGCGTGATTTTGTAGGTCGGTTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAg                                                              <  1:448329/149‑1 (MQ=255)
                                                                                          ccTGAAATAACGTGATTTTGAAGCTCGTAGAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCTAGACAACGCAATAGCTTGCTTTCTCGTGATTTTGTATGTCGGATCAgccc                                                           >  2:419522/1‑146 (MQ=11)
                                                                                           cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGt                                                          <  1:864509/149‑1 (MQ=38)
                                                                                            tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                        >  1:166529/1‑86 (MQ=255)
                                                                                            tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                        <  2:166529/86‑1 (MQ=255)
                                                                                            tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAATGCGtt                                                         >  1:364510/1‑149 (MQ=255)
                                                                                              atctAGCGTTATGTTGTAGGTTGTATAATGCTTTTCCGCCGTATCCGACATGAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTa                                                       <  1:234984/146‑1 (MQ=11)
                                                                                              aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGa                                                                                                                                          <  2:399110/66‑1 (MQ=255)
                                                                                              aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                                                          >  1:399110/1‑66 (MQ=37)
                                                                                              aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTa                                                       <  2:637631/149‑1 (MQ=37)
                                                                                              aTTGAGCCTGAATTCGTAGGTCGGATAAGGTGTTTACGCCGTATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTa                                                       <  1:386495/149‑1 (MQ=11)
                                                                                               tttAGCGTGAGTTTGTAGGTCGGATAAGGCTTTTACGCCGTATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGTCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAt                                                      <  1:331321/149‑1 (MQ=11)
                                                                                                      tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                                   <  2:236059/65‑1 (MQ=40)
                                                                                                      tGATTTTGTAGGTCGGATAAGGCGTTTACATCCAGGCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGATCGGATAAGGCGTTTATGCCGCAt                                               >  1:437952/1‑149 (MQ=11)
                                                                                                              tAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGTCGCATTCGACATc                                       >  2:212359/1‑149 (MQ=255)
                                                                                                                tggcGGATAGGGCGTTTACGCCGCATCCGACGTCATCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCaa                                     <  1:225468/146‑1 (MQ=255)
                                                                                                                              tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATgcggcg                       >  1:273896/1‑149 (MQ=255)
                                                                                                                                   gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAtttt                                                                            >  1:365334/1‑91 (MQ=39)
                                                                                                                                    ccgcGTTCGGCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGTGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGATATCAACGCCTGATGCGACGCTTAAc                 <  2:21016/149‑1 (MQ=11)
                                                                                                                                      tcATGCGACATCAACGCCTGATCTAACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTATCGTGATTTCGTACGTCGGACAAGGCGTTTATGCCGCATCCTACATCAACGCCTGATGCGACGCTTAAcgc               <  1:760404/148‑1 (MQ=11)
                                                                                                                                      gggtCCGGCATCAACGCCTGATGCGACGCTTGCCGTGTGTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc               <  2:269641/146‑1 (MQ=255)
                                                                                                                                               aTCAACGCCTGATGCGACGCTTGCCGCGTTTCATCAGGCCCACGCCAGACAGCACAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCa      <  2:62615/149‑1 (MQ=255)
                                                                                                                                                tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGCCATCAACGCCTGAGACTACGCTTAACGCGCCTTATCAg     >  1:290954/1‑149 (MQ=11)
                                                                                                                                                 cAACGCCTGATGCACCGCTTGCCGCGTCTTATCAGTCCTACGCCAGACAGCGCAATAGTCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATAATAATTTAACGCCTGATGCGACGCTTAACGCGTCTTATCAgg    <  1:212359/149‑1 (MQ=11)
                                                                                                                                                  caCGCCTGATGCGACCCTTTCCGTGTCTTATCAGGCCTAAGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc   <  1:552511/148‑1 (MQ=255)
                                                                                                                                                  aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc   <  1:552587/149‑1 (MQ=32)
                                                                                                                                                   aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGCTGCGACGCTTACCGCGTCTTATCAGGcc  >  2:427038/1‑149 (MQ=9)
                                                                                                                                                    |                                                                                                                                                   
GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCC  >  NC_000913_3_bme_pgi/4295913‑4296208

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: