Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,794 G→A intergenic (‑561/‑2113) 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,7940GA96.4% 368.6 / 2.0 110intergenic (‑561/‑2113)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (57/49);  total (59/51)
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.06e-01

Only 100 of 110 total aligned reads displayed.

AGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGA  >  NC_000913_3_bme_pgi/4233648‑4233936
                                                                                                                                                  |                                                                                                                                              
aGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACaa                                                                                                                                             >  4:26070/1‑150 (MQ=255)
aGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACaa                                                                                                                                             >  3:659249/1‑150 (MQ=255)
   gTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGc                                                                                                                                          >  2:1145797/1‑150 (MQ=255)
   gTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGc                                                                                                                                          <  2:298326/150‑1 (MQ=255)
    tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCt                                                                                                                                         >  3:308809/1‑150 (MQ=255)
       acacTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGaa                                                                                                                                      <  2:323703/150‑1 (MQ=255)
          cTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGaaaaa                                                                                                                                   <  4:308809/150‑1 (MQ=255)
           tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGaa                                                                                                                                      >  2:392418/1‑146 (MQ=255)
           tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGaa                                                                                                                                      <  1:392418/146‑1 (MQ=255)
           tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAg                                                                                                                                  >  2:638539/1‑150 (MQ=255)
           tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAg                                                                                                                                  >  3:298889/1‑150 (MQ=255)
             aaCATTACGCTAACGGCACTAAAACCATCACATTTTTCTTTGACTGGTTCTACAATCTTCCAAATTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCa                                                                                                                                >  2:787187/1‑150 (MQ=255)
             aaCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCa                                                                                                                                >  3:560059/1‑150 (MQ=255)
               cATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCAtt                                                                                                                              >  2:539625/1‑150 (MQ=255)
                       tAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGaa                                                                                                                      <  3:26070/150‑1 (MQ=255)
                         aCGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATc                                                                                                                    <  2:294410/150‑1 (MQ=255)
                         aCGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATc                                                                                                                    <  4:172354/150‑1 (MQ=255)
                         aCGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc                                                                                                                    <  4:560059/150‑1 (MQ=255)
                            gCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGc                                                                                                                 >  2:206020/1‑150 (MQ=255)
                            gCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGc                                                                                                                 >  1:837944/1‑150 (MQ=255)
                                tAAAACCATCACATTTTTCTGTAACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTCAAAAAGCATTGGGCGGGATCGGccccc                                                                                                             >  2:746082/1‑149 (MQ=255)
                                    cccATCACATTTTTCTGTGACTGGTGCTACAATCTGCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgcg                                                                                                         <  3:354359/149‑1 (MQ=255)
                                          acaTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgc                                                                                                   <  3:438727/150‑1 (MQ=255)
                                            aTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgc                                                                                                 <  1:638539/150‑1 (MQ=255)
                                            aTTTTTCTGTGACTGGTGCGACAAGCTTCCAAAGTCACAATTCTCAAAACCAGCAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgc                                                                                                 <  1:746082/150‑1 (MQ=255)
                                               tttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGcgcgcg                                                                                              <  1:380010/150‑1 (MQ=255)
                                                   tgtgACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAt                                                                                          >  2:140289/1‑150 (MQ=255)
                                                         tGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcg                                                                                    <  1:130164/150‑1 (MQ=255)
                                                           gtgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                                                                  >  2:462648/3‑150 (MQ=255)
                                                            tgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGACGCTTGAAGCCACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCg                                                                                 <  1:1063313/149‑1 (MQ=255)
                                                               tACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                                                                                                                  <  3:88954/114‑1 (MQ=255)
                                                               tACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAcc                                                                              >  2:1083259/1‑150 (MQ=255)
                                                                aCAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTg                                                                                                                             >  1:421276/1‑102 (MQ=255)
                                                                aCAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTg                                                                                                                             <  2:421276/102‑1 (MQ=255)
                                                                  aaTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcg                                                                                                           <  4:87420/118‑1 (MQ=255)
                                                                  aaTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcg                                                                                                           >  3:87420/1‑118 (MQ=255)
                                                                  aaTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccg                                                                           <  1:1145797/150‑1 (MQ=255)
                                                                   aTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccgg                                                                          <  1:206020/150‑1 (MQ=255)
                                                                         cAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCAtt                                                                                                                              <  3:471630/92‑1 (MQ=255)
                                                                         cAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCAtt                                                                                                                              >  4:471630/1‑92 (MQ=255)
                                                                          aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTc                                                                   <  1:462648/150‑1 (MQ=255)
                                                                           aaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcgc                                                                                   <  4:523905/133‑1 (MQ=255)
                                                                           aaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcgc                                                                                   >  3:523905/1‑133 (MQ=255)
                                                                              tCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                                                               <  4:555456/150‑1 (MQ=255)
                                                                                  aaTTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGata                                                           >  2:194888/1‑150 (MQ=255)
                                                                                       tcAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGc                                                      >  1:646407/1‑150 (MQ=255)
                                                                                        cAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGct                                                     >  1:516549/1‑150 (MQ=255)
                                                                                         aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGccc                                                                                                               <  2:1120660/91‑1 (MQ=255)
                                                                                         aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGccc                                                                                                               >  1:1120660/1‑91 (MQ=255)
                                                                                         aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCgg                                                                                             >  4:226262/1‑109 (MQ=255)
                                                                                         aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCgg                                                                                             <  3:226262/109‑1 (MQ=255)
                                                                                          aaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctct                                                   >  1:945359/1‑150 (MQ=255)
                                                                                          aaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctct                                                   <  2:1066190/150‑1 (MQ=255)
                                                                                               agaagaGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTg                                              >  2:1163454/1‑150 (MQ=255)
                                                                                               agaagaGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTg                                              >  3:381660/1‑150 (MQ=255)
                                                                                                 aagaGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGga                                            >  3:279931/1‑150 (MQ=255)
                                                                                                   gagTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCt                                                                                                                                         <  1:161290/55‑1 (MQ=255)
                                                                                                   gagTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCt                                                                                                                                         >  2:161290/1‑55 (MQ=255)
                                                                                                      tATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtc                                       <  1:194888/150‑1 (MQ=255)
                                                                                                      tATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtc                                       <  2:1249714/150‑1 (MQ=255)
                                                                                                        ttGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggccg                                                                       <  4:68845/116‑1 (MQ=255)
                                                                                                        ttGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggccg                                                                       >  3:68845/1‑116 (MQ=255)
                                                                                                          gCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtcttct                                   <  2:48494/150‑1 (MQ=255)
                                                                                                           cTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGc                                                 <  2:1028189/135‑1 (MQ=255)
                                                                                                           cTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGc                                                 >  1:1028189/1‑135 (MQ=255)
                                                                                                            tAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGaaaaa                                                                                                                                   <  3:299474/52‑1 (MQ=255)
                                                                                                            tAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGaaaaa                                                                                                                                   >  4:299474/1‑52 (MQ=255)
                                                                                                            tAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGt                                 >  1:1074337/1‑150 (MQ=255)
                                                                                                            tAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTCCGGGCGACCTCCGCTCCGCCCGGATGATCTTCTGt                                 >  2:249531/1‑150 (MQ=255)
                                                                                                             aaTGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtctt                                     >  2:467542/1‑145 (MQ=255)
                                                                                                             aaTGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtctt                                     <  1:467542/145‑1 (MQ=255)
                                                                                                               tGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATc                                                                                                                    <  1:428661/64‑1 (MQ=255)
                                                                                                               tGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATc                                                                                                                    >  2:428661/1‑64 (MQ=255)
                                                                                                                   aaGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTg                          >  1:306862/1‑150 (MQ=255)
                                                                                                                        aCGCGACGAAGCTTGAAGCAACGGTCACAATGCTGAAAAAGTATTGGGCGGAATCGGCCCCGCGCGTTATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATTATCTTCTGTTCGACTGGTCgg                     <  2:1010557/150‑2 (MQ=255)
                                                                                                                        aCGCGACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggcc                                                                        >  2:920036/1‑99 (MQ=255)
                                                                                                                          gcgACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCgg                                                                 >  2:688491/1‑104 (MQ=255)
                                                                                                                          gcgACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCgg                                                                 <  1:688491/104‑1 (MQ=255)
                                                                                                                          gcgACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAg                   >  1:851116/1‑150 (MQ=255)
                                                                                                                           cgacgaAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGt                  >  2:377615/1‑150 (MQ=255)
                                                                                                                            gacgaAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTg                                                  >  1:443248/1‑117 (MQ=255)
                                                                                                                            gacgaAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTg                                                  <  2:443248/117‑1 (MQ=255)
                                                                                                                            gacgaAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTg                 >  2:674872/1‑150 (MQ=255)
                                                                                                                               gaAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGc                                                      >  1:980418/1‑110 (MQ=255)
                                                                                                                               gaAGCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGc                                                      <  2:980418/110‑1 (MQ=255)
                                                                                                                                  gCTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCgg             >  2:227240/1‑148 (MQ=255)
                                                                                                                                   cTTGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGtt          >  4:152055/1‑150 (MQ=255)
                                                                                                                                    ttGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTa         >  2:1192164/1‑150 (MQ=255)
                                                                                                                                    ttGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTa         >  1:713737/1‑150 (MQ=255)
                                                                                                                                    ttGAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGACGGGTTa         >  1:705680/1‑150 (MQ=255)
                                                                                                                                      gAAGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAac       >  1:776336/1‑150 (MQ=255)
                                                                                                                                        aGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGt                  >  2:745959/1‑137 (MQ=255)
                                                                                                                                        aGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGt                  <  1:745959/137‑1 (MQ=255)
                                                                                                                                        aGCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAacct     >  3:556594/1‑148 (MQ=255)
                                                                                                                                         gCAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAacgac    >  3:727371/1‑150 (MQ=255)
                                                                                                                                         gCAACGGTCACAAAGCTGAAAAAGCATGGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGTTCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAtcgac    <  4:298889/150‑1 (MQ=255)
                                                                                                                                          cAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtctt                                     <  2:382035/116‑1 (MQ=255)
                                                                                                                                          cAACGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtctt                                     >  1:382035/1‑116 (MQ=255)
                                                                                                                                           aaCGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAAcgacga  <  2:1074337/150‑1 (MQ=255)
                                                                                                                                           aaCGGTCACAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAAcgacga  <  1:164349/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
AGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGA  >  NC_000913_3_bme_pgi/4233648‑4233936

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