Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 2,665,522 A→G 6.8% T30A (ACC→GCC)  csiE → stationary phase inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi2,665,5220AG 6.8% 129.1 / 3.7 60T30A (ACC→GCC) csiEstationary phase inducible protein
Reads supporting (aligned to +/- strand):  ref base A (23/32);  new base G (2/2);  total (26/34)
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.29e-02

CGCCATTTCCTTGAGCAAACTTTAGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGAAGGTA  >  NC_000913_3_bme_pgi/2665374‑2665669
                                                                                                                                                    |                                                                                                                                                    
cGCCAGTTCCTTGAGCTAACTTCCGCTATTCTTATCAATTATGCTTATGGGAGCTCTACAGATGATGCCTACGGTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGa                                                                                                                                                      <  2:83959/149‑1 (MQ=255)
 gCCATTTCCTTGAGCAAACTTTAGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGCCGCTCTTTCAGCCGAGGTTGAc                                                                                                                                                     >  2:125013/1‑149 (MQ=255)
 gCCATTTCCGTGAGCAAACTTTTGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACTATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGAc                                                                                                                                                     <  2:766193/149‑1 (MQ=255)
         cTTGAGCAAACTTTAGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGccacca                                                                                                                                             >  1:595875/1‑149 (MQ=255)
         cTTGAGCAAACTTTAGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGGTGACCGccacca                                                                                                                                             >  2:630274/1‑149 (MQ=255)
         cTTGAGCAAACTTTAGCTATTCTTATCAATTATGCGTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCAACTGTCCTTTCGGCTCCCCAGCGCCACTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGCCCGacacca                                                                                                                                             >  2:217769/1‑149 (MQ=255)
           tGAGCAAACTTTAGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGTGTTGACCGCCACCATg                                                                                                                                           <  1:783568/149‑1 (MQ=255)
              gCAAACTTTAGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATCGCa                                                                                                                                        >  1:294447/1‑149 (MQ=255)
               cAAACTTTAGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCaa                                                                                                                                       >  1:367245/1‑149 (MQ=255)
                aaaCTTTAGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAAc                                                                                                                                      <  1:898084/149‑1 (MQ=255)
                       aGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGc                                                                                                                               >  1:572328/1‑149 (MQ=255)
                         cTATTCTTATCAATTATGCTTATGGGAGATCTACAGATTATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGa                                                                                                                             <  1:501024/149‑1 (MQ=255)
                         cTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCTGGGTTGACCGCCACCATGGCAACCTTCAGCGa                                                                                                                             >  1:750236/1‑149 (MQ=255)
                            ttCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCt                                                                                                                          >  1:779408/1‑149 (MQ=255)
                            ttCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCt                                                                                                                          <  1:738821/149‑1 (MQ=255)
                            ttCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCCCTCCCAGATCTTGCTCACGCCGCTTCAGCCGGCGTTGACCCCCCCCATGGCCACCCTCAGCgccct                                                                                                                          >  2:96900/1‑145 (MQ=255)
                                    aaTTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGgtg                                                                                                                  <  1:309297/149‑1 (MQ=255)
                                        aTGGTTATGGGTGATCCACAGAGCAGGCGTACATTTGCTTTAGCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTGGCTGACGCTCTTTCTGCCGGGGTTGACCGCCACCATGGTAATCTTCAGCGAGCTTAATGGTGTGga                                                                                                              <  2:72497/149‑1 (MQ=255)
                                        aTGCTTATGGGAGATCTACAGAGGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGAGCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGga                                                                                                              <  2:228923/149‑1 (MQ=255)
                                               tGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGTAACCTGCAGCGAGCTTAATGGTGTGgatgatgat                                                                                                      <  2:643460/149‑1 (MQ=255)
                                                gggAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGata                                                                                                     <  2:796411/149‑1 (MQ=255)
                                                        tACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTTCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCTAGt                                                                                             >  1:725058/1‑149 (MQ=255)
                                                          cAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATAATGATATTGCCAGTCt                                                                                           >  1:493231/1‑149 (MQ=255)
                                                             atgatgCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGTCAGTCTTgt                                                                                        >  1:740556/1‑148 (MQ=255)
                                                             atgatgCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGa                                                                                        <  1:505760/149‑1 (MQ=255)
                                                             aggctgCCTTGGCTGGCTGCAACATCTGGCCATTCGGGTCCCCAGCGCCGCTGACAGATCCTGCTGACGCTCTTTGAGCCGGGGTTGATCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTTGATGATGATATTGCCAGTCTTGa                                                                                        <  2:204890/145‑1 (MQ=255)
                                                                  gCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCa                                                                                   <  1:744155/149‑1 (MQ=255)
                                                                     tACGCTTGCTCCACCATCTGTCCTATCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGAGGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTGCAACGAGCTTAATGGGGCGGCTGATGCTGTTCCCAGTCTTGACATCAGCg                                                                                >  2:39923/1‑149 (MQ=255)
                                                                        gCTTGCTCCGCCAGTTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGaaa                                                                             <  1:876383/149‑1 (MQ=255)
                                                                             cTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGa                                                                        <  1:718524/149‑1 (MQ=255)
                                                                              tCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGGCCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGAc                                                                       >  1:116394/1‑149 (MQ=255)
                                                                              tCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGGCCGCCACCATGGCAACCTTCAGCGAGCTTAATGGAGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGAc                                                                       >  1:118736/1‑149 (MQ=255)
                                                                                    aTCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGgaga                                                                 <  1:601044/149‑1 (MQ=255)
                                                                                     tCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGGCCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGAt                                                                <  2:116394/149‑1 (MQ=255)
                                                                                     tCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGGCCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGAt                                                                <  2:118736/149‑1 (MQ=255)
                                                                                      cTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATc                                                               <  2:635567/149‑1 (MQ=255)
                                                                                      cTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATc                                                               <  2:294447/149‑1 (MQ=255)
                                                                                       tGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCTTCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATcc                                                              >  1:127380/1‑149 (MQ=255)
                                                                                         tCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCATCGAAACAGGACGGGAGATCCTg                                                            >  1:730946/1‑149 (MQ=255)
                                                                                            tttCGGCTCCCCAGCGCTGCTGCCCGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGTCACCATGGCAACCTTTAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTgcgc                                                         <  2:816779/149‑1 (MQ=255)
                                                                                                  cTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCTGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTatcat                                                   <  1:852052/149‑1 (MQ=255)
                                                                                                      ccAGTGCCGCTGCCAGATCTTGCTGACTCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCATCGAAACAGGACGGGAGATCCTGCGCTATCATCAAc                                               <  2:177357/149‑1 (MQ=255)
                                                                                                      ccAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAAc                                               >  1:7762/1‑149 (MQ=255)
                                                                                                      ccAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGAACGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAAc                                               >  2:4413/1‑149 (MQ=255)
                                                                                                         gggccgtTGCCAGATCTTGCTGATGCTCTTTTAGCCGGGGTTGACCGCCATCATGGCAACCTTCAGCGAGCTTAAGGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTca                                            <  1:415316/147‑1 (MQ=255)
                                                                                                            ccgcTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAACGAAACAGGAAGGGAGATCCTGCCCCATCCTCAACTcaccc                                         >  2:202779/1‑147 (MQ=255)
                                                                                                               cTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTga                                      <  1:892706/149‑1 (MQ=255)
                                                                                                               cTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTga                                      <  1:554786/149‑1 (MQ=255)
                                                                                                                tGCCAGATCTTGCTGACGCTCTTTCAGCCGGTTGTGATCTCCATCATTGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTgac                                     <  1:52486/149‑1 (MQ=255)
                                                                                                                     gATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTg                                <  1:468843/149‑1 (MQ=255)
                                                                                                                      aTCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGCGATCCTGCGCTATCATCAACTCACACTGACGACTgg                               >  1:24581/1‑149 (MQ=255)
                                                                                                                         ttGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTa                            <  1:832349/149‑1 (MQ=255)
                                                                                                                                 gCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAAATCGCACTGACGATTGGTTATGACGGTa                    >  1:556586/1‑149 (MQ=255)
                                                                                                                                     tttCAGCCGGGGTTCACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGATTATGACGGTCGCTa                >  1:778672/1‑149 (MQ=255)
                                                                                                                                           ccGGGGTTGACCCCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGt          <  2:742648/149‑1 (MQ=255)
                                                                                                                                               ggTTGACCGCCACCATGGCAACCTTCAGTGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGaa      <  2:308343/149‑1 (MQ=255)
                                                                                                                                                gTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGAAg     <  2:779408/149‑1 (MQ=255)
                                                                                                                                                gTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCAAGTCTTGATATCAGCGAAACAGGATGGGAGATCCTGCGCTATCATCAACTCAAACTGACGACTGGTTATGACGGTAGCTACAGGGTTGacg     >  2:838143/1‑147 (MQ=255)
                                                                                                                                                gTTGACCGCCAACAGGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGAAg     <  2:737159/149‑1 (MQ=255)
                                                                                                                                                 ttGACAGCCACCATGACAAAATTCAGCGAAATTAAAAGTATGAATGGATGTATAGCCAGGGTGAATAGCAGCGGGACAGGACGGGTGATCCTGCGCTATCATCAACTCACACTGACAACTGGTTATAACGGTAGCTAACGGGTTGaacg    >  2:81266/1‑147 (MQ=255)
                                                                                                                                                   gACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGAAGGTa  <  2:127380/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCCATTTCCTTGAGCAAACTTTAGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGATGATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGAAGGTA  >  NC_000913_3_bme_pgi/2665374‑2665669

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