Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,159,081 C→T H198Y (CAC→TAC)  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,0810CT98.1% 359.7 / ‑4.3 103H198Y (CAC→TAC) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (49/52);  total (49/54)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

AAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCC  >  NC_000913_3_bme_pgi/4158937‑4159226
                                                                                                                                                    |                                                                                                                                                 
aaaGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGt                                                                                                                                              <  4:541110/150‑1 (MQ=255)
         gAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGAt                                                                                                                                     <  3:317445/150‑1 (MQ=255)
          aaGCATTCATTGAAGTTCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGAGGGCCGATCTGGCAGAGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATc                                                                                                                                    <  2:342953/150‑1 (MQ=255)
           aGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCa                                                                                                                                   >  4:624088/1‑150 (MQ=255)
               ttcattGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCa                                                                                                                               >  1:1024306/1‑150 (MQ=255)
               ttcattGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCa                                                                                                                               >  1:358648/1‑150 (MQ=255)
                 cattGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAAt                                                                                                                             <  2:358648/150‑1 (MQ=255)
                  attGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATg                                                                                                                            >  1:900258/1‑150 (MQ=255)
                  attGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATg                                                                                                                            >  1:744431/1‑150 (MQ=255)
                  attGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATg                                                                                                                            >  1:497763/1‑150 (MQ=255)
                     gAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGt                                                                                                                         <  2:497763/150‑1 (MQ=255)
                      aaGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGtt                                                                                                                        <  4:690453/150‑1 (MQ=255)
                        gTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTc                                                                                                                      <  2:1340412/150‑1 (MQ=255)
                            cGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGtt                                                                                                                  <  3:717248/150‑1 (MQ=255)
                             gttgttTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGttt                                                                                                                 <  3:481647/150‑1 (MQ=255)
                                  ttGATGAGCCAATGTTGCTGGCTATCTATCTATTAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGtttt                                                                                                                <  4:742523/150‑1 (MQ=255)
                                   tgatgaGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGc                                                                                                           >  1:1438010/1‑150 (MQ=255)
                                    gatgaGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGcc                                                                                                          >  4:72476/1‑150 (MQ=255)
                                     atgaGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCg                                                                                                         >  4:81453/1‑150 (MQ=255)
                                         gCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCgggg                                                                                                      <  3:24391/149‑1 (MQ=255)
                                         gCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCgggg                                                                                                      >  4:24391/1‑149 (MQ=255)
                                         gCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTCTGAAGCCGGGGc                                                                                                     <  2:744431/150‑1 (MQ=255)
                                            aaTGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGa                                                                                                  >  2:983501/1‑150 (MQ=255)
                                            aaTGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGa                                                                                                  <  3:624088/150‑1 (MQ=255)
                                            aaTGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGa                                                                                                  >  4:187889/1‑150 (MQ=255)
                                               gTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGtt                                                                                                                                             <  3:427394/104‑1 (MQ=255)
                                               gTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGtt                                                                                                                                             >  4:427394/1‑104 (MQ=255)
                                                 tGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAg                                                                                             <  3:222875/150‑1 (MQ=255)
                                                  gCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGa                                                                                            <  1:983613/150‑1 (MQ=255)
                                                    tGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGa                                                                                            >  2:1163248/1‑148 (MQ=255)
                                                    tGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGa                                                                                            <  1:1163248/148‑1 (MQ=255)
                                                    tGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATa                                                                                          <  3:81453/150‑1 (MQ=255)
                                                      g‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATaca                                                                                        <  2:1061018/150‑1 (MQ=255)
                                                            tatctatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATacacac                                                                                     >  1:1589655/1‑150 (MQ=255)
                                                             atctatGAAGAGCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATacacac                                                                                     <  2:694574/150‑1 (MQ=255)
                                                                tatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATaca                                                                                        <  3:612771/144‑1 (MQ=255)
                                                                tatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATaca                                                                                        >  4:612771/1‑144 (MQ=255)
                                                                 atGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGCAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTcc                                                                                 >  1:1386891/1‑150 (MQ=255)
                                                                  tGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGaa                                                                                                             <  1:1302096/121‑1 (MQ=255)
                                                                  tGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGaa                                                                                                             >  2:1302096/1‑121 (MQ=255)
                                                                  tGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGa                                                                                               <  2:1433347/135‑1 (MQ=255)
                                                                  tGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGa                                                                                               >  1:1433347/1‑135 (MQ=255)
                                                                  tGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATa                                                                                          >  4:393315/1‑140 (MQ=255)
                                                                  tGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATa                                                                                          <  3:393315/140‑1 (MQ=255)
                                                                   gAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCcgc                                                                               <  3:717830/150‑1 (MQ=255)
                                                                       agcACCCGTGGGGAAACCTCGAATGCGTACCGATGGCCGATCTAGCAGGGGAAAAACTGCTGATGCGGTAAGATGGTCACGGTGTGCGCGATCAGGCAATGGGTTGCTATTTTGAAGCCGAGGCAGATGAAGATACACACTTCTGCGCGa                                                                           <  2:1255366/148‑1 (MQ=255)
                                                                        tCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGAc                                                                          >  4:833277/1‑150 (MQ=255)
                                                                           cccGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAg                                                                       <  1:255975/150‑1 (MQ=255)
                                                                           cccGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAg                                                                       <  1:983501/150‑1 (MQ=255)
                                                                            ccGTGGGCGAACCGCTAATGCGTACCAATTTCCTATCTGTCAGATAAAAAATTTCTTATTCTTGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGc                                                                      <  3:581432/150‑1 (MQ=255)
                                                                            ccGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCggg                                                                                                       <  4:139597/117‑1 (MQ=255)
                                                                            ccGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCggg                                                                                                       >  3:139597/1‑117 (MQ=255)
                                                                              gTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCt                                                                    <  1:1097013/150‑1 (MQ=255)
                                                                              gTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCt                                                                    <  2:1024306/150‑1 (MQ=255)
                                                                              gTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCGGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGGTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTTCGCGCGACCAGCCt                                                                    >  4:280004/1‑150 (MQ=255)
                                                                                gggCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTgg                                                                  <  1:959387/150‑1 (MQ=255)
                                                                                   cGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACAGACTTCCGCGCGACCAGCCTGGaaa                                                               <  4:519058/150‑1 (MQ=255)
                                                                                     aaCCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGAt                                                                                           <  3:18695/120‑1 (MQ=255)
                                                                                     aaCCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGAt                                                                                           >  4:18695/1‑120 (MQ=255)
                                                                                      aCCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATaca                                                                                        <  1:313314/122‑1 (MQ=255)
                                                                                      aCCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATaca                                                                                        >  2:313314/1‑122 (MQ=255)
                                                                                       ccGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGc                                                                                                     >  1:282019/1‑108 (MQ=255)
                                                                                       ccGCGAATGCGTACAGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGc                                                                                                     <  2:282019/108‑1 (MQ=255)
                                                                                        cgcgAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTg                                                          >  4:320772/1‑150 (MQ=255)
                                                                                         gcgAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTgc                                                         >  1:1299672/1‑150 (MQ=255)
                                                                                         gcgAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTgc                                                         <  3:72476/150‑1 (MQ=255)
                                                                                          cgAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTgcgc                                                                                                                                        >  2:133982/1‑70 (MQ=255)
                                                                                          cgAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTgcgc                                                                                                                                        <  1:133982/70‑1 (MQ=255)
                                                                                           gAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTgcgc                                                       <  3:320772/150‑1 (MQ=255)
                                                                                           gAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTgcgc                                                       <  1:1185120/150‑1 (MQ=255)
                                                                                                  tACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACAtgg                                                >  1:441147/1‑150 (MQ=255)
                                                                                                      gATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGa                                                                 >  2:527519/1‑129 (MQ=255)
                                                                                                      gATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGa                                                                 <  1:527519/129‑1 (MQ=255)
                                                                                                        tGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTggcgg                                          <  2:223041/150‑1 (MQ=255)
                                                                                                        tGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTggcgg                                          >  2:1462066/1‑150 (MQ=255)
                                                                                                         ggCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTggcggc                                         >  3:423557/1‑150 (MQ=255)
                                                                                                          gCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCa                                        >  3:323674/1‑150 (MQ=255)
                                                                                                          gCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCa                                        <  4:469467/150‑1 (MQ=255)
                                                                                                           ccGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGt                                                                                                                   >  2:96422/1‑74 (MQ=255)
                                                                                                           ccGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGt                                                                                                                   <  1:96422/74‑1 (MQ=255)
                                                                                                           ccGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAg                                       >  1:836337/1‑150 (MQ=255)
                                                                                                            cGATCTGGCAGGGGAAAACCTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAgg                                      <  3:187889/150‑1 (MQ=255)
                                                                                                              aTCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTa                                    <  2:900258/150‑1 (MQ=255)
                                                                                                               tCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAg                                   >  3:679102/1‑150 (MQ=255)
                                                                                                                   gCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCggg                               >  3:438112/1‑150 (MQ=255)
                                                                                                                   gCAGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCggg                               >  2:98525/1‑150 (MQ=255)
                                                                                                                     aGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGAt                             >  4:399442/1‑150 (MQ=255)
                                                                                                                     aGGGGAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGAt                             >  3:802324/1‑150 (MQ=255)
                                                                                                                        ggAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGACGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCAc                          >  3:949840/1‑150 (MQ=255)
                                                                                                                        ggAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGc                                  >  3:688665/1‑142 (MQ=255)
                                                                                                                        ggAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGc                                  <  4:688665/142‑1 (MQ=255)
                                                                                                                         gAAAAACTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACt                         <  4:423557/150‑1 (MQ=255)
                                                                                                                          aaaaaCTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACtt                        >  3:33758/1‑150 (MQ=255)
                                                                                                                          aaaaaCTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACtt                        >  2:283799/1‑150 (MQ=255)
                                                                                                                           aaaaCTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACttt                       >  2:1313203/1‑150 (MQ=255)
                                                                                                                           aaaaCTGCTGATGCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACttt                       >  2:698677/1‑150 (MQ=255)
                                                                                                                                      ggCTGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAgcgc            <  4:258188/149‑1 (MQ=255)
                                                                                                                                         tGGAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctgg         >  3:483634/1‑150 (MQ=255)
                                                                                                                                          ggAAGATGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctggc        <  1:914499/150‑1 (MQ=255)
                                                                                                                                                tGGTTACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTgcc  <  3:399442/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGG‑‑‑‑CTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCC  >  NC_000913_3_bme_pgi/4158937‑4159226

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