Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,159,121 C→A A211D (GCC→GAC)  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,1210CA92.2% 282.7 / 12.5 90A211D (GCC→GAC) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base A (40/43);  total (43/47)
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.12e-01

TGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGC  >  NC_000913_3_bme_pgi/4158975‑4159266
                                                                                                                                                  |                                                                                                                                                 
tgaGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACgg                                                                                                                                                <  3:465238/150‑1 (MQ=255)
 gaGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACggg                                                                                                                                               >  2:587402/1‑150 (MQ=255)
    ccAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCg                                                                                                                                            >  4:605529/1‑150 (MQ=255)
    ccAATGTTGCTCGCTATCTATGAAGATCACCCGTGGTCGAACCGCGAATGAGTACCTATGGCCGACCGGGCAGGGGAAAGTCTGCTGAGGGTGGTAGATGGTGACTGTTTGCGCGAGGAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCg                                                                                                                                            >  2:893699/1‑150 (MQ=255)
      aaTGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGa                                                                                                                                          <  1:196387/150‑1 (MQ=255)
        tGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATg                                                                                                                                        >  1:395674/1‑150 (MQ=255)
          ttGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGaa                                                                                                                                      <  2:140842/150‑1 (MQ=255)
          ttGCTGACTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGaa                                                                                                                                      <  1:809428/150‑1 (MQ=255)
           tGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAg                                                                                                                                     <  3:459837/150‑1 (MQ=255)
           tGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAg                                                                                                                                     <  2:137572/150‑1 (MQ=255)
            gCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGa                                                                                                                                    <  4:162403/150‑1 (MQ=255)
              tGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATa                                                                                                                                  <  3:649978/150‑1 (MQ=255)
                  tatctatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATacaca                                                                                                                              <  1:893699/150‑1 (MQ=255)
                   atctatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATacacac                                                                                                                             >  2:1059599/1‑150 (MQ=255)
                    tctatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACt                                                                                                                            <  2:236340/150‑1 (MQ=255)
                      tatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTc                                                                                                                          >  4:387964/1‑150 (MQ=255)
                        tGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCcg                                                                                                                        >  1:809752/1‑150 (MQ=255)
                             aTCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCGGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGa                                                                                                                   <  2:247033/150‑1 (MQ=255)
                               cACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGAcc                                                                                                                 >  1:450493/1‑150 (MQ=255)
                                  ccGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGa                                                                                                                   <  4:167055/145‑1 (MQ=255)
                                  ccGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGa                                                                                                                   >  3:167055/1‑145 (MQ=255)
                                  ccGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGc                                                                                                              >  2:563586/1‑150 (MQ=255)
                                    gTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGAGGATGAAGATACACACTTCCGCGCGACCAGCCt                                                                                                            <  2:684769/150‑1 (MQ=255)
                                    gTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAg                                                                                                                                     <  1:605787/125‑1 (MQ=255)
                                    gTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAg                                                                                                                                     >  2:605787/1‑125 (MQ=255)
                                     tGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTg                                                                                                           <  1:1051648/150‑1 (MQ=255)
                                       ggCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGa                                                                                                         >  2:864076/1‑150 (MQ=255)
                                                  aaTGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCa                                                                                              >  1:899002/1‑150 (MQ=255)
                                                       gTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACgggg                                                                                                                                              >  4:611599/1‑97 (MQ=255)
                                                       gTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACgggg                                                                                                                                              <  3:611599/97‑1 (MQ=255)
                                                         aCCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACAtggt                                                                                       >  2:326349/1‑150 (MQ=255)
                                                              tGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTggcgg                                                                                  <  4:127610/150‑1 (MQ=255)
                                                                gCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCa                                                                                <  3:387964/150‑1 (MQ=255)
                                                                 ccGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAg                                                                               <  1:326349/150‑1 (MQ=255)
                                                                  cGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCACGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAgg                                                                              >  4:379547/1‑150 (MQ=255)
                                                                       tGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCg                                                                         >  1:721264/1‑150 (MQ=255)
                                                                       tGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCg                                                                         >  3:425976/1‑150 (MQ=255)
                                                                        ggCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGaa                                                                                                                                      >  2:144427/1‑88 (MQ=255)
                                                                        ggCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGaa                                                                                                                                      <  1:144427/88‑1 (MQ=255)
                                                                         gCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGCAAct                                                                                                     >  1:797820/1‑120 (MQ=255)
                                                                         gCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAAct                                                                                                     <  2:797820/120‑1 (MQ=255)
                                                                          cAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGa                                                                      >  2:290702/1‑150 (MQ=255)
                                                                           aGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAGCATGGTGGCGGCAGGTAGCGGGAt                                                                     >  3:631758/1‑150 (MQ=255)
                                                                              ggAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCAc                                                                  >  4:579873/1‑150 (MQ=255)
                                                                              ggAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCAc                                                                  <  2:687675/150‑1 (MQ=255)
                                                                                   aaCTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGGTCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTAc                                                             >  4:387369/1‑150 (MQ=255)
                                                                                   aaCTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTAc                                                             <  2:9386/150‑1 (MQ=255)
                                                                                        ctgATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTgc                                                                                                 >  1:806896/1‑109 (MQ=255)
                                                                                        ctgATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTgc                                                                                                 <  2:806896/109‑1 (MQ=255)
                                                                                         tgATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAg                                                       >  3:85035/1‑150 (MQ=255)
                                                                                          gATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAgc                                                      >  1:79531/1‑150 (MQ=255)
                                                                                           aTGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTACCGGGATCACTTTACTGCCcgcg                                                     >  2:1006034/1‑150 (MQ=255)
                                                                                            tGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAgcgc                                                    <  3:429858/150‑1 (MQ=255)
                                                                                            tGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAgcgc                                                    >  1:596184/1‑150 (MQ=255)
                                                                                              cTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTGACTGCCAGCgctg                                                  >  2:885443/1‑150 (MQ=255)
                                                                                                 gAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGa                                                                                                                   <  3:420828/82‑1 (MQ=255)
                                                                                                 gAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGa                                                                                                                   >  4:420828/1‑82 (MQ=255)
                                                                                                  aaGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctggctg                                              >  4:326431/1‑150 (MQ=255)
                                                                                                     aTGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTgc                                           <  1:1059599/150‑1 (MQ=255)
                                                                                                      tGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTgcc                                          <  2:899002/150‑1 (MQ=255)
                                                                                                      tGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTgcc                                          <  3:56302/150‑1 (MQ=255)
                                                                                                         tCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACtt                                                                >  2:822008/1‑125 (MQ=255)
                                                                                                         tCACTGTTTGCGCGATCAGGCAATGGGTTACTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACtt                                                                <  1:822008/125‑1 (MQ=255)
                                                                                                             tGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAg                                   <  3:579873/150‑1 (MQ=255)
                                                                                                              gTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAgc                                  <  1:572613/150‑1 (MQ=255)
                                                                                                                ttGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAgcgc                                <  2:809752/150‑1 (MQ=255)
                                                                                                                ttGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAgcgc                                <  2:395674/150‑1 (MQ=255)
                                                                                                                   cgcgATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACt                                                            >  1:223359/1‑119 (MQ=255)
                                                                                                                   cgcgATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACt                                                            <  2:223359/119‑1 (MQ=255)
                                                                                                                      gATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAAcgc                          >  2:179963/1‑150 (MQ=255)
                                                                                                                        tCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGa                        <  1:563586/150‑1 (MQ=255)
                                                                                                                          aGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGCTg                      <  1:587402/150‑1 (MQ=255)
                                                                                                                          aGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATg                      >  1:815449/1‑150 (MQ=255)
                                                                                                                           ggCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCTAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATgg                     >  1:979769/1‑150 (MQ=255)
                                                                                                                            gCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCACGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATggg                    <  3:379547/150‑1 (MQ=255)
                                                                                                                             cAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGCTACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATgggg                   <  4:146512/150‑1 (MQ=255)
                                                                                                                             cAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATgggg                   <  2:242946/150‑1 (MQ=255)
                                                                                                                              aaTGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGgt                  <  2:458322/150‑1 (MQ=255)
                                                                                                                              aaTGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGgt                  >  1:182729/1‑150 (MQ=255)
                                                                                                                               aTGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGgtt                 <  1:290702/150‑1 (MQ=255)
                                                                                                                                    tttCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTa            >  1:1087117/1‑150 (MQ=255)
                                                                                                                                      tCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATc          <  2:5520/150‑1 (MQ=255)
                                                                                                                                        tGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGAGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTg        >  2:589573/1‑150 (MQ=255)
                                                                                                                                          ttttGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGcc      <  1:864076/150‑1 (MQ=255)
                                                                                                                                          ttttGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGcc      >  3:408758/1‑150 (MQ=255)
                                                                                                                                          ttttGAAGACGGGGCGGATGAAGATACACACGTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCGGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGcc      <  2:701415/150‑1 (MQ=255)
                                                                                                                                            ttGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGt    >  2:774817/1‑150 (MQ=255)
                                                                                                                                             tGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATg                      >  2:127747/1‑131 (MQ=255)
                                                                                                                                             tGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATg                      <  1:127747/131‑1 (MQ=255)
                                                                                                                                              gAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGc  <  1:589573/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGC  >  NC_000913_3_bme_pgi/4158975‑4159266

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: