Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,235,408 Δ1 bp intergenic (‑2175/‑499) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,235,4070A.98.3% 252.2 / ‑2.9 59intergenic (‑2174/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (34/23);  minor base T (0/1);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_bme_pgi/4235259‑4235550
                                                                                                                                                    |                                                                                                                                               
tCAAGGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                                                                 <  2:718489/149‑1 (MQ=255)
   aTGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGatca                                                                                                                                              >  2:632684/1‑146 (MQ=255)
     gcgctCGATCAAGTGGTCGTGGAACTGGGTAAGTAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGACGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTTTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑atcatc                                                                                                                                           <  2:361158/145‑1 (MQ=255)
        ttCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCCTCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgta                                                                                                                                        >  1:255890/1‑148 (MQ=255)
        ttCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGACGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTACCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGC‑TGGATCAtcgcc                                                                                                                                        >  1:273026/1‑147 (MQ=255)
            aTCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCTCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCCTCGata                                                                                                                                    >  2:527806/1‑149 (MQ=255)
            aTCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGAAAGGAAGGGGCGAGTGGTCGCGAGGCTTGCACGCGAGGGCTTGTCGCGCATCTGCATGCGCGAAGGCAAGCTTG‑ATCCTCGTCGata                                                                                                                                    >  2:304072/1‑149 (MQ=255)
               aaTGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGt                                                                                                                                 <  2:489895/149‑1 (MQ=255)
                   gggCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGc                                                                                                                             >  2:651696/1‑149 (MQ=255)
                    ggggTGGAACGGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGcc                                                                                                                            <  2:698088/146‑1 (MQ=255)
                              tGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAATGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                                                                  >  2:163535/1‑149 (MQ=255)
                                   aaGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                                                             >  1:506121/1‑149 (MQ=255)
                                       aaCTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                         <  1:600988/149‑1 (MQ=255)
                                            gCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGATGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                    <  1:454407/149‑1 (MQ=255)
                                                 ggAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGAATACGGCac                                                                                               >  2:777090/1‑148 (MQ=255)
                                                     cTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGa                                                                                           <  2:102677/149‑1 (MQ=255)
                                                            cGGTTGGATCGGGCAAGGAAGGGGGGAGTGGTCGCGACGCTGCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                    <  1:229416/149‑1 (MQ=255)
                                                                    tCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAGGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGCTGcgccgc                                                                            >  1:242551/1‑149 (MQ=255)
                                                                      gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGCCGCgg                                                                          >  1:594193/1‑149 (MQ=255)
                                                                      gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCgg                                                                          >  1:760591/1‑149 (MQ=255)
                                                                          aaggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                      >  1:312538/1‑149 (MQ=255)
                                                                            ggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                                                                    <  1:651696/149‑1 (MQ=255)
                                                                               attgccGAGGGGTCGCGACGGTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                 <  2:628431/144‑1 (MQ=37)
                                                                                ggggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGCCGCGGTCTCGTCTTa                                                                >  1:175192/1‑149 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                          >  2:391542/1‑120 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                          <  1:391542/120‑1 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                             >  1:599474/1‑149 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAAGCGGTCGCGTCTTATCa                                                             >  1:26190/1‑149 (MQ=255)
                                                                                     gAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg                                                           >  1:480902/1‑149 (MQ=255)
                                                                                        tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGGGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAGCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt                                                        <  2:484943/149‑1 (MQ=255)
                                                                                        tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                              <  2:830816/111‑1 (MQ=255)
                                                                                        tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                              >  1:830816/1‑111 (MQ=255)
                                                                                          gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                      >  1:633901/1‑149 (MQ=255)
                                                                                           tCGCGACGCTTCCCCGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                     <  2:212246/149‑1 (MQ=255)
                                                                                           tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                     >  2:178437/1‑149 (MQ=255)
                                                                                              cgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                  <  1:51860/149‑1 (MQ=255)
                                                                                               gACGCTTCCACGCTGTGGGTGGTCGCGCAGCTGCAGGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                 <  2:757165/149‑1 (MQ=25)
                                                                                               gACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                     >  1:508818/1‑81 (MQ=255)
                                                                                               gACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                     <  2:508818/81‑1 (MQ=255)
                                                                                                 cGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                                                                  <  2:572124/130‑1 (MQ=255)
                                                                                                 cGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                                                                  >  1:572124/1‑130 (MQ=255)
                                                                                                  gCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                        <  2:36709/123‑1 (MQ=255)
                                                                                                  gCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                        >  1:36709/1‑123 (MQ=255)
                                                                                                   cTTCCACGCAAGGGCTTGTTGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCAGCAACCGAGGCCTGATGCGACGCGGGCGCGTCTTATAAGGCCTTCAGGTCGGTg                                             >  2:706451/1‑149 (MQ=37)
                                                                                                    ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                                                                 <  2:366092/49‑1 (MQ=255)
                                                                                                    ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                                                                 >  1:366092/1‑49 (MQ=255)
                                                                                                           ccaGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTT‑TATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                     <  1:232122/147‑1 (MQ=255)
                                                                                                                  ttGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                               >  1:59554/1‑148 (MQ=255)
                                                                                                                  ttGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                               <  2:59554/148‑1 (MQ=255)
                                                                                                                    gTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                   <  2:292727/94‑1 (MQ=255)
                                                                                                                    gTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                   >  1:292727/1‑94 (MQ=255)
                                                                                                                      cgcgCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAtt                          >  2:655687/1‑149 (MQ=255)
                                                                                                                         gcATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGTCGAAATGTACCTCGTAGTCgg                       >  2:385087/1‑149 (MQ=25)
                                                                                                                           aTCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                     >  1:463015/1‑149 (MQ=255)
                                                                                                                              tGCATGCGCGAAGGAAAGCTTG‑ATCAGCGACGGTAACTAAGCAGAATAAGGGGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGCCGCGGTCTCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTTTTCGGCAAtt                  >  2:392876/1‑149 (MQ=25)
                                                                                                                                 aTGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGACGCGTTTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt               <  1:163535/149‑1 (MQ=255)
                                                                                                                                   gcgcGAAGGAAAGCTTG‑ATCACCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGACTGATGCGACGCGGCCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAc             <  2:514439/149‑1 (MQ=255)
                                                                                                                                     gcgAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt           >  2:789423/1‑149 (MQ=255)
                                                                                                                                      cgAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATACCAAGCGGGCGCCGATTATGAACACCACGGATCGATGCCGGTAAATCCATCGCGTTCGCCAATCAATACaaa          >  2:789693/1‑147 (MQ=16)
                                                                                                                                          ggAAAGCTTG‑ATGATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTTGGCAATTAATACATAGCAc      <  2:368496/149‑1 (MQ=255)
                                                                                                                                          ggAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTACCAGGCCTACAGGTCGATGCCGATATCAACATCCTATTCGCCACTTAATACAAAacac      >  2:329027/1‑149 (MQ=25)
                                                                                                                                           gAAAGCTTG‑ATCATCGTCGATATGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCAATGCCTAAGAGGGCGCTGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGGATTAGGCAATTAATACATAGCACg     >  2:719756/1‑149 (MQ=25)
                                                                                                                                            aaaGCTTG‑ATCATCGTCGATATGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACga    >  2:638713/1‑149 (MQ=255)
                                                                                                                                              aGCTTG‑ATCATCGTTGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgatt  <  1:632684/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_bme_pgi/4235259‑4235550

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