Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,8000AG25.3% 153.7 / 55.3 92E1180E (GAA→GAG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (26/42);  new base G (13/10);  total (39/53)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTAC  >  NC_000913_3_bme_pgi/526652‑526949
                                                                                                                                                    |                                                                                                                                                     
cgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAat                                                                                                                                                      >  1:694154/1‑150 (MQ=255)
  cAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                                    >  1:663325/1‑150 (MQ=255)
   aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatg                                                                                                                                                   <  1:640399/150‑1 (MQ=255)
     ggCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                                 <  1:731131/149‑1 (MQ=255)
      gCGGGCTGACGGTGGAGCAACTGGCCCGCCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTCTCCCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCATCGAAGACGGAAATACGGCGTGGCGCGGGGAGTatgatg                                                                                                                                                <  1:522451/150‑1 (MQ=11)
      gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgatg                                                                                                                                                <  2:926841/150‑1 (MQ=255)
         ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAAt                                                                                                                                             >  2:16379/1‑150 (MQ=255)
            tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                                          <  1:752542/150‑1 (MQ=255)
               cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCa                                                                                                                                       <  2:7644/150‑1 (MQ=255)
               cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCa                                                                                                                                       <  2:763205/150‑1 (MQ=255)
               cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatg                                                                                                                                                   <  2:265073/138‑1 (MQ=255)
               cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatg                                                                                                                                                   >  1:265073/1‑138 (MQ=255)
               cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCa                                                                                                                                       <  2:138424/150‑1 (MQ=255)
               cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCa                                                                                                                                       >  1:122851/1‑150 (MQ=255)
               cGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCa                                                                                                                                       <  2:832671/150‑1 (MQ=255)
                 gTGGAGCAACTTGCCAGACAGGGGGAGCCGGAATACACACCGGCGCGAAAAGCGCATCTTTATCACGGCGACCCCCGGTGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                                     <  1:852672/150‑1 (MQ=17)
                 gTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                                     >  2:783046/1‑150 (MQ=255)
                 gTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                                     >  2:83068/1‑150 (MQ=255)
                  tGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAcc                                                                                                                                    <  2:593678/150‑1 (MQ=255)
                   ggAGCAACTGGCCAGACAGGTTGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCa                                                                                                                                   <  1:240976/150‑1 (MQ=255)
                    gAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAg                                                                                                                                  >  2:795466/1‑150 (MQ=255)
                     aGCAACTGGCCAGCCAGGTGGAGCCGGAATACACACCGGCGCGCAAAGCTCATCTTGATCACTGCGACAACCGGGGACTGCCGCTGGCGCTTATCAGCGACGACGGCAATACGGCGTGGCGCGGGGAGGATGATGATTGGGGCAACCAGc                                                                                                                                 >  2:768356/1‑150 (MQ=11)
                     aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGc                                                                                                                                 >  1:748944/1‑150 (MQ=255)
                     aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAcc                                                                                                                                    <  1:276959/147‑1 (MQ=255)
                     aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAcc                                                                                                                                    >  2:276959/1‑147 (MQ=255)
                     aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGc                                                                                                                                 <  1:277521/150‑1 (MQ=255)
                     aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAACCACATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCACCGAACACGGCAATACGGCGTGGCGAGGGGAGTATGATGAATGCGGCACCCatc                                                                                                                                 >  2:437668/1‑148 (MQ=11)
                     aGCAAAAGGCCAGACAGGTGGAGACGGAATAACCACCGGCGCGAAAAGCTCATCTTTATCACTGCCACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAAAATGATGAAAGGGGCAACCAGc                                                                                                                                 <  1:842216/150‑1 (MQ=17)
                      gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCt                                                                                                                                >  2:353600/1‑150 (MQ=255)
                      gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCt                                                                                                                                <  1:581011/150‑1 (MQ=255)
                         aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGCCCCCCGGGGACTGCCGCTGGCGCTTATAAGCGAACACGGCAATACGGCGTGGAGCGCGGAAGATGATGTATGGGGCAACCAGCTTaa                                                                                                                             >  1:322885/1‑150 (MQ=17)
                         aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCTATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTaa                                                                                                                             >  2:37693/1‑150 (MQ=255)
                           tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAAtt                                                                                                                           >  2:443089/1‑149 (MQ=255)
                           tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATg                                                                                                                           <  1:708418/150‑1 (MQ=255)
                              ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGCCCACCGGGGACTGCCGCTGGAGCTTATAAGCGAAGACGGCAAAACGGCGTGGAGCGCGGAATAAGATGAATGGGGCAACCAGCATAATgagg                                                                                                                        >  1:73130/1‑150 (MQ=17)
                              ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                        >  2:567079/1‑150 (MQ=255)
                              ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                        <  2:205135/150‑1 (MQ=255)
                               cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                       >  1:329154/1‑150 (MQ=255)
                               cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTAATATTTATCACTGCGCCCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGGAGCCAGCTTAaggaggg                                                                                                                       >  2:407435/1‑149 (MQ=11)
                                aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCGTAATgaggag                                                                                                                      >  1:875440/1‑150 (MQ=255)
                                 gACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGgaga                                                                                                                     <  1:290753/150‑1 (MQ=255)
                                 gACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGgaga                                                                                                                     <  2:663325/150‑1 (MQ=255)
                                  aCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGaa                                                                                                                    <  1:775141/150‑1 (MQ=255)
                                  aCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGaa                                                                                                                    <  1:588230/150‑1 (MQ=255)
                                     ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                                 >  2:582672/1‑150 (MQ=255)
                                      gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                                <  2:154505/150‑1 (MQ=255)
                                       tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGTATGGGGCAACCAGCTTAATGAGGAGAACCGGc                                                                                                               >  2:388779/1‑150 (MQ=11)
                                       tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGCCAATACGGCGGGGCGCGGGGAGAATGATGAATGGGGCAACCAGCTTAATGAGGAGACCCCGc                                                                                                               >  2:617537/1‑150 (MQ=11)
                                       gggAGCCGGAATACACCCCGGCGCGAAAAGCTCATCTTTATCCCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGc                                                                                                               <  1:513782/149‑1 (MQ=11)
                                          aGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGAGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCACCCAGCTTAATGAGGAGAACCCGcatc                                                                                                            >  1:278655/1‑150 (MQ=255)
                                           gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatca                                                                                                           >  1:493150/1‑150 (MQ=255)
                                           gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatca                                                                                                           >  1:268290/1‑150 (MQ=255)
                                             cGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGGGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtg                                                                                                         <  1:647696/150‑1 (MQ=255)
                                               gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgtg                                                                                                       <  1:228474/150‑1 (MQ=255)
                                                aaTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGGATGATGAATGGGGCACCCAGCGTAAGGACGAGAACCCGcatcacctgc                                                                                                      >  2:739957/1‑145 (MQ=11)
                                                 aTACAGACCGGCGCGAAATGATCATATTTATCAATGCGACCAGCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCTATACGGCGTGGAGCGCGGAATATGATGAATGGGGCATCCAGCTTAATGAGGAGAACCCGCATCATGTGTa                                                                                                     <  1:510334/150‑1 (MQ=17)
                                                 aTACACACCGGCGCGAAAAGCGCCTCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCATATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTa                                                                                                     <  1:97631/150‑1 (MQ=255)
                                                  tACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGCCCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctgcac                                                                                                    <  1:867845/150‑8 (MQ=11)
                                                     acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAg                                                                                                 <  2:302053/150‑1 (MQ=255)
                                                        cccgcgcGAAACGCTCATCTATATCACTGCGACCACCGGGGACTGCCGCTCGCGCTTATCAGCGAAGTCGGCAATACGGCGGGGAGCGCGGATTATGATGAATGGGGCATCCATCTTTATGAGGAGAACCGGCATCATGTGTATCAGCCg                                                                                              <  1:608485/147‑1 (MQ=255)
                                                          ggCGCGAGAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTa                                                                                            >  1:714126/1‑150 (MQ=255)
                                                           gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTAc                                                                                           <  2:694154/150‑1 (MQ=255)
                                                                 aaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTg                                                                                     <  1:795466/150‑1 (MQ=255)
                                                                   aGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGcc                                                                                   >  2:314360/1‑150 (MQ=255)
                                                                     cTCATCTTTATCACTGCGACCACCGGTGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCACATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                                 <  2:906417/150‑1 (MQ=255)
                                                                     cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                                 >  1:330900/1‑150 (MQ=255)
                                                                     cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                                 <  2:116015/150‑1 (MQ=255)
                                                                     cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCACATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                                 <  2:455881/150‑1 (MQ=255)
                                                                     cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCACATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                                 <  2:468396/150‑1 (MQ=255)
                                                                     cTCATCTTTATCACTGCGACCACCGGCGACTGACGCTGCCTCTTATCTGCGAAGAAGGCAAGACGGCGTGGATCGCGGAATATGATGAATGGGGCAAGCAGCGTCATGAGCAGATCCCGAATGATATGTATCAGACGAACAGACTCCCAg                                                                                 >  2:908402/1‑150 (MQ=255)
                                                                     cTCACCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCACATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                                 <  2:776638/150‑1 (MQ=255)
                                                                        aTCTTTATCACTGCGACCACCTGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTACCAGGgc                                                                              <  2:503728/150‑1 (MQ=255)
                                                                        aTCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgc                                                                              <  2:329154/150‑1 (MQ=255)
                                                                          cTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                                                                            <  1:83068/150‑1 (MQ=255)
                                                                                 aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGGGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatg                                                                     >  2:51185/1‑150 (MQ=255)
                                                                                     cGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                                 >  1:463296/1‑150 (MQ=255)
                                                                                         caccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGaa                                                             >  1:433749/1‑150 (MQ=255)
                                                                                          accGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGGGTATCCGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAAt                                                            >  2:623736/1‑150 (MQ=255)
                                                                                                gACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGc                                                      <  1:16379/150‑1 (MQ=255)
                                                                                                  cTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTg                                                    <  2:122851/150‑1 (MQ=255)
                                                                                                  cTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTg                                                    <  1:582672/150‑1 (MQ=255)
                                                                                                     ccGCTGGAGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTAc                                                 <  2:278655/150‑1 (MQ=255)
                                                                                                                  tCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCAcc                                    <  2:268290/150‑1 (MQ=255)
                                                                                                                   cAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCg                                   >  1:877623/1‑150 (MQ=255)
                                                                                                                    aGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCgg                                  >  1:70711/1‑150 (MQ=255)
                                                                                                                     gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGTATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGt                                 <  1:58275/150‑1 (MQ=255)
                                                                                                                     gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGt                                 <  1:26352/150‑1 (MQ=255)
                                                                                                                                  tACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGtt                    >  1:599407/1‑150 (MQ=255)
                                                                                                                                  tACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGCCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGtt                    >  1:148017/1‑150 (MQ=255)
                                                                                                                                      gCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAg                <  2:565954/150‑1 (MQ=255)
                                                                                                                                          ggAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGc            <  1:927756/150‑1 (MQ=255)
                                                                                                                                            atcgcgGAATATGCTCAATGCGGCAACCAGCTTAATGAGGAGAACCCGCATGATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGCGGAAGCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCgg          <  2:704049/148‑1 (MQ=17)
                                                                                                                                                  gAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATAtt    <  2:714126/150‑1 (MQ=255)
                                                                                                                                                  gAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCCTCACCGGTAATACGATCCGTTGCAGGGGCGGTATAtt    >  2:481346/1‑150 (MQ=255)
                                                                                                                                                    atatGATGAATGGGGCAACCAGCTTACTGAGGAGAACCCGCATCATGTGTATCAGCCGTATCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTAc  <  1:593453/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTAC  >  NC_000913_3_bme_pgi/526652‑526949

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