Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7960CG32.0% 110.2 / 66.3 76A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (21/30);  new base G (8/16);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGT  >  NC_000913_3_bme_pgi/526648‑526942
                                                                                                                                                    |                                                                                                                                                  
cTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGGCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgcgc                                                                                                                                                    <  2:743900/149‑1 (MQ=255)
    cgcAGTGGGGGCGGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGaa                                                                                                                                                <  2:310454/149‑1 (MQ=255)
     gcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACAACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAat                                                                                                                                               <  1:238939/149‑1 (MQ=255)
      caggtcGGGCTGACGGTGGAGCAACTGGCTAGACAGGTGGAGCCGGAATACACACCGGCGCTAAATGCTCATCTTTATCACTTCTACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACTGCGTGGAGCGCGGAata                                                                                                                                              <  2:366498/144‑1 (MQ=17)
       aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAgtat                                                                                                                                             <  1:680296/149‑1 (MQ=255)
        ggccGGGCCGACGGTGAAGCGACTGGCCAGACAGGTGGAGCCGAAATACACATCGGCGCGAAATGCTCATCTTAAGCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatg                                                                                                                                            <  2:529682/146‑1 (MQ=14)
         tGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACCAGCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                           <  1:194009/149‑1 (MQ=255)
             ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGaa                                                                                                                                       <  1:548249/149‑1 (MQ=255)
             ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGCACTGCCCCTGGCGCTTATCCGCGAAGACGGCCATACGCCCTGGGCCGGGGAAGATgaagaa                                                                                                                                       >  2:362055/1‑149 (MQ=14)
              gCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                                                                      <  2:530322/149‑1 (MQ=255)
                tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                                                    >  2:786990/1‑149 (MQ=255)
                tGACGGTGGAGCAACGGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACAACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                    <  2:621154/149‑1 (MQ=255)
                  aCGGTGGAGCAACTGGCCAGACAGGTAGAGCCGGAATACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                                  <  1:532336/149‑1 (MQ=11)
                    ggTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCgggg                                                                                                                                                  <  2:642182/131‑1 (MQ=255)
                    ggTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCgggg                                                                                                                                                  >  1:642182/1‑131 (MQ=255)
                       ggAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAccc                                                                                                                             <  2:454089/149‑1 (MQ=255)
                       ggAGCAACTGGCCAAACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCCGGCGCTTATCAGCGAAGACGGCAATACGACGTGGCGCGGGGAGTATGATGAATGGGGCAAcc                                                                                                                             >  2:520504/1‑149 (MQ=11)
                              cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGGTTaa                                                                                                                      >  1:473981/1‑149 (MQ=255)
                               tgcccATACAGGTGGAGCCGGAATACACTCCGGCGCGAAAAGCTCATCTTTATTACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAAt                                                                                                                     <  1:291428/146‑1 (MQ=255)
                               tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCCCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAAt                                                                                                                     <  2:434545/149‑1 (MQ=11)
                               tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAAt                                                                                                                     <  2:552799/149‑1 (MQ=255)
                               tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAAt                                                                                                                     <  1:473845/149‑1 (MQ=255)
                                 gtcacaCAGGGGGGGCGGGAATACACACCGGTGCGATTATCTCATCTTTATCACTGCGTCCACCGGGGACTGCCGCTGGCGCTTATCATCGAAGACGGTAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                   <  2:647701/144‑1 (MQ=17)
                                  ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                  <  2:600357/149‑1 (MQ=255)
                                  ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                  <  1:809534/149‑1 (MQ=255)
                                  ccAGACAGGTGGAGCCGGAATACACAACGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                  <  2:463339/149‑1 (MQ=255)
                                   cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                 <  1:786990/149‑1 (MQ=255)
                                    gaaCAGGAGGAGCCGGAATACACACCGGCGCCATAAGCTCAGCTCTATCACGGCGACCACCGGGGACTGCCGCCGGCGCTTATCAGCGAAGACAGCAATACGGCGTGGCGCGGGGAGTATGTTGAATGGGGCAACCAGCTTAATgagga                                                                                                                <  2:648698/147‑1 (MQ=11)
                                    aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGTATGGGGCAACCAGCTTAATgagga                                                                                                                <  1:723561/149‑1 (MQ=255)
                                    aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                <  2:114864/149‑1 (MQ=255)
                                      aCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCACTACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGgaga                                                                                                              <  1:553746/149‑1 (MQ=255)
                                       cAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGaa                                                                                                             <  2:436012/149‑1 (MQ=255)
                                        aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAc                                                                                                            <  1:542579/149‑1 (MQ=255)
                                         gggggAGCCGGAAAACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAcc                                                                                                           <  1:725735/146‑1 (MQ=11)
                                          gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCGACCAGATGATTGGGGAGAAccc                                                                                                          >  2:566977/1‑149 (MQ=11)
                                           tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAATACGGCAATACGGCGTGGCGCGGGGAGGATGATGGATGAGGCAACCAGCTTAATGCGGCGAAcgcg                                                                                                         >  2:279172/1‑149 (MQ=11)
                                            ggAGCGGGAATACACACCGGCGCGAAAAGCTCATCTTTATCGCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGc                                                                                                        <  1:561728/149‑1 (MQ=11)
                                              aGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGCATATGATGAATCGGGAAACCAGCTTAATAAGGAGAACCCGcat                                                                                                      >  2:582898/1‑149 (MQ=255)
                                               gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGGGGAGAACCCGcatc                                                                                                     >  1:349386/1‑149 (MQ=255)
                                                ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAATACTGCAATACGTCGTGGGGCGCGGCGTCTGATGAGTGGGGCAACCAGCTTAATGAGGAGAACCCGcatca                                                                                                    >  2:96/1‑149 (MQ=14)
                                                ccGGAATACACACCGGCGCGAAAAGCTCAGCTGTATCACTGCGACCACCGGGGACTGCCGCAGGCGCCTATCAACGAAGACCGCAATACGGCGTGGAGGGCGGAATTGGAAGAATGGGGCAAACAGCTTAAGGAGGAGAACCCGcatca                                                                                                    >  2:678890/1‑149 (MQ=17)
                                                 cGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGCTGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcat                                                                                                   >  1:67723/1‑149 (MQ=255)
                                                   gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                 <  1:550324/133‑1 (MQ=255)
                                                   gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                 >  2:550324/1‑133 (MQ=255)
                                                   gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATAAGCGGTGATGGGCTGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgg                                                                                                 >  1:193826/1‑148 (MQ=17)
                                                     aTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATgtgt                                                                                               >  2:800820/1‑149 (MQ=255)
                                                       gcacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAt                                                                                             <  1:629330/148‑1 (MQ=255)
                                                       acacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGGGGAGTATGATGAATGGGGCAAACAGCTTAATGAGGAGAAACCGAATAAcgtgcaa                                                                                             >  1:256163/1‑146 (MQ=14)
                                                         acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGCACCCGc                                                                                                        >  2:510042/1‑136 (MQ=255)
                                                         acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGc                                                                                                        <  1:510042/136‑1 (MQ=255)
                                                            ccGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGcc                                                                                        >  1:268176/1‑149 (MQ=255)
                                                                   gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                                 >  1:330167/1‑149 (MQ=255)
                                                                       aGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGc                                                                             <  2:388097/149‑1 (MQ=255)
                                                                         cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCa                                                                           <  1:100025/149‑1 (MQ=255)
                                                                         cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCa                                                                           >  1:419491/1‑149 (MQ=255)
                                                                            aTCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                         >  2:289099/1‑116 (MQ=255)
                                                                            aTCTTTATCACTGCCACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                         <  1:289099/116‑1 (MQ=255)
                                                                                ttATCGCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATTAGGAGAAACCGCATCCTGTCTATCAGCCGTACCGTCTGCCAGGgcagc                                                                    >  1:738407/1‑149 (MQ=17)
                                                                                 tATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACGTGTATCAGCCGTACCGTCTGCCACGgcagca                                                                   <  1:381202/149‑1 (MQ=255)
                                                                                  aTCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcag                                                                  >  1:86489/1‑149 (MQ=255)
                                                                                   tCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                           <  1:159874/75‑1 (MQ=255)
                                                                                   tCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                           >  2:159874/1‑75 (MQ=255)
                                                                                    cagtGCGGCCAGGGGGGGCGGCCGGTGGCGGTTATCAGCGAAGACGGCAAGACGGCGTGGAGCGCGGAATAGGATGAATGGGGCAACCAGCTTAATGAGGGGAACCCGAATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcagca                                                                <  2:230423/146‑1 (MQ=17)
                                                                                            ccaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGCTGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGCTGAgg                                                        >  1:447058/1‑149 (MQ=255)
                                                                                               ccGGGGAATGCCGCTGGCGCTGATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGACTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAAt                                                     <  2:754091/149‑1 (MQ=255)
                                                                                                 ggggACTTCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCa                                                   <  2:112674/149‑1 (MQ=255)
                                                                                                      attcggCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCt                                              <  2:419491/144‑1 (MQ=255)
                                                                                                       tGCCGCTGGCGCTTATCAGCGAAGTCGGCAATACGGCGTGGAGCGCGGAATATGATGAAt                                                                                                                                      >  1:725107/1‑60 (MQ=255)
                                                                                                       tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAAt                                                                                                                                      <  2:725107/60‑1 (MQ=255)
                                                                                                          cGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCACATCATGTGTATCACCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTAc                                          <  1:735106/149‑1 (MQ=255)
                                                                                                                 gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGCGGAGCACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATGAGGGCTGTACTATAAcc                                   >  2:242464/1‑149 (MQ=255)
                                                                                                                 gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                   <  1:833262/149‑1 (MQ=255)
                                                                                                                 gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                   <  1:48304/149‑1 (MQ=255)
                                                                                                                         gCGAAGACGGCAATACGGCGTGGAGAGAGTAAGATGATGAATGGGGCAACCAGCTTAATGAGGGGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAAGCCAGCAGCATGATGAGGAACCAGGGCTGTACTATAACCGTTACCgg                           >  2:308106/1‑149 (MQ=17)
                                                                                                                          cGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGt                          >  1:719820/1‑149 (MQ=255)
                                                                                                                                 ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGa                   >  1:393916/1‑149 (MQ=255)
                                                                                                                                  gCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGAt                  <  1:800820/149‑1 (MQ=255)
                                                                                                                                          gCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCa          <  1:582898/149‑1 (MQ=255)
                                                                                                                                               gAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGc     >  2:739860/1‑149 (MQ=255)
                                                                                                                                                  cgcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCggg  >  1:67136/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGT  >  NC_000913_3_bme_pgi/526648‑526942

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