Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 4,233,793 C→T 9.7% intergenic (‑560/‑2114) 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,233,7930CT 9.7% 145.6 / 4.0 64intergenic (‑560/‑2114)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (31/23);  new base T (3/3);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

TAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGAGA  >  NC_000913_3_bme_pgi/4233646‑4233938
                                                                                                                                                   |                                                                                                                                                 
tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTAACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAAAGACACGACGAAGCTTGAAGCAACGGACg                                                                                                                                                  >  2:465089/1‑149 (MQ=255)
    ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTGGCaag                                                                                                                                              >  1:490450/1‑148 (MQ=255)
      tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                                                                                                                            >  2:42952/1‑149 (MQ=255)
      tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                                                                                                                            >  1:808976/1‑149 (MQ=255)
           acTCAACATTTCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGTCGAGTCTTGAAGCAACGGTCGCAAAGCTGaaa                                                                                                                                       <  1:438352/149‑1 (MQ=255)
             tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                                                                                                                     >  1:412687/1‑149 (MQ=255)
             tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAACCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                                                                                                                     <  2:748361/149‑1 (MQ=255)
               aaCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGc                                                                                                                                   >  1:556646/1‑149 (MQ=255)
                aCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCa                                                                                                                                  <  1:718394/149‑1 (MQ=255)
                 cATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAt                                                                                                                                 <  1:42952/149‑1 (MQ=255)
                   ttACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACAAAGCTTAAAGCAACGGTCGCAAAGCTGAGAAAGCATTg                                                                                                                               >  2:321231/1‑149 (MQ=255)
                         tAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAATGGGCGGa                                                                                                                         >  1:650273/1‑149 (MQ=255)
                           actgCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGTAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTTGCAAAGCTGAAAAAGCATTGGGCGGAAt                                                                                                                       <  1:349210/146‑1 (MQ=255)
                           aCGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTTGCAAAGCTGAAAAAGCATTGGGCGGAAt                                                                                                                       >  1:554390/1‑149 (MQ=255)
                           aCGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAAt                                                                                                                       >  1:214802/1‑149 (MQ=255)
                              gCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCgg                                                                                                                    >  2:687977/1‑149 (MQ=255)
                                 ccaaaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGATGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGccc                                                                                                                 <  1:273915/147‑1 (MQ=255)
                                  tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                                                                                                                >  2:512626/1‑149 (MQ=255)
                                  tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                                                                                                                >  1:717666/1‑149 (MQ=255)
                                  tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                                                                                                                >  1:895520/1‑149 (MQ=255)
                                  tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGCCGGAATCGGcccc                                                                                                                >  2:570933/1‑149 (MQ=255)
                                  tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACACGACGAAGCTTAAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGAcccc                                                                                                                >  2:270259/1‑149 (MQ=255)
                                    aaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgc                                                                                                              >  1:235848/1‑149 (MQ=255)
                                    aaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgc                                                                                                              >  1:607053/1‑149 (MQ=255)
                                        catcatATTGTTATGTAACTGTTGTTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTACTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa                                                                                                          <  1:482755/149‑1 (MQ=255)
                                        cATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTTGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa                                                                                                          >  1:577472/1‑149 (MQ=255)
                                        cATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa                                                                                                          <  1:512626/149‑1 (MQ=255)
                                              aTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTTCTAATGGCAAGAGACACAACGAAGCATGAAACAACAATCGCAAAGCTGAAAAAGCATTGGGCGAAATCGGCCCCGCGCGGTATccg                                                                                                    >  1:10863/1‑149 (MQ=255)
                                                         aCTGGTGCTACAATCTTCCAAAGTTACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAAGCGACGAAGCTTGAAGAAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCTGACCCGCGCGGTATGGTCTCGGCATTCa                                                                                         >  2:103827/1‑149 (MQ=255)
                                                         aCTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCa                                                                                         <  1:734600/149‑1 (MQ=255)
                                                                 tACAATCTTCCAAAGTCACAATTATCAAAATCAGAAGAGTATTGCTCATGGCAAAAGACGCGACGAAGATAGAAGCAACGGTAGCACAGCTGAAAAACCATTGCGCGGAATCGGCCCCGCGCGCTCTGCGCGCGGCATTCAGCCACGAc                                                                                 >  2:205742/1‑149 (MQ=255)
                                                                  aCAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACACGACGAAGCTTAAAACAACGGTAGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCCCGATATGCGCGCGGCATTCATGGCCAAcc                                                                                >  1:18205/1‑149 (MQ=255)
                                                                        ttCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGGGATATGCGAGCGGCATTCAGCGCAGACCccggac                                                                          >  2:577019/1‑147 (MQ=255)
                                                                           cAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTAGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGt                                                                       >  1:274002/1‑149 (MQ=255)
                                                                                 caaaaTTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGATGCGACGAAGCTTGAAGCAACGGTTGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                                                                 <  1:527432/146‑1 (MQ=255)
                                                                                   caATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAAGCGGCCGCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCAGGCGa                                                               <  2:650273/149‑1 (MQ=255)
                                                                                    aaTTCTAAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAATAGCCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGat                                                              <  1:169808/149‑1 (MQ=255)
                                                                                     aTTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGata                                                             >  1:656584/1‑149 (MQ=255)
                                                                                         tcAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGTCATTCg                                                         >  1:610264/1‑149 (MQ=255)
                                                                                         tcAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCGGAAAAAGCATTGGGCGGAATCGGCACCGCGCGATATGCGCGCGGCATTCAGGCCAGACCCCGGCCGGTCAGGGGGAAAttgg                                                         >  2:201597/1‑147 (MQ=255)
                                                                                          taaaaTCAGTACAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATGGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGc                                                        <  1:321231/148‑1 (MQ=255)
                                                                                           aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGct                                                       <  2:17724/149‑1 (MQ=255)
                                                                                               tCAGAAGAGTATTTCTAATGGTAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGc                                                   <  1:143949/149‑1 (MQ=255)
                                                                                                  gaagaGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGCCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTg                                                <  2:835195/149‑1 (MQ=255)
                                                                                                   aagaGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTgg                                               >  1:289215/1‑149 (MQ=255)
                                                                                                    agagTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGga                                              <  1:126289/149‑1 (MQ=255)
                                                                                                      agTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGCCCCCGGCCGGTTGGGGCGCTATTCACCCTGCGTGgatg                                            >  2:268258/1‑149 (MQ=255)
                                                                                                        tATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTGGCAAAGCTGAAAAAGCATTGAGAGGAATCGGCCCCGCGCGATATGCACGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCGGCCTGgctgat                                          >  2:137381/1‑149 (MQ=255)
                                                                                                        tATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCCCTGCCTGgatgct                                          >  2:654782/1‑147 (MQ=255)
                                                                                                            tctgattgCAAGAGACGCCTCGAGGCTTGAAGCAACGGTCGAAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCTCGGCATTCAGCGCCGTCCCCGGGCGGTTCGGTCGATATTCGCTCTGCCTGGATGAtcttc                                      <  2:166018/142‑1 (MQ=255)
                                                                                                             cTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCGGGATGAtcgtcc                                     >  2:523061/1‑148 (MQ=255)
                                                                                                             cTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCGGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGatcatcc                                     >  2:523020/1‑148 (MQ=255)
                                                                                                                   gCAAGAGACGCGACGAAGCTTGAAGCAACGGTTGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGa                               <  1:354388/149‑1 (MQ=255)
                                                                                                                          aacagactaAGCTTGAAGAAACGGGCTCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCg                        >  2:647313/5‑149 (MQ=255)
                                                                                                                           cgcgACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCAACCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGGCGa                       >  1:558880/1‑149 (MQ=255)
                                                                                                                            gcgACGAAGCTTGAAGCAACGGTTGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGaa                      >  1:231400/1‑149 (MQ=255)
                                                                                                                            gcgACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGaa                      <  1:570933/149‑1 (MQ=255)
                                                                                                                            gcgACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAAGCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGCCCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGaa                      <  2:745635/149‑1 (MQ=255)
                                                                                                                            gacacgaAGTTTAAAGCAACGGTCGCAAAGCTGAAAAAGCTTTGGGCGGAATCGGCCCCGCGCGATATGCGGGCGGAATTCAGCGCCGACCCCGGCCGGTTCGGGCGCTCTTTCTGCTGCCTGGATGATCTTCTGTTCGACTGGTCGaa                      <  1:436630/146‑1 (MQ=255)
                                                                                                                                        gAAGCAATGGTCGCAAAGTTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCTACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTaa          <  2:802615/149‑1 (MQ=255)
                                                                                                                                        gAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCTGAATCGGCCCCGCGCGATATGCGCGCGTCATTCAGCGCCGACCCCGGCAGCTTGCGTCTATATTCGCTGTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTaa          <  1:577019/149‑1 (MQ=255)
                                                                                                                                           gCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAacga       >  2:886212/1‑149 (MQ=255)
                                                                                                                                              aCGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAAcgacga    <  1:596508/149‑1 (MQ=255)
                                                                                                                                                ggTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTAGTCTTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACgaga  <  1:397367/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGAGA  >  NC_000913_3_bme_pgi/4233646‑4233938

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: