Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,159,127 C→A A213E (GCG→GAG)  oxyR → oxidative and nitrosative stress transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,159,1270CA98.1% 534.5 / 1.7 156A213E (GCG→GAG) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base A (87/66);  total (90/66)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

Only 100 of 156 total aligned reads displayed.

CAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAG  >  NC_000913_3_bme_pgi/4158980‑4159272
                                                                                                                                                   |                                                                                                                                                 
cAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAGACTGCTGGTGCTGGAAGATGGTCACTGTTTGCGCGCTCAGGCGATGGGGTTCTGTTTTGAAGCCGGGGAgg                                                                                                                                                 >  4:912502/1‑150 (MQ=255)
cAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAgg                                                                                                                                                 <  4:803330/150‑1 (MQ=255)
  aTGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGTGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGAt                                                                                                                                               >  2:1283668/1‑150 (MQ=255)
      tGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAg                                                                                                                                           <  4:226120/150‑1 (MQ=255)
        cTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGAt                                                                                                                                         >  2:377674/1‑150 (MQ=255)
         tGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATa                                                                                                                                        >  1:384889/1‑150 (MQ=255)
          ggCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATac                                                                                                                                       >  2:1296408/1‑150 (MQ=255)
              atctatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGAt                                                                                                                                               >  1:293188/1‑138 (MQ=255)
              atctatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGAt                                                                                                                                               <  2:293188/138‑1 (MQ=255)
               tctatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGAGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACt                                                                                                                                  >  1:953608/1‑150 (MQ=255)
                 tatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTATTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTc                                                                                                                                >  3:325130/1‑150 (MQ=255)
                   tGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGATGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGACGCCGGGGAGGATGAAGATACAAACTTCcg                                                                                                                              >  4:958398/1‑150 (MQ=255)
                     aaGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCcgcg                                                                                                                            <  2:1473124/150‑1 (MQ=255)
                     aaGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCcgcg                                                                                                                            >  1:17063/1‑150 (MQ=255)
                         tCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGAc                                                                                                                        <  3:999924/150‑1 (MQ=255)
                          cACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATacac                                                                                                                                     >  2:1424146/1‑136 (MQ=255)
                           aCCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCa                                                                                                                      >  4:299798/1‑150 (MQ=255)
                           aCCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCa                                                                                                                      <  2:953608/150‑1 (MQ=255)
                            cccGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAg                                                                                                                     <  1:105006/150‑1 (MQ=255)
                            cccGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAg                                                                                                                     <  2:465908/150‑1 (MQ=255)
                                    cGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTgg                                                                                                                <  4:220730/147‑1 (MQ=255)
                                     gAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCt                                                                                                                  >  3:235645/1‑144 (MQ=255)
                                     gAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAAc                                                                                                            >  4:875534/1‑150 (MQ=255)
                                        ccGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAActct                                                                                                         >  2:288086/1‑150 (MQ=255)
                                          gcgAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTgc                                                                                                       >  2:228742/1‑150 (MQ=255)
                                              aTGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCaa                                                                                                   >  3:946194/1‑150 (MQ=255)
                                                gCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACa                                                                                                 >  2:149723/1‑150 (MQ=255)
                                                gCGTACCGATGGCCGATCCAGCTACGCCAAAACTACGGGTGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGACGATACCCACCTCCTCGCAACCAGCCTGGAAACTCTGCGCAACa                                                                                                 >  2:1658053/1‑150 (MQ=255)
                                                   tACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACAtgg                                                                                              >  4:571186/1‑150 (MQ=255)
                                                   tACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACAtgg                                                                                              >  3:178768/1‑150 (MQ=255)
                                                       gATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTggc                                                                                          >  3:861116/1‑150 (MQ=255)
                                                        aTGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACCTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTggcg                                                                                         >  2:235569/1‑150 (MQ=255)
                                                         tGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTTTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTggcgg                                                                                        <  1:72248/150‑1 (MQ=255)
                                                         tGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTggcgg                                                                                        <  3:918407/150‑1 (MQ=255)
                                                         tGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTggcgg                                                                                        <  2:53760/150‑1 (MQ=255)
                                                            ccGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAg                                                                                     >  1:1261711/1‑150 (MQ=255)
                                                            ccGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACCTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAg                                                                                     <  1:235569/150‑1 (MQ=255)
                                                             cGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCAGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAgg                                                                                    >  1:1014042/1‑150 (MQ=255)
                                                                    gCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCcg                                                                                                                              >  1:693663/1‑101 (MQ=255)
                                                                      aGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGAt                                                                           >  1:1423466/1‑150 (MQ=255)
                                                                      aGGGGAAAAACTGCTGATGCTGGAAGATGGTCACCGTTTGCGCGATCACGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGCGAt                                                                           >  2:792035/1‑150 (MQ=255)
                                                                         ggAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAgg                                                                                                                                                 >  3:915699/1‑77 (MQ=255)
                                                                         ggAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGa                                                                                                                         <  2:1600651/101‑1 (MQ=255)
                                                                         ggAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGa                                                                                                                         >  1:1600651/1‑101 (MQ=255)
                                                                         ggAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCAc                                                                        >  3:1006015/1‑150 (MQ=255)
                                                                          gAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACt                                                                       >  1:1638849/1‑150 (MQ=255)
                                                                          gAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACt                                                                       >  4:166494/1‑150 (MQ=255)
                                                                             aaaCTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTa                                                                    <  2:385661/150‑1 (MQ=255)
                                                                                 tgctgATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGc                                                                <  4:808122/150‑1 (MQ=255)
                                                                                  gctgATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGcc                                                               >  1:737297/1‑150 (MQ=255)
                                                                                   ctgATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCa                                                              <  3:166494/150‑1 (MQ=255)
                                                                                    tgATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCAc                                                                        >  3:376292/1‑139 (MQ=255)
                                                                                    tgATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAg                                                             <  3:603091/150‑1 (MQ=255)
                                                                                     gATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAgc                                                            <  3:875534/150‑1 (MQ=255)
                                                                                     gATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAgc                                                            <  1:423343/150‑1 (MQ=255)
                                                                                       tGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAgcgc                                                          >  4:819996/1‑150 (MQ=255)
                                                                                          tGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTa                                                                                  <  3:171029/123‑1 (MQ=255)
                                                                                          tGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTa                                                                                  >  4:171029/1‑123 (MQ=255)
                                                                                          tGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctgg                                                       >  2:626485/1‑150 (MQ=255)
                                                                                          tGGAAGATGGCCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCAGGGGAGGATGAAGATAAACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctgg                                                       >  2:1396150/1‑150 (MQ=255)
                                                                                           ggAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctgg                                                       <  4:357670/149‑1 (MQ=255)
                                                                                            gAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAAct                                                                                                           <  1:652886/96‑1 (MQ=255)
                                                                                            gAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAAct                                                                                                           >  2:652886/1‑96 (MQ=255)
                                                                                            gAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctggct                                                     <  1:1573687/150‑1 (MQ=255)
                                                                                            gAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctggct                                                     <  1:1364478/150‑1 (MQ=255)
                                                                                               gATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCtgtg                                                  >  1:1382718/1‑150 (MQ=255)
                                                                                                aTGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTgc                                                 <  3:569609/150‑1 (MQ=255)
                                                                                                 tGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTgcc                                                >  4:661491/1‑150 (MQ=255)
                                                                                                   gTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTAc                                                                   >  4:973483/1‑129 (MQ=255)
                                                                                                   gTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTAc                                                                   <  3:973483/129‑1 (MQ=255)
                                                                                                    tCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTgccgcc                                             >  2:120176/1‑150 (MQ=255)
                                                                                                     cACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTa                                                                    >  1:1465740/1‑126 (MQ=255)
                                                                                                       cTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGa                                          >  3:117938/1‑150 (MQ=255)
                                                                                                        tGTTTGCGCGATCAGGCAATGGGGTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCCGGCGGTGCCGCCGGAg                                         >  4:917573/1‑150 (MQ=255)
                                                                                                          tttGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAgcg                                       >  4:555971/1‑150 (MQ=255)
                                                                                                                  atcgggCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCAc                                                                        <  4:448783/105‑1 (MQ=255)
                                                                                                                  aTCGGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCAc                                                                        >  3:448783/1‑109 (MQ=255)
                                                                                                                  aTCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGTATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAAcgcg                               >  3:761891/1‑150 (MQ=255)
                                                                                                                  aTCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGCAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCACAcgcg                               >  2:1187200/1‑150 (MQ=255)
                                                                                                                    cAGGCAATGAGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctggc                                                      <  3:472994/125‑1 (MQ=255)
                                                                                                                    cAGGCAATGAGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctggc                                                      >  4:472994/1‑125 (MQ=255)
                                                                                                                     aGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAgcgc                                                          >  3:709299/1‑120 (MQ=255)
                                                                                                                     aGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAgcgc                                                          <  4:709299/120‑1 (MQ=255)
                                                                                                                       gCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCtgtg                                                  >  2:715230/1‑126 (MQ=255)
                                                                                                                       gCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATggg                          >  3:538221/1‑150 (MQ=255)
                                                                                                                       gCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATggg                          >  1:184467/1‑150 (MQ=255)
                                                                                                                          aTGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGgtt                       >  4:26/1‑150 (MQ=255)
                                                                                                                            gggTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGgttgt                     <  1:416320/150‑1 (MQ=255)
                                                                                                                             ggTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGgttgtt                    <  1:1296408/150‑1 (MQ=255)
                                                                                                                               tttCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTa                  >  2:529131/1‑150 (MQ=255)
                                                                                                                                 tCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATc                <  4:36994/150‑1 (MQ=255)
                                                                                                                                  cTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCt               <  3:107123/150‑1 (MQ=255)
                                                                                                                                     ttttGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGAt                             >  1:1637479/1‑133 (MQ=255)
                                                                                                                                     ttttGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGAt                             <  2:1637479/133‑1 (MQ=255)
                                                                                                                                      tttGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCg           <  3:661491/150‑1 (MQ=255)
                                                                                                                                         gAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGc        <  3:26/150‑1 (MQ=255)
                                                                                                                                            gCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGGGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGGTGTTTATCTGCCGTGCAtt     >  4:479917/1‑150 (MQ=255)
                                                                                                                                            gCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCAtt     >  4:946083/1‑150 (MQ=255)
                                                                                                                                               ggggAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACTTTTTTGCGGCAGGTGGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAg  <  1:1396150/150‑1 (MQ=255)
                                                                                                                                               ggggAGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGGTGTTTATCTGCCGTGCATTAAg  >  4:948992/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAG  >  NC_000913_3_bme_pgi/4158980‑4159272

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: