Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,188,509 A→C T1054P (ACC→CCC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,188,5090AC96.1% 619.9 / 7.8 178T1054P (ACC→CCC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base C (84/87);  total (87/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

ACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGT  >  NC_000913_3_bme_pgi/4188411‑4188604
                                                                                                  |                                                                                               
aCCCGCACA‑CATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGccc                                                                                               <  2:1343617/100‑1 (MQ=255)
   cGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCgg                                                                                             <  2:3998209/100‑1 (MQ=255)
   cGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCg                                                                                              >  1:3028258/1‑99 (MQ=255)
     cacaCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGgtc                                                                                           <  2:535253/100‑1 (MQ=255)
        acCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGgtctgt                                                                                        >  2:2500719/1‑100 (MQ=255)
          cATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTAACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCAGgtctgtct                                                                                      >  1:4726760/1‑100 (MQ=255)
          cATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCGTTACTGACATGATCGCCGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGgtctgtct                                                                                      <  2:1500683/100‑1 (MQ=255)
          cATGCCGGTTATCACCGAAGTAAGCGGTCTTGTACGCTTTACTGACATGACCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCNCGgcctgtct                                                                                      >  2:314048/1‑100 (MQ=255)
              ccGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCg                                                                                   <  2:1938788/99‑1 (MQ=255)
                ggTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCtg                                                                                <  1:2467368/100‑1 (MQ=255)
                ggTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCCACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCg                                                                                   >  2:4301436/1‑97 (MQ=255)
                  ttATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCtggt                                                                              >  2:896891/1‑100 (MQ=255)
                    aTCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCtggtg                                                                             <  1:3919038/99‑1 (MQ=255)
                      cACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGtt                                                                          >  1:1046123/1‑100 (MQ=255)
                       aCCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTc                                                                         <  1:2959655/100‑1 (MQ=255)
                       aCCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTc                                                                         <  1:4066882/100‑1 (MQ=255)
                       aCCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTc                                                                         >  2:1306496/1‑100 (MQ=255)
                         cGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTg                                                                       <  2:4134047/100‑1 (MQ=255)
                         cGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTg                                                                       <  1:3395007/100‑1 (MQ=255)
                         cGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTg                                                                       <  1:2877755/100‑1 (MQ=255)
                         cGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTg                                                                       <  2:512673/100‑1 (MQ=255)
                         cGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTg                                                                       <  2:512704/100‑1 (MQ=255)
                          gAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTgg                                                                      <  2:1384475/100‑1 (MQ=255)
                             gTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGAtt                                                                   >  2:1441204/1‑100 (MQ=255)
                             gTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGAtt                                                                   <  1:4305395/100‑1 (MQ=255)
                              tAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTc                                                                  >  1:37700/1‑100 (MQ=255)
                               aaGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTcc                                                                 >  2:230244/1‑100 (MQ=255)
                                aGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGAt                                                                    <  2:350903/96‑1 (MQ=255)
                                aGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTcc                                                                 <  1:4105177/99‑1 (MQ=255)
                                  cGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCg                                                                >  1:1783735/1‑98 (MQ=255)
                                  cGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCg                                                                >  1:2198375/1‑98 (MQ=255)
                                    gTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGGCCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAg                                                             >  1:3958559/1‑99 (MQ=255)
                                    gTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGa                                                            <  2:2452947/100‑1 (MQ=255)
                                     ttttGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGaa                                                           >  1:841127/1‑100 (MQ=255)
                                      tttGTACGCTTTACTGGCATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAAc                                                          >  1:2679264/1‑100 (MQ=255)
                                       ttGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAAc                                                          <  2:271145/99‑1 (MQ=255)
                                         gTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTa                                                       >  2:2067508/1‑100 (MQ=255)
                                         gTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTa                                                       >  1:2745406/1‑100 (MQ=255)
                                            cGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCg                                                    >  1:3858228/1‑100 (MQ=255)
                                              cTTTACTGACATGATCGACGGTCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCa                                                  <  2:702349/100‑1 (MQ=255)
                                              cTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCg                                                    >  2:2033172/1‑98 (MQ=255)
                                                 tACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGACTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAggt                                               >  1:4621851/1‑100 (MQ=255)
                                                  aCTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAggt                                               <  1:1682839/99‑1 (MQ=255)
                                                   cTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAggtgg                                             >  2:1075005/1‑100 (MQ=255)
                                                   cTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAggtgg                                             >  1:1080246/1‑100 (MQ=255)
                                                     gACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAggtggt                                            >  1:4566903/1‑99 (MQ=255)
                                                     gACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAggtggt                                            >  2:168206/1‑99 (MQ=255)
                                                     gACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCCGGTGGTTCTGGATTCCGCAGAACGTACCGCAggtggt                                            >  2:3556516/1‑99 (MQ=255)
                                                       cATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTa                                           >  2:3183780/1‑98 (MQ=255)
                                                        aTGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAg                                        <  1:92687/100‑1 (MQ=255)
                                                         tGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTa                                           >  2:4695753/1‑96 (MQ=255)
                                                         tGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTa                                           >  1:695686/1‑96 (MQ=255)
                                                          gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTa                                           >  2:3389403/1‑95 (MQ=255)
                                                          gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTa                                           >  1:592522/1‑95 (MQ=255)
                                                           aTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTCCCGCAGGTGGTa                                           >  2:3881742/1‑94 (MQ=255)
                                                           aTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGa                                       <  2:1126329/98‑1 (MQ=255)
                                                           aTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATc                                     <  1:649970/100‑1 (MQ=255)
                                                            tCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATc                                     >  1:34434/1‑99 (MQ=255)
                                                              gACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTg                                   >  2:3347002/1‑99 (MQ=255)
                                                              gACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTg                                   <  1:2935173/99‑1 (MQ=255)
                                                              gACGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGc                                  >  2:4223355/1‑100 (MQ=255)
                                                                cGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGt                                <  2:2578692/100‑1 (MQ=255)
                                                                cGGCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGt                                >  2:2998034/1‑100 (MQ=255)
                                                                 ggCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCg                                 <  1:3677276/98‑1 (MQ=255)
                                                                  gCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTcc                              >  1:3624395/1‑100 (MQ=255)
                                                                  gCCAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTcc                              <  2:1253825/100‑1 (MQ=255)
                                                                   ccAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCg                             <  2:2272145/100‑1 (MQ=255)
                                                                   ccAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCg                             >  1:3165640/1‑100 (MQ=255)
                                                                    cggcccATTACGCGACAGACCGACGAACTGCCCGGTCTGTCTTCGCTAGTGCTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGCTCTGCGTCCgg                            <  1:734891/96‑1 (MQ=255)
                                                                     aGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGc                           <  2:3271596/100‑1 (MQ=255)
                                                                     aGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGc                           <  2:837338/100‑1 (MQ=255)
                                                                      gACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCa                          >  2:4449655/1‑100 (MQ=255)
                                                                       aCCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCAc                         <  1:3649546/100‑1 (MQ=255)
                                                                       aCCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCAc                         <  1:4236554/100‑1 (MQ=255)
                                                                         cATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACt                        >  1:2082655/1‑99 (MQ=255)
                                                                          aTTACGCGTCAGACCGACGTACTGCCCGGTCTGCGTTCGCTGGTGGTTCTGGATTGCGCAGAACCTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGa                      <  2:3953748/100‑1 (MQ=255)
                                                                           ttACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGaa                     <  1:2421846/100‑1 (MQ=255)
                                                                           ttACGCGTCAGACCGACGAACTCCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTCACGATCTGCGTCCGGCACTGa                      >  2:1649350/1‑99 (MQ=255)
                                                                             aCGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTCAAGATCTGCGTCCGGCa                          >  2:4098362/1‑93 (MQ=255)
                                                                              cgcgTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGaaaa                   <  1:4010453/99‑1 (MQ=255)
                                                                              cgcgTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAAt                  >  1:1890536/1‑100 (MQ=255)
                                                                              cgcgTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAAt                  <  2:1126638/100‑1 (MQ=255)
                                                                               gcgTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATc                 <  2:1160745/100‑1 (MQ=255)
                                                                                cgTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTCGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCg                <  1:3582109/100‑1 (MQ=255)
                                                                                cgTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATc                 <  2:474129/99‑1 (MQ=255)
                                                                                  tCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGtt              <  1:1771818/100‑1 (MQ=255)
                                                                                    aGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTg             <  1:744594/99‑1 (MQ=255)
                                                                                    aGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTg             <  1:3555556/99‑1 (MQ=255)
                                                                                    aGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTg             <  2:2281406/99‑1 (MQ=255)
                                                                                    aGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTg             <  2:1088193/99‑1 (MQ=255)
                                                                                    aGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGa            >  1:1271130/1‑100 (MQ=255)
                                                                                     gACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGAt           >  2:1543173/1‑100 (MQ=255)
                                                                                       ccGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGAt           >  2:3758816/1‑98 (MQ=255)
                                                                                       ccGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATg          >  2:3660985/1‑99 (MQ=255)
                                                                                        cgacgaACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCt        >  2:404050/1‑100 (MQ=255)
                                                                                           cgaACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAg     <  2:2354726/100‑1 (MQ=255)
                                                                                             aaCTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAggg   <  1:3761089/100‑1 (MQ=255)
                                                                                             aaCTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAggg   <  1:1605523/100‑1 (MQ=255)
                                                                                             aaCTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGGACCGCAGGTGGTAAAGAGCTGCGTCCGGCACTGAAAATCGTTGATGCTCAggg   <  2:4306280/100‑1 (MQ=255)
                                                                                               cTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGAGTGCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGt  <  2:3434174/99‑1 (MQ=255)
                                                                                                  |                                                                                               
ACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGT  >  NC_000913_3_bme_pgi/4188411‑4188604

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