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% 259.6 / ‑3.5 58intergenic (‑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 . (29/28);  minor base G (0/1);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTG  >  NC_000913_3_bme_pgi/4235259‑4235551
                                                                                                                                                    |                                                                                                                                                
tCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                                                                  <  2:698913/149‑1 (MQ=255)
 cAATGCCTTCGATCAATGGGGCGTGGAACCGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGat                                                                                                                                                 <  1:321474/149‑2 (MQ=255)
  aaTGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGatc                                                                                                                                                <  1:559678/149‑3 (MQ=255)
     gCCGTCGATCAATGGGGCGGGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGatcatc                                                                                                                                            <  1:590093/149‑1 (MQ=255)
       cTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcgt                                                                                                                                          >  2:339356/1‑149 (MQ=255)
          cGATCAATGGGGCGTGGAACTGGGTTAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGa                                                                                                                                       <  1:822585/149‑1 (MQ=255)
           gATCAAGGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGat                                                                                                                                      <  2:334960/149‑1 (MQ=255)
               aaTGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGt                                                                                                                                  <  1:339356/149‑1 (MQ=255)
                 gggggCGTGTAACTGGGTAAGGAACCGGCGACGGAAATGTGGCCGGTTGTTTCGGCCAAAGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAg                                                                                                                                <  1:444850/148‑1 (MQ=255)
                      cGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCgg                                                                                                                           >  1:491500/1‑149 (MQ=255)
                        tGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGAt                                                                                                                         >  2:356950/1‑149 (MQ=255)
                            aCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAgg                                                                                                                     <  1:717358/149‑1 (MQ=255)
                            aCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAgg                                                                                                                     >  1:456794/1‑149 (MQ=255)
                                ggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                 <  2:456794/149‑1 (MQ=255)
                                     ggAACGGGCGACGGAACTTTGGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTGGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                            <  2:62555/149‑1 (MQ=255)
                                       aaCTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCACGGGCTTGTCGCGCATCTGCATGCGCACAGGACACCTTGATCATCGTCTGTATGTACGCCGGAGAAGGCGTTAAcgccg                                                                                                          >  1:798291/1‑149 (MQ=255)
                                        aCTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGACGAGGCGAGGGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAagcagc                                                                                                         >  1:120278/1‑149 (MQ=255)
                                            gCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCAtac                                                                                                     >  2:260789/1‑147 (MQ=255)
                                                  gAACTTTTGCCGGTTGTAGCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                               <  2:500583/149‑1 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCgct                                                                                           >  2:694500/1‑147 (MQ=255)
                                                        ttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTACACGCATGGGTTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                         <  1:818638/149‑1 (MQ=255)
                                                           ccGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGACGTTCACGCAGCATCCGGCAACCGCTGCCTg                                                                                      >  2:152989/1‑149 (MQ=255)
                                                            cGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                     >  1:111323/1‑149 (MQ=255)
                                                                    tggggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                             <  1:425257/147‑1 (MQ=255)
                                                                    tCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                             <  1:152989/149‑1 (MQ=255)
                                                                     cGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgcg                                                                            <  1:439504/149‑1 (MQ=255)
                                                                          aaggaaggGGCGAGTGGTCGGCACGGTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                       <  2:533655/149‑1 (MQ=255)
                                                                          aaggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                       <  2:823060/149‑1 (MQ=255)
                                                                            ggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGGCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                                                                     <  1:316333/149‑1 (MQ=255)
                                                                              aaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                                                                   >  1:465309/1‑149 (MQ=255)
                                                                               aggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                  <  1:356950/149‑1 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                              >  1:508469/1‑149 (MQ=255)
                                                                                           tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                           <  1:545183/144‑1 (MQ=255)
                                                                                           tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                           >  2:545183/1‑144 (MQ=255)
                                                                                           tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                      >  1:511410/1‑149 (MQ=255)
                                                                                             gcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAgg                                                    >  1:661915/1‑149 (MQ=255)
                                                                                               gACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                  <  2:508469/149‑1 (MQ=255)
                                                                                                 cGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                                >  2:571778/1‑149 (MQ=255)
                                                                                                  gCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGAt                                               >  2:758141/1‑149 (MQ=255)
                                                                                                    ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                             >  1:218735/1‑149 (MQ=255)
                                                                                                       cACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                          <  1:674741/149‑1 (MQ=255)
                                                                                                        aCGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                         >  1:184287/1‑149 (MQ=255)
                                                                                                         cGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGGCGATGCCGata                                        >  2:657913/1‑149 (MQ=255)
                                                                                                          ccAAGGGCTTGTCGCGCATCTGCATTCCCGAAGGAAAGCTTGATCCTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                       <  1:754578/148‑1 (MQ=25)
                                                                                                            aaGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                     >  2:632910/1‑149 (MQ=255)
                                                                                                             aGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTa                                    <  1:758141/149‑1 (MQ=255)
                                                                                                                  ttGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACACCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATTTGTAAATCg                               >  1:1191/1‑149 (MQ=37)
                                                                                                                    ggcgcgCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGCAGCCTGAATAAGGTGTGCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTa                             <  1:193455/147‑1 (MQ=37)
                                                                                                                     tCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                   <  2:216078/126‑1 (MQ=255)
                                                                                                                     tCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                   >  1:216078/1‑126 (MQ=255)
                                                                                                                     tCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAt                            <  2:724400/149‑1 (MQ=255)
                                                                                                                         gcATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                        <  1:567230/149‑1 (MQ=255)
                                                                                                                         gcATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAGGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                        <  1:113160/149‑1 (MQ=255)
                                                                                                                             cTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAt                    <  1:543224/149‑1 (MQ=255)
                                                                                                                              tGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                   >  1:387533/1‑149 (MQ=255)
                                                                                                                                 aTGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt                <  1:588921/149‑1 (MQ=255)
                                                                                                                                 aTGCGCGAAGGAAAGCTTGATCATCGTCGATATGGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGCTATGTACATCGTATTCGGCAATTAAt                >  2:600524/1‑149 (MQ=255)
                                                                                                                                  tGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGGTATGTACATCGTATTCGGCAATTAATa               >  1:439193/1‑149 (MQ=255)
                                                                                                                                     gcgAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt            >  1:74475/1‑149 (MQ=255)
                                                                                                                                         aGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCAACGCGAATGCCGTGTATCAGGCCTACAGGTCAACGCCGATATGTACATCGTATTCGGCAATTAATACATAGCa        <  1:785539/149‑1 (MQ=37)
                                                                                                                                               gCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattg  >  1:141499/1‑149 (MQ=255)
                                                                                                                                               gCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattg  >  1:68927/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTG  >  NC_000913_3_bme_pgi/4235259‑4235551

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: