Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,689 C→T intergenic (‑456/‑2218) 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,6890CT99.2% 474.1 / ‑4.6 131intergenic (‑456/‑2218)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (73/56);  total (74/57)
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.34e-01

Only 100 of 133 total aligned reads displayed.

TATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCG  >  NC_000913_3_bme_pgi/4233543‑4233833
                                                                                                                                                  |                                                                                                                                                
tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATtaca                                                                                                                                               <  3:943653/150‑1 (MQ=255)
tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATtaca                                                                                                                                               <  1:121970/150‑1 (MQ=255)
 ataaataaAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACAt                                                                                                                                              >  3:829959/1‑150 (MQ=255)
  taaataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACAtt                                                                                                                                             >  2:264909/1‑150 (MQ=255)
    aataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACAtttt                                                                                                                                           >  2:353586/1‑150 (MQ=255)
     ataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACAttttt                                                                                                                                          <  3:413148/150‑1 (MQ=255)
     ataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCAATGAGCTAAGGGTTTACACTCAACATTAAGCTAACGGCACTAAAACCATTACAttttt                                                                                                                                          >  1:1484318/1‑150 (MQ=255)
      tacaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGTTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTc                                                                                                                                         <  1:1274824/147‑1 (MQ=255)
        aaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCtg                                                                                                                                       <  1:1491770/150‑1 (MQ=255)
         aGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCtgt                                                                                                                                      <  3:961706/150‑1 (MQ=255)
           ggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGa                                                                                                                                    <  1:1116153/150‑1 (MQ=255)
           ggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGa                                                                                                                                    >  1:1213078/1‑150 (MQ=255)
              gAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTg                                                                                                                                 >  2:1155084/1‑150 (MQ=255)
               aGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTCCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTgg                                                                                                                                <  1:1013725/150‑1 (MQ=255)
                 cGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTggtg                                                                                                                              >  1:932895/1‑148 (MQ=255)
                        cttGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAAt                                                                                                                       >  4:839681/2‑150 (MQ=255)
                            tCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTc                                                                                                                   <  1:1045268/150‑1 (MQ=255)
                              aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCa                                                                                                                 >  1:487982/1‑150 (MQ=255)
                              aaCGATGGGGTCAGGCGGAGTTTTCATCCACTTGTGGCGGTCAGTAGTTGAGCGAATAAAAGCTTCACTGCGCTAAGGGGTTACACTCAACGTCACGCTAACGGCACTAAAACCATTACATCTTTCTGTGACTGTTGCTACAATCTTCCa                                                                                                                 >  2:1407138/1‑150 (MQ=255)
                                 gATGGGGGCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAg                                                                                                              >  4:403413/1‑150 (MQ=255)
                                  aTGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGt                                                                                                             <  2:487982/150‑1 (MQ=255)
                                      ggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTcaca                                                                                                         >  1:626390/1‑150 (MQ=255)
                                         cATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAAtt                                                                                                      <  1:69885/150‑1 (MQ=255)
                                           tGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAACCCCTTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtct                                                                                                    <  1:1133322/150‑1 (MQ=255)
                                            gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                                   >  2:453309/1‑150 (MQ=255)
                                            gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                                   <  2:1459435/150‑1 (MQ=255)
                                             cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATc                                                                                                                      >  2:1352926/1‑130 (MQ=255)
                                             cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATc                                                                                                                      <  1:1352926/130‑1 (MQ=255)
                                              ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaa                                                                                                 >  1:1593200/1‑150 (MQ=255)
                                              ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaa                                                                                                 >  3:851803/1‑150 (MQ=255)
                                               gATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaaa                                                                                                <  2:1593200/150‑1 (MQ=255)
                                                   tttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                            >  4:549723/1‑150 (MQ=255)
                                                    ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCag                                                                                           >  2:774701/1‑150 (MQ=255)
                                                    ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCag                                                                                           >  1:455564/1‑150 (MQ=255)
                                                    tgcATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCag                                                                                           <  2:1031308/148‑1 (MQ=255)
                                                     tCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCaga                                                                                          <  1:264909/150‑1 (MQ=255)
                                                        tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                       <  4:1018140/150‑1 (MQ=255)
                                                        tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                       >  3:647277/1‑150 (MQ=255)
                                                            cTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTTCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAt                                                                                   >  4:46756/1‑150 (MQ=255)
                                                               cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                >  1:1541561/1‑150 (MQ=255)
                                                               cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                <  2:304563/150‑1 (MQ=255)
                                                               cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                >  1:1084457/1‑150 (MQ=255)
                                                                tGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCt                                                                               <  3:344195/150‑1 (MQ=255)
                                                                  gCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTggt                                                                                                                               <  3:270525/100‑2 (MQ=255)
                                                                  gCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTggt                                                                                                                               >  4:270525/1‑99 (MQ=255)
                                                                   cGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAAt                                                                            <  2:1512025/150‑1 (MQ=255)
                                                                   cGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAAt                                                                            <  3:49649/150‑1 (MQ=255)
                                                                    ggTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATg                                                                           <  4:156955/150‑1 (MQ=255)
                                                                      tCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGc                                                                         >  1:485974/1‑150 (MQ=255)
                                                                        agtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTa                                                                                                                           <  1:250538/98‑1 (MQ=255)
                                                                           agtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAagag                                                                    <  2:485974/150‑1 (MQ=255)
                                                                            gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAgaga                                                                   >  3:550030/1‑150 (MQ=255)
                                                                            gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAgaga                                                                   <  1:1247300/150‑1 (MQ=255)
                                                                                aGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAAAATTACGCTAACGGAACTAAAACCATTACATTTTTCTGTGAATGGTGCTACAATCTTCCAAAGTCACAATTCTCAAACTCAGACGACTATTGCAAATGACCACAGAcgct                                                               >  4:874605/1‑149 (MQ=255)
                                                                                  cTAATAAATGCTTCACTGCGCTAAGGGACTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGGATTGCTAATGGCAAGAGACGctac                                                             >  1:1104860/1‑150 (MQ=255)
                                                                                   taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGACGAGTATTGCTAATGGCAAGAGACGctacg                                                            >  1:774133/1‑150 (MQ=255)
                                                                                   taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacg                                                            >  3:754375/1‑150 (MQ=255)
                                                                                   taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacg                                                            >  3:332576/1‑150 (MQ=255)
                                                                                     agaaaTGCTTCACTGCGCGAAGGGTTTACTCTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACCATCTTCCAAAGTCACAATTCTCAAAAGCAGTAGAGTATTGCTAATGGCAAGAGACGCTACGaa                                                          <  1:1436926/148‑1 (MQ=255)
                                                                                         aTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCtt                                                      <  1:1518886/150‑1 (MQ=255)
                                                                                            cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                   <  1:518526/150‑1 (MQ=255)
                                                                                            cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                   <  2:527762/150‑1 (MQ=255)
                                                                                             ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                  >  3:600930/1‑150 (MQ=255)
                                                                                             ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                  <  2:1112064/150‑1 (MQ=255)
                                                                                                  tGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACg                                             >  3:229868/1‑150 (MQ=255)
                                                                                                       tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTa                                                                                                                           >  2:838345/1‑67 (MQ=255)
                                                                                                       tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTa                                                                                                                           <  1:838345/67‑1 (MQ=255)
                                                                                                       tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                        >  2:466516/1‑150 (MQ=255)
                                                                                                       tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                        >  1:1175400/1‑150 (MQ=255)
                                                                                                         aGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCaa                                      >  3:665860/1‑150 (MQ=255)
                                                                                                          gggTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGc                                                                         >  1:1327401/1‑114 (MQ=255)
                                                                                                          gggTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGc                                                                         <  2:1327401/114‑1 (MQ=255)
                                                                                                           ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                    <  1:1155084/150‑1 (MQ=255)
                                                                                                            gTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGc                                   <  4:829959/150‑1 (MQ=255)
                                                                                                            gTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGc                                   >  2:1027921/1‑150 (MQ=255)
                                                                                                            gTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGc                                   <  2:455564/150‑1 (MQ=255)
                                                                                                               tACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAagag                                                                                      <  2:715186/96‑1 (MQ=255)
                                                                                                                   cTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                            >  2:1377389/1‑150 (MQ=255)
                                                                                                                    tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTc                                                                                                                   >  2:1095780/1‑62 (MQ=255)
                                                                                                                    tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTc                                                                                                                   <  1:1095780/62‑1 (MQ=255)
                                                                                                                    tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                           >  4:44735/1‑150 (MQ=255)
                                                                                                                    tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                           >  1:106510/1‑150 (MQ=255)
                                                                                                                    tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                           >  2:952738/1‑150 (MQ=255)
                                                                                                                      aaCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATc                                                                                             <  1:1037603/82‑1 (MQ=255)
                                                                                                                      aaCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATc                                                                                             >  2:1037603/1‑82 (MQ=255)
                                                                                                                      aaCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgcta                                                              <  2:1404255/113‑3 (MQ=255)
                                                                                                                       aCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAt                        >  1:634360/1‑150 (MQ=255)
                                                                                                                         aTTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                      <  1:1282582/150‑1 (MQ=255)
                                                                                                                         aTTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                      >  4:40074/1‑150 (MQ=255)
                                                                                                                             cGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                  >  3:341311/1‑150 (MQ=255)
                                                                                                                              gCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAag                                                                      >  3:908103/1‑97 (MQ=255)
                                                                                                                                  aCGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc             <  4:229868/150‑1 (MQ=255)
                                                                                                                                    ggCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                >  4:553006/1‑81 (MQ=255)
                                                                                                                                    ggCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg           <  1:131156/150‑1 (MQ=255)
                                                                                                                                      cACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcc         >  2:952009/1‑150 (MQ=255)
                                                                                                                                         tAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                       <  4:109162/85‑1 (MQ=255)
                                                                                                                                         tAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                       >  3:109162/1‑85 (MQ=255)
                                                                                                                                         tAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcg      >  3:281588/1‑150 (MQ=255)
                                                                                                                                           aaaCCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcc         >  2:1492702/1‑145 (MQ=255)
                                                                                                                                             aCCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgcg  >  4:229167/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCG  >  NC_000913_3_bme_pgi/4233543‑4233833

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: