Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,772 G→T intergenic (‑539/‑2135) 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,7720GT99.4% 640.8 / NA 176intergenic (‑539/‑2135)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (98/77);  total (99/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

TAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAA  >  NC_000913_3_bme_pgi/4233626‑4233918
                                                                                                                                                  |                                                                                                                                                   
taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacg                                                                                                                                                  >  2:1062177/1‑150 (MQ=255)
 aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacga                                                                                                                                                 >  4:682640/1‑150 (MQ=255)
 aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacga                                                                                                                                                 >  2:869797/1‑150 (MQ=255)
   taaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAg                                                                                                                                               <  3:849911/150‑1 (MQ=255)
    aaaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGc                                                                                                                                              <  4:535497/150‑1 (MQ=255)
     aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACCCTACGAAGCt                                                                                                                                             >  2:1503653/1‑150 (MQ=255)
      aTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCtt                                                                                                                                            <  2:1402608/150‑1 (MQ=255)
         cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                                                                                         <  4:214875/150‑1 (MQ=255)
         cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                                                                                         <  2:289915/150‑1 (MQ=255)
          ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                                                                                        >  1:1081577/1‑150 (MQ=255)
          ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                                                                                        >  4:582171/1‑150 (MQ=255)
            cACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCa                                                                                                                                      >  1:1144138/1‑150 (MQ=255)
                   cTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGc                                                                                                                                              <  4:228000/135‑1 (MQ=255)
                   cTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGc                                                                                                                                              >  3:228000/1‑135 (MQ=255)
                        ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                                                                                                          >  3:893378/1‑150 (MQ=255)
                          tttACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGc                                                                                                                         >  3:649581/1‑149 (MQ=255)
                          tttACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCt                                                                                                                        >  2:421419/1‑150 (MQ=255)
                                cTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                                                                                                  <  1:365641/150‑1 (MQ=255)
                                 tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                                                                                                 >  1:914705/1‑150 (MQ=255)
                                  cAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGc                                                                                                                <  2:1609745/150‑1 (MQ=255)
                                   aaCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCa                                                                                                               <  3:313190/150‑1 (MQ=255)
                                   aaCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCa                                                                                                               <  4:631980/150‑1 (MQ=255)
                                     cATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGAAAAGATAAAAAAGCAtt                                                                                                             <  2:1081577/150‑1 (MQ=255)
                                       ttACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTgg                                                                                                           <  1:1062177/150‑1 (MQ=255)
                                        tACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTggg                                                                                                          <  1:869797/150‑1 (MQ=255)
                                           gCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCTCAAAGCGGAAAAAGCATTGGGCgg                                                                                                       >  4:105151/1‑150 (MQ=255)
                                           gCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCgg                                                                                                       >  4:963901/1‑150 (MQ=255)
                                           gCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCgg                                                                                                       >  1:371871/1‑150 (MQ=255)
                                               aCGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGATCAAAAAGCTGAAAAAGCATTGGGCGGAATc                                                                                                   <  2:1785333/150‑1 (MQ=255)
                                                 ggCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                                                                                                 >  3:96518/1‑150 (MQ=255)
                                                 ggCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                                                                                                 >  2:1156015/1‑150 (MQ=255)
                                                  gCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGTGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGc                                                                                                >  2:764059/1‑150 (MQ=255)
                                                     cTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                                                                                         <  1:1160311/106‑1 (MQ=255)
                                                     cTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                                                                                         >  2:1160311/1‑106 (MQ=255)
                                                     cTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                                                                                             >  1:477177/1‑150 (MQ=255)
                                                      tAAAACCATTACATTTTTTTTTTTTTTTCTCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAATAGACTCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcg                                                                                            >  2:92296/1‑150 (MQ=255)
                                                      tAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcg                                                                                            >  3:811107/1‑150 (MQ=255)
                                                      tAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcg                                                                                            >  1:794668/1‑150 (MQ=255)
                                                        aaaCCATTACATTTTTCTTTTACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAATTCAGAAGAGTATTGCTAATGTCAATAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcg                                                                                          >  4:1094105/1‑150 (MQ=255)
                                                                  aTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgc                                                                               >  4:434258/1‑150 (MQ=255)
                                                                        cTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCa                                                                         >  2:704728/1‑150 (MQ=255)
                                                                         tgtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAAt                                                                                                    <  2:40978/123‑1 (MQ=255)
                                                                         tgtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAt                                                                        >  2:779042/1‑150 (MQ=255)
                                                                          gtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAtt                                                                       >  1:379315/1‑150 (MQ=255)
                                                                           tgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTc                                                                      >  4:648899/1‑150 (MQ=255)
                                                                           tgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTc                                                                      <  3:963901/150‑1 (MQ=255)
                                                                                ggCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcgc                                                                 <  2:477161/150‑1 (MQ=255)
                                                                                 gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                                                                                                          <  1:581563/93‑1 (MQ=255)
                                                                                 gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                                                                                                          >  2:581563/1‑93 (MQ=255)
                                                                                 gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                                                >  2:285740/1‑150 (MQ=255)
                                                                                   gcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGa                                                              <  4:893378/150‑1 (MQ=255)
                                                                                            ttCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                                                                                        >  1:823058/1‑68 (MQ=255)
                                                                                            ttCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                                                                                        <  2:823058/68‑1 (MQ=255)
                                                                                            ttCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGa                                                                                                                      >  3:121500/1‑86 (MQ=255)
                                                                                             tCCAAAGAAACAATTCTCAAAATCAGAAGAGTATTCCTCATGGCAAACCACGCTACGACGCTTAAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggccgg                                                    >  4:352692/1‑150 (MQ=255)
                                                                                                aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCANGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTc                                                 >  1:243500/1‑150 (MQ=255)
                                                                                                aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggc                                                       <  4:623356/144‑1 (MQ=255)
                                                                                                aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggc                                                       >  3:623356/1‑144 (MQ=255)
                                                                                                aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGtca                                                 >  3:499937/1‑148 (MQ=255)
                                                                                                aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTc                                                 >  4:1026839/1‑150 (MQ=255)
                                                                                                    tCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                                             >  1:135994/1‑150 (MQ=255)
                                                                                                    tCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                                             <  3:187239/150‑1 (MQ=255)
                                                                                                     cacaATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCTCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCg                                            >  1:994897/1‑150 (MQ=255)
                                                                                                      acaATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCgg                                           >  4:126753/1‑149 (MQ=255)
                                                                                                      acaATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGa                                           >  1:1492350/1‑150 (MQ=255)
                                                                                                      acaATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGa                                           >  4:287089/1‑150 (MQ=255)
                                                                                                        aaTTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGata                                         >  1:697741/1‑150 (MQ=255)
                                                                                                          ttCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATAtt                                       >  3:93383/1‑150 (MQ=255)
                                                                                                          ttCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATAtt                                       >  3:301592/1‑150 (MQ=255)
                                                                                                            ctcAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCg                                     >  4:447179/1‑150 (MQ=255)
                                                                                                               aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGACGCAACGGTCGCAAAGCTGAAAAAGCCTTGCGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATAGTCGctc                                  >  2:1182741/1‑150 (MQ=255)
                                                                                                               aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctc                                  >  3:34370/1‑150 (MQ=255)
                                                                                                               aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctc                                  >  4:266676/1‑150 (MQ=255)
                                                                                                                   tCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCATATCGCGCgggg                                                                         >  2:63438/1‑103 (MQ=255)
                                                                                                                    cAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGCCAGGTTCGGGCGATATTCGCTCTGCCt                             <  4:783603/150‑1 (MQ=255)
                                                                                                                     agaagaGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTg                            >  1:1304598/1‑150 (MQ=255)
                                                                                                                      gaagaGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTgg                           >  2:1394278/1‑150 (MQ=255)
                                                                                                                       aagaGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGCATATCGCGCggg                                                                          <  1:1229489/102‑2 (MQ=255)
                                                                                                                       aagaGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCt                             <  1:978814/147‑1 (MQ=255)
                                                                                                                       aagaGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCt                             >  2:978814/1‑147 (MQ=255)
                                                                                                                       aagaGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGga                          >  3:720767/1‑150 (MQ=255)
                                                                                                                           gTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGc                                                                                                >  1:462729/1‑77 (MQ=255)
                                                                                                                              ttGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtctt                   <  2:552624/150‑1 (MQ=255)
                                                                                                                              ttGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtctt                   <  3:737277/150‑1 (MQ=255)
                                                                                                                               tGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                                                                                                        >  1:1735024/1‑65 (MQ=255)
                                                                                                                               tGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                                                                                                        <  2:1735024/65‑1 (MQ=255)
                                                                                                                               tGCTAACGGACAGAGACGCTACGACGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATCTTCGCTCTGCCTGGATGAtcttc                  <  2:1339481/150‑1 (MQ=255)
                                                                                                                                 cTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCCGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgc                                                                                 >  1:100995/1‑85 (MQ=255)
                                                                                                                                 cTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCCGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgc                                                                                 <  2:100995/85‑1 (MQ=255)
                                                                                                                                 cTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGTCGCCTCGCTAGATGCGCGCTGCATGCAGCGCCGACCCCGGCCTGTTCGGGCGATAGTCGCTCTGCCTGGATTAtcttcta                >  1:1339481/1‑149 (MQ=255)
                                                                                                                                  tAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                                              >  4:23455/1‑119 (MQ=255)
                                                                                                                                  tAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGt               <  3:1026839/150‑1 (MQ=255)
                                                                                                                                  tAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGt               >  2:1144672/1‑150 (MQ=255)
                                                                                                                                  tAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGt               >  4:161371/1‑150 (MQ=255)
                                                                                                                                      ggCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATGGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGa           <  3:687321/150‑1 (MQ=255)
                                                                                                                                      ggCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGa           >  1:158098/1‑150 (MQ=255)
                                                                                                                                        cAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACt         <  3:253125/150‑1 (MQ=255)
                                                                                                                                         aaGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTg        >  2:614349/1‑150 (MQ=255)
                                                                                                                                              aCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGa   <  2:135994/150‑1 (MQ=255)
                                                                                                                                               cgctACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGaa  <  1:1394278/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAA  >  NC_000913_3_bme_pgi/4233626‑4233918

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: