Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 4,233,708 C→T 100% intergenic (‑475/‑2199) 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,7080CT98.6% 187.5 / ‑5.4 64intergenic (‑475/‑2199)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (41/22);  minor base G (0/1);  total (41/23)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCA  >  NC_000913_3_bme_pgi/4233566‑4233851
                                                                                                                                              |                                                                                                                                               
ggTTGTCAACGATGGGGTCATGCGGATATTTCATCCACTCCTGGCGGACAGTAGTTCAGCTAATCAATGCTTCACTGCGCTAAGGGTGTACACCCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTTACTGGTGCTACa                                                                                                                                           >  2:288406/1‑149 (MQ=255)
     tCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACACAACATTACGCTAACGGCACTAAAACAATCACCTTTTTATGTGACTGGTGCTACAATCtt                                                                                                                                      >  1:17470/1‑149 (MQ=255)
          gagggggTCATGCGGATTTTTCATCCACTCCTTGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCaaa                                                                                                                                 <  1:105588/146‑1 (MQ=255)
          gATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCTCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCaaa                                                                                                                                 >  1:886278/1‑149 (MQ=255)
           aTGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTCGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAg                                                                                                                                <  1:794694/149‑1 (MQ=255)
                gTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGAGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGTTACAAACTTCCAAAGTcaca                                                                                                                           >  1:564478/1‑149 (MQ=255)
                 tCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTCCACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGAATGGTGCTACAATATTCCAAAGTCACaa                                                                                                                          >  1:775780/1‑149 (MQ=255)
                      cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTTCGCTAAGGGTTTACACTTAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                                                     <  1:87584/149‑1 (MQ=255)
                      cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGGTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                                                     >  2:481656/1‑149 (MQ=255)
                               cATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCCCAATTCTCCAAATCaca                                                                                                            >  2:72993/1‑147 (MQ=255)
                                aTCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaa                                                                                                           <  1:58096/149‑1 (MQ=255)
                                 tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaag                                                                                                          >  2:718394/1‑149 (MQ=255)
                                  ccACTCCTGGCGGGTAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                                         <  2:605034/149‑1 (MQ=255)
                                     ctacTGGCGGTCAGTATTTTATCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                                                                                      <  1:358184/146‑1 (MQ=255)
                                     cTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                                                                                      <  2:445737/149‑1 (MQ=255)
                                      tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTGAGGGTTTACACTAAAGATAACGATAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTACAAAGTCACAATTCTCAAAATCAGAAGAGTAt                                                                                                     >  1:169599/1‑149 (MQ=255)
                                       ccTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTAACAATTCTCAAAATCAGAAGAGTAtt                                                                                                    >  1:17724/1‑149 (MQ=255)
                                       ccTGGCCGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAATGTTTTACACTTAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAtt                                                                                                    <  1:147911/149‑1 (MQ=255)
                                             ggTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTCTTAATTCCCAATCTTGCTGTCGCTTGGACTCACATCCTGCAAAGTCACAATTCTCACAAACAGAAGAGTATTGCTAAt                                                                                              <  2:865603/149‑1 (MQ=255)
                                               tCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAAACATCAAAAACTTCTGTAGCTGTTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAAAAGAGTATTGCTAATgg                                                                                            >  2:546135/1‑149 (MQ=255)
                                                      ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAAACTTCCAAAGTCACAATACACAAAATCAGAAGAGTATTGCTAATGGCAAgaga                                                                                     >  1:582495/1‑149 (MQ=255)
                                                        cAGCTAATCATTTGTTCTCTGGGCTAAGGGTTTACACTCAACATTACGCTAACGGGATTAAAACCATCACATTTTTCTGTGATTGGTGCTACAAGCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcg                                                                                   <  2:96647/149‑1 (MQ=255)
                                                          gCTAATAAATGCTTCACGGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATAAAAAGAGTCTTGCTAATGGCAAGACAcgcg                                                                                 >  2:527375/1‑149 (MQ=255)
                                                               taaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGaa                                                                            >  1:870024/1‑149 (MQ=255)
                                                                aaaTGCTTCGCTGCGCTGAGGGTTTACACTAAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGACCCCCTACGCGAAGCTCACAATTCTTAAAATCAGAAGAGTATTGCTAATTGCAAGAGACGCGACGAAg                                                                           <  2:412703/149‑1 (MQ=255)
                                                                 aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCGCTAAAACCATAACATTTTTCTGTGACTGGTGCTACGATCAGCTCACGGTGGTCGCCGTATAATTAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGc                                                                          >  2:864454/1‑149 (MQ=255)
                                                                  aTGCCTCACTGCGCTAAGAGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAAGTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCt                                                                         <  2:648806/149‑1 (MQ=255)
                                                                      ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGTATTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCATAATTCTCAAAAGCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGaa                                                                     <  1:324940/149‑1 (MQ=255)
                                                                      ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACACTTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTGGaa                                                                     >  1:521807/1‑149 (MQ=255)
                                                                       tCACTGCGCTAAGTGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAg                                                                    <  1:737386/149‑1 (MQ=255)
                                                                         aCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCCAAATCAGAAGAGTATTGCTAATGTCAAGAGACGCGACGAAGCTTGAATCa                                                                  >  2:527432/1‑149 (MQ=255)
                                                                          cTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGAAGAACCTTGAAGCac                                                                 >  1:506557/1‑148 (MQ=255)
                                                                           tGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAAc                                                                >  1:661112/1‑149 (MQ=255)
                                                                              gcTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAACCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGt                                                             >  2:354388/1‑149 (MQ=255)
                                                                                tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTAACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAAAGACACGACGAAGCTTGAAGCAACGGACg                                                           >  2:465089/1‑149 (MQ=255)
                                                                                    ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTGGCaag                                                       >  1:490450/1‑148 (MQ=255)
                                                                                      tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                                     >  1:808976/1‑149 (MQ=255)
                                                                                      tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                                     >  2:42952/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                         >  1:895520/1‑149 (MQ=255)
                                                                                                                  tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                         >  2:512626/1‑149 (MQ=255)
                                                                                                                  tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                         >  1:717666/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)
                                                                                                                                              |                                                                                                                                               
GGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCA  >  NC_000913_3_bme_pgi/4233566‑4233851

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