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

CGCCATTTCCTTGAGCAAACTTTAGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGAAGGTA  >  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)
                                               tGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGTAACCTGCAGCGAGCTTAATGGTGTGgat‑gatgat                                                                                                      <  2:643460/149‑1 (MQ=255)
                                                gggAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGata                                                                                                     <  2:796411/149‑1 (MQ=255)
                                                        tACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTTCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCTAGt                                                                                             >  1:725058/1‑149 (MQ=255)
                                                          cAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑AATGATATTGCCAGTCt                                                                                           >  1:493231/1‑149 (MQ=255)
                                                             atgatgCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGTCAGTCTTgt                                                                                        >  1:740556/1‑148 (MQ=255)
                                                             atgatgCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGa                                                                                        <  1:505760/149‑1 (MQ=255)
                                                             aggctgCCTTGGCTGGCTGCAACATCTGGCCATTCGGGTCCCCAGCGCCGCTGACAGATCCTGCTGACGCTCTTTGAGCCGGGGTTGATCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTTGAT‑GATGATATTGCCAGTCTTGa                                                                                        <  2:204890/145‑1 (MQ=255)
                                                                  gCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCa                                                                                   <  1:744155/149‑1 (MQ=255)
                                                                     tACGCTTGCTCCACCATCTGTCCTATCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGAGGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTGCAACGAGCTTAATGGGGCGGCT‑GATGCTGTTCCCAGTCTTGACATCAGCg                                                                                >  2:39923/1‑149 (MQ=255)
                                                                        gCTTGCTCCGCCAGTTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGaaa                                                                             <  1:876383/149‑1 (MQ=255)
                                                                             cTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGa                                                                        <  1:718524/149‑1 (MQ=255)
                                                                              tCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGGCCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGAc                                                                       >  1:116394/1‑149 (MQ=255)
                                                                              tCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGGCCGCCACCATGGCAACCTTCAGCGAGCTTAATGGAGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGAc                                                                       >  1:118736/1‑149 (MQ=255)
                                                                                    aTCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGgaga                                                                 <  1:601044/149‑1 (MQ=255)
                                                                                     tCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGGCCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGAt                                                                <  2:116394/149‑1 (MQ=255)
                                                                                     tCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGGCCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGAt                                                                <  2:118736/149‑1 (MQ=255)
                                                                                      cTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATc                                                               <  2:635567/149‑1 (MQ=255)
                                                                                      cTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATc                                                               <  2:294447/149‑1 (MQ=255)
                                                                                       tGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCTTCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATcc                                                              >  1:127380/1‑149 (MQ=255)
                                                                                         tCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCATCGAAACAGGACGGGAGATCCTg                                                            >  1:730946/1‑149 (MQ=255)
                                                                                            tttCGGCTCCCCAGCGCTGCTGCCCGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGTCACCATGGCAACCTTTAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTgcgc                                                         <  2:816779/149‑1 (MQ=255)
                                                                                                  cTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCTGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTatcat                                                   <  1:852052/149‑1 (MQ=255)
                                                                                                      ccAGTGCCGCTGCCAGATCTTGCTGACTCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCATCGAAACAGGACGGGAGATCCTGCGCTATCATCAAc                                               <  2:177357/149‑1 (MQ=255)
                                                                                                      ccAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAAc                                               >  1:7762/1‑149 (MQ=255)
                                                                                                      ccAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGAACGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAAc                                               >  2:4413/1‑149 (MQ=255)
                                                                                                         gggccgtTGCCAGATCTTGCTGATGCTCTTTTAGCCGGGGTTGACCGCCATCATGGCAACCTTCAGCGAGCTTAAGGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTca                                            <  1:415316/147‑1 (MQ=255)
                                                                                                            ccgcTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAACGAAACAGGAAGGGAGATCCTGCCCCATCCTCAACTcaccc                                         >  2:202779/1‑147 (MQ=255)
                                                                                                               cTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTga                                      <  1:892706/149‑1 (MQ=255)
                                                                                                               cTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTga                                      <  1:554786/149‑1 (MQ=255)
                                                                                                                tGCCAGATCTTGCTGACGCTCTTTCAGCCGGTTGTGATCTCCATCATTGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTgac                                     <  1:52486/149‑1 (MQ=255)
                                                                                                                     gATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTg                                <  1:468843/149‑1 (MQ=255)
                                                                                                                      aTCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGCGATCCTGCGCTATCATCAACTCACACTGACGACTgg                               >  1:24581/1‑149 (MQ=255)
                                                                                                                         ttGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTa                            <  1:832349/149‑1 (MQ=255)
                                                                                                                                 gCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAAATCGCACTGACGATTGGTTATGACGGTa                    >  1:556586/1‑149 (MQ=255)
                                                                                                                                     tttCAGCCGGGGTTCACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGATTATGACGGTCGCTa                >  1:778672/1‑149 (MQ=255)
                                                                                                                                           ccGGGGTTGACCCCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGt          <  2:742648/149‑1 (MQ=255)
                                                                                                                                               ggTTGACCGCCACCATGGCAACCTTCAGTGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGaa      <  2:308343/149‑1 (MQ=255)
                                                                                                                                                gTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGAAg     <  2:779408/149‑1 (MQ=255)
                                                                                                                                                gTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCAAGTCTTGATATCAGCGAAACAGGATGGGAGATCCTGCGCTATCATCAACTCAAACTGACGACTGGTTATGACGGTAGCTACAGGGTTGacg     >  2:838143/1‑147 (MQ=255)
                                                                                                                                                gTTGACCGCCAACAGGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGAAg     <  2:737159/149‑1 (MQ=255)
                                                                                                                                                 ttGACAGCCACCATGACAAAATTCAGCGAAATTAAAAGTATGAATGGATG‑TATAGCCAGGGTGAATAGCAGCGGGACAGGACGGGTGATCCTGCGCTATCATCAACTCACACTGACAACTGGTTATAACGGTAGCTAACGGGTTGaacg    >  2:81266/1‑147 (MQ=255)
                                                                                                                                                   gACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGAAGGTa  <  2:127380/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCCATTTCCTTGAGCAAACTTTAGCTATTCTTATCAATTATGCTTATGGGAGATCTACAGATGATGCCTACGCTTGCTCCACCATCTGTCCTTTCGGCTCCCCAGCGCCGCTGCCAGATCTTGCTGACGCTCTTTCAGCCGGGGTTGACCGCCACCATGGCAACCTTCAGCGAGCTTAATGGTGTGGAT‑GATGATATTGCCAGTCTTGATATCAGCGAAACAGGACGGGAGATCCTGCGCTATCATCAACTCACACTGACGACTGGTTATGACGGTAGCTACCGGGTTGAAGGTA  >  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: ‑