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

ATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTAT  >  NC_000913_3_bme_pgi/3214511‑3214806
                                                                                                                                                     |                                                                                                                                                  
aTGCCTGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCAGAGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGATGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGaa                                                                                                                                                    >  4:25212/1‑150 (MQ=255)
aTGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGaa                                                                                                                                                    >  3:299351/1‑150 (MQ=255)
aTGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGaa                                                                                                                                                    >  2:487400/1‑150 (MQ=255)
aTGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGaa                                                                                                                                                    >  1:470518/1‑150 (MQ=255)
aTGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTAGAGCTGCCGCTGGATTCTGCGACCACCGaa                                                                                                                                                    >  1:388371/1‑150 (MQ=255)
aTGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAACCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGaa                                                                                                                                                    <  3:245212/150‑1 (MQ=255)
           cAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtgc                                                                                                                                         >  3:259835/1‑150 (MQ=255)
            aaGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGTTGCTGGATTCGGCGACCACCGAGAGCCtgcgtgcg                                                                                                                                        <  4:262050/150‑1 (MQ=255)
             aGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCgg                                                                                                                                       >  3:193353/1‑150 (MQ=255)
                  cGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCCGCAACg                                                                                                                                  >  1:320214/1‑150 (MQ=255)
                       aGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGTGTTGTCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCacga                                                                                                                             <  4:261544/150‑1 (MQ=255)
                             gAAGATCGCCAAAGGGCCAATCGCCATGGAAACGCCGATCGGTGATGATGACGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTTGTGCTTTTGAAGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTgct                                                                                                                       <  4:116489/150‑1 (MQ=255)
                             gAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgct                                                                                                                       >  4:330847/1‑150 (MQ=255)
                             gAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgct                                                                                                                       >  3:253994/1‑150 (MQ=255)
                                 aTCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAATCCTGCGTGCGGCAACGCACGACGTgctggct                                                                                                                   >  1:587524/1‑150 (MQ=255)
                                 aTCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggct                                                                                                                   >  1:587533/1‑150 (MQ=255)
                                   cGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTTGTTTTTTTGAAGGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGctggctgg                                                                                                                 <  4:58478/150‑1 (MQ=255)
                                      cAAAGAGCCAATCTCCATGGCAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATAGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCCGCGTGCGGCAAAGCACCACGTGCTGGCTGGccc                                                                                                              >  4:172909/1‑149 (MQ=255)
                                       aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATGTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTg                                                                                                             >  4:308455/1‑150 (MQ=255)
                                        aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGa                                                                                                            <  4:290876/150‑1 (MQ=255)
                                            gCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcgc                                                                                                        >  1:179102/1‑150 (MQ=255)
                                                aTCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGt                                                                                                    <  4:193353/150‑1 (MQ=255)
                                                   tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGTGCTGTCGCATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGaa                                                                                                 <  2:76301/150‑1 (MQ=255)
                                                   tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGaa                                                                                                 <  2:559620/150‑1 (MQ=255)
                                                   tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGaa                                                                                                 <  2:467036/150‑1 (MQ=255)
                                                   tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGaa                                                                                                 <  1:339740/150‑1 (MQ=255)
                                                    ccATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGTGTTTTTGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAg                                                                                                <  2:159864/150‑1 (MQ=255)
                                                    ccATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAg                                                                                                <  3:85558/150‑1 (MQ=255)
                                                      aTGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCa                                                                                              >  1:318311/1‑150 (MQ=255)
                                                      aTGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCa                                                                                              >  2:542638/1‑150 (MQ=255)
                                                              gCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGTCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCt                                                                                      <  4:339573/150‑1 (MQ=255)
                                                              gCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCTGCAAGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCt                                                                                      <  4:24959/150‑1 (MQ=255)
                                                              gCCGATCGGTGATGATGAAGACTCGCATCTGGGGGATTTCATCGAGGATACCACCCACGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCCGCGAGCGCCCAAGCACGACGTGCAGCCTAGCCCGACCGCGCCTGAAGAAGAAGTTcg                                                                                      >  2:8184/1‑149 (MQ=255)
                                                               ccGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTg                                                                                     >  2:284449/1‑150 (MQ=255)
                                                                cGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCAGCCTCGAGCTGTCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGc                                                                                    <  4:176320/150‑1 (MQ=255)
                                                                         atgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGGTCTGCGTATGCGtt                                                                           >  2:333650/1‑150 (MQ=255)
                                                                          tgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTAACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGttt                                                                          <  1:291204/150‑1 (MQ=255)
                                                                           gatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGCGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTc                                                                         >  4:27613/1‑150 (MQ=255)
                                                                            atcaaGATCCGCATCTGGGCGATTTCACCGAGGATACCACCCCCGAGCTGCCGCTGGATGCTGCGACCACCGAAAGCCTGCGTGCGGCAACTCACGACGTGCTGGAGTGCCTGACCGCGCGTGAAGCAAAAGCTCTGCGTATGCGTTTCg                                                                        <  1:346197/147‑1 (MQ=255)
                                                                               aaGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTa                                                                     >  4:53376/1‑150 (MQ=255)
                                                                               aaGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTa                                                                     >  3:313442/1‑150 (MQ=255)
                                                                                   ttCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGa                                                                 <  3:69461/150‑1 (MQ=255)
                                                                                     cGCCTCTGGGGGATTTCAGCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCNTGAAGCAAAAGTNCTGCGTATGCGTTTCGGTATCGata                                                               <  1:333650/150‑1 (MQ=255)
                                                                                       cATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATg                                                             >  2:319373/1‑150 (MQ=255)
                                                                                           tGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAaca                                                         <  4:313442/150‑1 (MQ=255)
                                                                                           tGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAaca                                                         >  1:103217/1‑150 (MQ=255)
                                                                                                  tttCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCTGCAACGCACGACGTGCACGCAGGCCTGACCGCAGGAGAAGACAAAGTTCTGCGAATACGTTTAGGCACCGATACGACCACAGAcaa                                                  >  2:367218/1‑148 (MQ=255)
                                                                                                    tCATCGAGGATACCACCCTCGAGCTGCCGATGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTaca                                                <  4:162755/150‑1 (MQ=255)
                                                                                                       tCGAGGATACCACCCTCGAGCTGCCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGc                                             <  2:269027/150‑1 (MQ=255)
                                                                                                        cGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCt                                            <  2:1508/150‑1 (MQ=255)
                                                                                                         gAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTg                                           <  1:319373/150‑1 (MQ=255)
                                                                                                            gATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGTGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGgaa                                        <  3:86840/150‑1 (MQ=255)
                                                                                                                ccaccCTCGAGCTGCCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGaagaag                                    <  2:238828/150‑1 (MQ=255)
                                                                                                                 caccCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGt                                   <  2:550834/150‑1 (MQ=255)
                                                                                                                  accCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTg                                  <  3:231905/150‑1 (MQ=255)
                                                                                                                   cccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTgg                                 <  3:174588/150‑1 (MQ=255)
                                                                                                                    ccTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGTCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACCCCGACTACACGCTGGCAGAAGTggg                                >  2:579095/1‑150 (MQ=255)
                                                                                                                        gAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTaaa                            <  3:330078/150‑1 (MQ=255)
                                                                                                                         aGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAAc                           <  2:320214/150‑1 (MQ=255)
                                                                                                                           cTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAg                         >  1:63787/1‑150 (MQ=255)
                                                                                                                                   ggATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGt                 >  3:355681/1‑150 (MQ=255)
                                                                                                                                      ttCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTAc              <  3:53376/150‑1 (MQ=255)
                                                                                                                                        cTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTAccc            >  1:126301/1‑150 (MQ=255)
                                                                                                                                        cTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTAccc            >  3:49850/1‑150 (MQ=255)
                                                                                                                                         tGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCtggctggc                                                                                                                >  2:361273/1‑49 (MQ=255)
                                                                                                                                         tGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCtggctggc                                                                                                                <  1:361273/49‑1 (MQ=255)
                                                                                                                                           cGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCcgcg         <  3:330847/150‑1 (MQ=255)
                                                                                                                                            gACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGa        <  3:327737/150‑1 (MQ=255)
                                                                                                                                            gACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGa        >  3:136745/1‑150 (MQ=255)
                                                                                                                                               caccGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACg     >  4:25830/1‑150 (MQ=255)
                                                                                                                                                accGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGt    >  1:409321/1‑150 (MQ=255)
                                                                                                                                                  cGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTAt  <  2:54596/150‑1 (MQ=255)
                                                                                                                                                  cGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACGCGGTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTAt  >  4:360731/1‑150 (MQ=255)
                                                                                                                                                  cGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACGCGATGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTAt  >  4:360827/1‑150 (MQ=255)
                                                                                                                                                    aaaGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTc                      >  2:306167/1‑128 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTTACCCGCGAACGTAT  >  NC_000913_3_bme_pgi/3214511‑3214806

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