Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi3,214,6600AG20.8% 371.0 / 22.7 173E538E (GAA→GAG) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base A (69/68);  new base G (0/36);  total (69/104)
Fisher's exact test for biased strand distribution p-value = 7.60e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 175 total aligned reads displayed.

CATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTA  >  NC_000913_3_bme_pgi/3214564‑3214758
                                                                                                |                                                                                                  
cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGTTGCCGATGGATTCGGCGACCACCGagag                                                                                                  <  2:1062276/99‑1 (MQ=255)
cATGGAAACGCCGATCGGCGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGc                                                                                                 >  1:321877/1‑100 (MQ=255)
   ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtg                                                                                              >  1:1102060/1‑100 (MQ=255)
    gAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgc                                                                                             >  1:1151012/1‑100 (MQ=255)
    gAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgc                                                                                             >  1:2356482/1‑100 (MQ=255)
    gAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtg                                                                                              <  1:843493/99‑1 (MQ=255)
    gAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtg                                                                                              <  1:3623208/99‑1 (MQ=255)
         gCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGTCGAAGGATTCGGCGACCACCGAGAGCCtgcgtgcg                                                                                        <  2:2096184/100‑1 (MQ=255)
           cGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGc                                                                                      <  1:1034265/100‑1 (MQ=255)
           cGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAACTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGc                                                                                      <  1:3179790/100‑1 (MQ=255)
              tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAAc                                                                                   <  1:236163/100‑1 (MQ=255)
                ggTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACg                                                                                  <  2:3498666/99‑1 (MQ=255)
                ggTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACg                                                                                  >  2:1938355/1‑99 (MQ=255)
                 gTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCa                                                                                >  2:2109546/1‑100 (MQ=255)
                    atgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCacg                                                                              <  1:547255/99‑1 (MQ=255)
                    atgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCacg                                                                              <  1:1194154/99‑1 (MQ=255)
                     tgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGc                                                                                                 >  2:1203449/1‑79 (MQ=255)
                     tgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCacgac                                                                            <  1:2442891/100‑1 (MQ=255)
                       atgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTTGTTTTTTTTAATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCacgacg                                                                           <  2:1417797/99‑1 (MQ=255)
                       atgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGTGCTGCCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCacgacg                                                                           <  2:1973032/99‑1 (MQ=255)
                          aaGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgc                                                                        <  1:3649525/99‑1 (MQ=255)
                           aGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGTGTTGTCGCATGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTgctg                                                                      <  2:252702/100‑1 (MQ=255)
                             aTTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGTGTTGTCGATGGATTCGGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggc                                                                    <  2:2617194/100‑1 (MQ=255)
                             aTTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGTGCTGCCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTgctggc                                                                    <  2:509154/100‑1 (MQ=255)
                                cGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggctg                                                                  <  1:4376553/99‑1 (MQ=255)
                                cGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGctggctgg                                                                 >  1:4639454/1‑100 (MQ=255)
                                 gCATCTGGGGGATTTCATCGAGGATACCACTCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCtggctggc                                                                <  1:4404288/100‑1 (MQ=255)
                                 gCATCTGGGGGATTTCATCGAGGATACCACGCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCtggctggc                                                                <  1:4405581/100‑1 (MQ=255)
                                 gCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGctggctgg                                                                 <  1:2174185/99‑1 (MQ=255)
                                 gCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGctggctgg                                                                 <  1:1313568/99‑1 (MQ=255)
                                  cATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGcc                                                               <  1:523422/100‑1 (MQ=255)
                                  cATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGcc                                                               <  1:2294257/100‑1 (MQ=255)
                                      tGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGAc                                                           >  1:269977/1‑100 (MQ=255)
                                       gggggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGAcc                                                          >  2:2590918/1‑100 (MQ=255)
                                        ggggATTTCATCGAGGATTCCACCCTCTATCTTCTGATGGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCacgacg                                                                           <  2:1535395/82‑1 (MQ=255)
                                        ggggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGAc                                                           <  1:634166/98‑1 (MQ=255)
                                          ggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcg                                                       <  1:950634/100‑1 (MQ=255)
                                            aTTTCATCGAGGATACCACCCTCGTGTTGCCGCAGGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcgcg                                                     <  2:307766/100‑1 (MQ=255)
                                            aTTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcgcg                                                     <  1:3258719/100‑1 (MQ=255)
                                            aTTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcgcg                                                     >  1:960855/1‑100 (MQ=255)
                                            aTTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcgcg                                                     >  1:1857955/1‑100 (MQ=255)
                                            aTTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgcgcg                                                     <  1:2235701/100‑1 (MQ=255)
                                            aTTTCATCGAGGAAACCACCCTCGAGCTGCCGCTTGATTCTGCGACCACCGAAAGCCTGCGTGCGCCAACGCac                                                                               >  2:1376440/1‑74 (MQ=255)
                                             tttCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGt                                                    <  2:767020/100‑1 (MQ=255)
                                               tCATCGAGGATACCACCCTCGTGTTGCCGCAGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGt                                                    <  2:274572/98‑1 (MQ=255)
                                                 aTCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAg                                                >  1:1228327/1‑100 (MQ=255)
                                                 aTCGAGGATACCACCCTCGAGCTGCCGCTGGATTCGGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAg                                                <  2:1223171/100‑1 (MQ=255)
                                                  tCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAg                                                >  1:1750189/1‑99 (MQ=255)
                                                  tCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGGGCCGGCTGGCCTGACCGCGCGTGAAg                                                >  2:4411830/1‑99 (MQ=255)
                                                    gAGGATACCACCCTCGGGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCa                                              <  1:564927/99‑1 (MQ=255)
                                                     aGGATACCACTCTCGAGCTGTCGAAGGATGCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCaaa                                            <  2:4730400/100‑1 (MQ=255)
                                                      ggATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCaaaa                                           <  1:3325992/100‑1 (MQ=255)
                                                           ccaccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCacgacg                                                                           >  2:1778492/1‑63 (MQ=255)
                                                           ccaccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTc                                       <  1:222932/99‑1 (MQ=255)
                                                           ccaccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCt                                      >  1:1686903/1‑100 (MQ=255)
                                                            caccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTg                                     <  2:411436/100‑1 (MQ=255)
                                                            caccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTg                                     >  1:4216490/1‑100 (MQ=255)
                                                            caccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTg                                     <  1:1460649/100‑1 (MQ=255)
                                                             accCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGc                                    <  1:2955694/100‑1 (MQ=255)
                                                             accCTCGAGCTGCCGAATGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGc                                    <  2:584480/100‑1 (MQ=255)
                                                              cccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCg                                   <  2:3882588/100‑1 (MQ=255)
                                                              cccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCg                                   >  1:1517610/1‑100 (MQ=255)
                                                               ccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCg                                   >  1:4741142/1‑99 (MQ=255)
                                                               ccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGt                                  <  1:3789834/100‑1 (MQ=255)
                                                                   gAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTAtn                               <  1:2898583/99‑2 (MQ=255)
                                                                   gAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATg                               <  1:3100493/99‑1 (MQ=255)
                                                                     gCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGCGCGGCGACGCACGACGCGCTGGCTGGCCTGACCGCGCGTGAAGCAACAGTTCTGCGTATGCGt                            >  2:2270005/1‑100 (MQ=255)
                                                                       tGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGttt                          >  1:444373/1‑100 (MQ=255)
                                                                       tGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGttt                          >  1:1257153/1‑100 (MQ=255)
                                                                         ccGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGTTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTc                         <  1:1642246/99‑1 (MQ=255)
                                                                          cGCTGGATTCTGCGACCACCGAAAGCCTGCGTGTGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCgg                       <  1:2091255/100‑1 (MQ=255)
                                                                          cGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCgg                       >  1:3527177/1‑100 (MQ=255)
                                                                            cTGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTc                         <  2:552914/96‑1 (MQ=255)
                                                                            cTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTa                     <  1:4301975/100‑1 (MQ=255)
                                                                              ggATTCTGCGACCACCGAAAGCCTGCGTGCGGCGACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCg                                   >  2:3995137/1‑84 (MQ=255)
                                                                               gATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTc                         <  2:2859038/93‑1 (MQ=255)
                                                                               gATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATc                   <  1:3202235/99‑1 (MQ=255)
                                                                               gATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCg                  <  1:2312396/100‑1 (MQ=255)
                                                                               gATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTAt                    <  2:3814767/98‑1 (MQ=255)
                                                                                aTTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCg                                   <  2:656717/82‑1 (MQ=255)
                                                                                aTTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGa                 >  1:253632/1‑100 (MQ=255)
                                                                                aTTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGa                 <  1:4164303/100‑1 (MQ=255)
                                                                                   cTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGata               <  2:161663/99‑1 (MQ=255)
                                                                                   cTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGatat              >  2:4164327/1‑100 (MQ=255)
                                                                                   cTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGata               >  1:3588733/1‑99 (MQ=255)
                                                                                      cGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATg             <  2:2351756/98‑1 (MQ=255)
                                                                                      cGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGaa           <  1:2598889/100‑1 (MQ=255)
                                                                                        accaccGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAaca         <  2:2686979/100‑1 (MQ=255)
                                                                                         ccaccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAacac        >  1:1538794/1‑100 (MQ=255)
                                                                                         ccaccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAacac        >  1:2594639/1‑100 (MQ=255)
                                                                                          caccGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACAcc       <  2:2672322/100‑1 (MQ=255)
                                                                                          caccGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACAcc       <  2:1821828/100‑1 (MQ=255)
                                                                                           accGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACAcc       >  2:3279712/1‑99 (MQ=255)
                                                                                           accGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACAcc       >  2:3128010/1‑99 (MQ=255)
                                                                                           accGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACAcc       >  1:4045236/1‑99 (MQ=255)
                                                                                           accGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACAcc       >  2:1456489/1‑99 (MQ=255)
                                                                                           accGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCg      <  1:3392256/100‑1 (MQ=255)
                                                                                            ccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGa     >  2:2556363/1‑100 (MQ=255)
                                                                                            ccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGa     >  1:2962545/1‑100 (MQ=255)
                                                                                               aaaGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTa  <  1:3876267/100‑1 (MQ=255)
                                                                                                |                                                                                                  
CATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTA  >  NC_000913_3_bme_pgi/3214564‑3214758

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: ‑