Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7960CG26.9% 179.2 / 69.4 104A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (28/48);  new base G (9/19);  total (37/67)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
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 107 total aligned reads displayed.

TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTA  >  NC_000913_3_bme_pgi/526649‑526943
                                                                                                                                                   |                                                                                                                                                   
ttGCGCAGGGCGGGCTGACGGTGGAGCAACTGGCCAGCCAGGTGGCGCCGGAATACACACCGGCGCGCAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCgg                                                                                                                                                   <  2:745357/150‑1 (MQ=255)
gtGCGCAGGGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCgg                                                                                                                                                   <  1:331296/149‑1 (MQ=255)
    gcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAata                                                                                                                                               >  1:222894/1‑150 (MQ=255)
     cAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGTTGGAGCCCGAATACACACCGGCGCGAAAAGCTAAACTTTATCACTGCGACCACCGGGGACTGCCGCTCGAGCTTATCAGCGAAGACGGCAATACGGCGTCGAGCGCGGAatat                                                                                                                                              >  1:362523/1‑150 (MQ=17)
     cAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGCGGAatat                                                                                                                                              >  1:319276/1‑150 (MQ=255)
     cAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                              >  1:763339/1‑150 (MQ=255)
     cAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                              >  1:769753/1‑150 (MQ=255)
     cAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATAAGCGAAGACGGCAAAACGGCGGGGCGCGGGGaggat                                                                                                                                              >  2:724343/1‑146 (MQ=11)
      cgggCGGGCAGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACCCACCGGCGCGAAAAGCGCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatg                                                                                                                                             <  2:905660/147‑1 (MQ=255)
      cgTGCGGGCTGACGGTGGAGCAACTGGCCAGACCGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatg                                                                                                                                             <  1:258702/149‑1 (MQ=255)
      aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatg                                                                                                                                             <  1:159874/150‑1 (MQ=255)
       gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGTATACACACCGGCGCGAAAAGCTCATCTTTATCACTGAGACCACCGGGGACTGACGCTGGCGCTTATCAGCGACGACCGCACAACGGCGGGGAGAGCGGActatca                                                                                                                                            >  1:698297/1‑148 (MQ=14)
       gTGCGGGCTGACGGCGGAGCAACTGGCCAGACAGGCGGAGCCGGAATACACACCGGCGCGAATATCTCATCTTTAGCACAGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                            <  1:646113/150‑1 (MQ=14)
        tGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGCCCCCCGGGGACTGCCGCTGGCTCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                           >  2:313727/1‑150 (MQ=255)
        tGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGGGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                           <  1:79388/150‑1 (MQ=255)
        tGCGGGCTGACGGTGCAGCAACTGGCCAGACAGCTGGAGGCGGAATACACACCGTCGCGAAAAGCTCATCTTAATCACTGCGACCACCGGGGACTCCCGCTGGGGCTTATCAGCGAAGACGGCTATACTGCGAGGCGCGGGGAGTatgat                                                                                                                                           <  1:663897/150‑1 (MQ=255)
        tGCGGGCGGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACCCCGGCGCGTAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTTTCAGCGAAGAGGGCAATACGGCGTGGCGCGGGGAGTatgat                                                                                                                                           <  1:662954/150‑1 (MQ=11)
        tGCAGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCCCCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                           >  2:218261/1‑150 (MQ=255)
         gCGGGCTTAAGGTGGAGCAACTGGCCAGACTGGAGGAGCCTGAATAGACACCGGCGCGAAAAGATCATCTTTATCACTGCGACCACCGGTGACTGCCGCTGGAGCTTATCAGCTAAGACTGCAATACGGCGTGGAGCGCGGAATatgatg                                                                                                                                          <  1:67956/150‑1 (MQ=14)
         gCGGGCTGACGGTGGCGCAACTGGCCAGACCGGTGGAGCCGGACTACCCACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgatg                                                                                                                                          <  2:934670/150‑1 (MQ=255)
         gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgatg                                                                                                                                          >  1:665177/1‑150 (MQ=255)
         gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgatg                                                                                                                                          >  2:830869/1‑150 (MQ=255)
            ggCTGACGGTGGAGCAACTGGCCAGACCGGTGGAGCCTGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                                                                       <  2:342173/150‑1 (MQ=11)
            ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                                                                       <  2:150166/150‑1 (MQ=255)
            ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                                                                       <  2:755289/150‑1 (MQ=255)
              cTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGACGTGGAGCGCGGAATATGATGAATgg                                                                                                                                     >  1:873526/1‑150 (MQ=255)
               tGACGGTGGATCACCTTTCCAGACAGGTGTAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACTGCAATACGGCGTGGAGCGCGGAATATGATGAATggg                                                                                                                                    <  1:783075/150‑1 (MQ=17)
               ggACGGTGGAGCAACTCGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                                    <  1:443565/149‑1 (MQ=255)
                gACGGTGGAGCAACTGGCCCGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgggg                                                                                                                                   <  2:155359/150‑1 (MQ=255)
                 tcGGTGGAGCAACTGGCCAGACAGGTTGAGCCGGAATACACACCGGCTCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGc                                                                                                                                  <  1:595271/149‑1 (MQ=11)
                  cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCa                                                                                                                                 <  1:819933/150‑1 (MQ=255)
                  cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCa                                                                                                                                 <  1:245840/150‑1 (MQ=255)
                      ggAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCa                                                                                                                             >  2:923251/1‑150 (MQ=255)
                       gAGCAACTGGCCAGACAGGTGGAGACGGAATACACACCGGCGCGAAAAGCTCATATTAATCACTGCGACCACCGGGGACTGCCGCTCGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGCTGACTGGGGCAACCAg                                                                                                                            >  1:656677/1‑150 (MQ=11)
                        aGCAACTGGCCAGACAGGTGGAGCCGGAATACACTCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGGAGGATGAATGGGGCAACCAGc                                                                                                                           >  2:781160/1‑150 (MQ=11)
                         gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGCAGACGGCAATACGGCGTGGCGAGGGGAGTATGATGAATGGGGCAAGCAGCt                                                                                                                          >  1:625075/1‑150 (MQ=11)
                          catcTGGCCAGACAGTTGGAGCCGGAATACACAGCGGCGCGAAAAGCTCAGCTTTATCCCTGCGAACACCTGTGACTGCGGCTGGGGGTTATAAGCGAAGACGGCTATACGGCGTGGCGCGGTTAGTATGATGAATGGGGCAACCAGCtt                                                                                                                         <  1:389495/147‑1 (MQ=11)
                          cAACTGGCCATACTGGTGGAGCCGGAATACACACCGGCGCGTAAAGCTCATCTTTATCCCTGGGACCTCCGGGGTCTGCCGCTGGCGCTTTTCCGCGAAGTCGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCATCCAGCtt                                                                                                                         <  1:546421/150‑1 (MQ=11)
                          cAACTGGCCAGACTGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGGGCGGAATATGATGAATGGGGCAACCAGCtt                                                                                                                         <  2:925679/150‑1 (MQ=255)
                           aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTa                                                                                                                        <  2:43761/150‑1 (MQ=255)
                           aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTa                                                                                                                        >  1:239176/1‑150 (MQ=255)
                              tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATg                                                                                                                     <  2:110699/150‑1 (MQ=255)
                               tgCCAGACAGGTTGAGCCGGAATACACACCGGCGCGTAAAGATCATCTTTATAACTGCGCCCACCGGGGACTGCAGCTGGCGCTTTTCAGCGAAGACGGCAATAAGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                    <  2:625618/149‑1 (MQ=14)
                               ggacataCAGGTGGATCAGGAATACACACCGGCGCGAAAAGCTCATCTTTATCTCTGCGACCTCCGGGGACTGCGGATGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                    <  1:300768/144‑1 (MQ=17)
                               ggCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGTCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATga                                                                                                                    >  1:831474/1‑150 (MQ=255)
                                 ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTTCGACCTCCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                  <  1:721324/150‑1 (MQ=11)
                                  cAGTCAGGTGGAGCCGGAATACACAGCTGCGCGTAACGCTCATCTTTCTCTCTGCGCCCTCCGGGGACTGCCGCTGGCGGTTATCAGCGAAGTGGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                 <  2:77201/150‑1 (MQ=11)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                 <  1:53123/150‑1 (MQ=255)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                 <  1:494382/150‑1 (MQ=255)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                 <  2:317771/150‑1 (MQ=255)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                 <  2:445242/150‑1 (MQ=255)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                 <  1:808825/150‑1 (MQ=255)
                                  cAGACAGGTGCTGGCGGAATACACACCGGCGCGAAAAGCTCATCTATATCACAGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGTTGAATGGGGCAACCTGCTTAATgagga                                                                                                                 <  1:629537/150‑1 (MQ=11)
                                   aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgaggag                                                                                                                <  2:918505/150‑1 (MQ=255)
                                    gACAGGTGGAGCCGGAATACCCACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGgaga                                                                                                               >  2:310372/1‑150 (MQ=255)
                                       aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACAGCGACCACCGGGGACTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAcc                                                                                                            <  1:205489/150‑1 (MQ=255)
                                       aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCAGCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAcc                                                                                                            <  1:208719/150‑1 (MQ=255)
                                           ggAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                        <  2:367478/150‑1 (MQ=255)
                                           ggAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                        <  2:808541/150‑1 (MQ=255)
                                                  gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCTCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctg                                                                                                 <  1:853813/150‑5 (MQ=11)
                                                    aTACACACCGGCGCGAAATGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTa                                                                                               <  1:226417/150‑1 (MQ=255)
                                                     tACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGCCCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAt                                                                                              <  2:258681/150‑1 (MQ=255)
                                                          acCCGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGcc                                                                                         >  1:838496/1‑150 (MQ=255)
                                                              gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTAc                                                                                     <  2:672522/150‑1 (MQ=255)
                                                                   aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCt                                                                                <  1:479546/150‑1 (MQ=255)
                                                                   aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCCTGGAGCGCGGAATATCATGAATGGGGCAACCAGCTAAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCt                                                                                >  2:462304/1‑150 (MQ=255)
                                                                     aaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGc                                                                              <  2:70783/150‑1 (MQ=255)
                                                                       tcgcATCTTTATCACTGCGACCAGCGGGGACTGCCGCTGGCGCTTATAAGCGAAGACGGCTATACGGCGTGGAGCGCGGATTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCa                                                                            <  1:856737/147‑1 (MQ=17)
                                                                        cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGAGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                           <  1:416674/150‑1 (MQ=255)
                                                                                tATCACAGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGAAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcag                                                                   <  2:758560/150‑1 (MQ=255)
                                                                                    aCTGCGCCCACCGGGGCCTGCCGCTGGCGCTGATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGCAAGGGGCAACCAGCTTAATGAGGAGCACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGTGCAGCCGCatg                                                               >  2:354570/1‑150 (MQ=14)
                                                                                     atGCGACCACCGGTGACGTCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatga                                                              <  2:529419/149‑1 (MQ=255)
                                                                                         gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                          <  2:856801/150‑1 (MQ=255)
                                                                                          accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAgg                                                         <  2:239176/150‑1 (MQ=255)
                                                                                          accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAgg                                                         <  1:5985/150‑1 (MQ=255)
                                                                                                      tGCCGCTGGCGCTTATCAGCGAAGACGTCAATACGGCGTGGAGCGCGTAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGt                                             <  2:144715/150‑1 (MQ=255)
                                                                                                      tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGt                                             <  2:148332/150‑1 (MQ=255)
                                                                                                        ccGCTGGCGCTTATCAGAGAAGACGTCAATACGGCGTGGAGCGCGGATTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTAc                                           <  2:873526/150‑1 (MQ=255)
                                                                                                               cgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                    >  2:796174/1‑150 (MQ=255)
                                                                                                                gcTTATGAGCGAAGACGGCAATACTCCGTGGAGCTCGTAATATGATGATTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGGATCAGAATGATGAGTAATCAGGGCTGTTCTATAAcgg                                   <  2:362523/150‑3 (MQ=17)
                                                                                                                gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCg                                   <  2:319276/150‑1 (MQ=255)
                                                                                                                gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCg                                   <  2:769753/150‑1 (MQ=255)
                                                                                                                      cAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAAGTGTATCAGCCGTACCGTCTGCCAGTGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCg                             >  2:246911/1‑150 (MQ=255)
                                                                                                                       aGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCgg                            <  1:830869/150‑1 (MQ=255)
                                                                                                                         cGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGta                          >  1:906248/1‑150 (MQ=255)
                                                                                                                         cGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGta                          >  1:558388/1‑150 (MQ=255)
                                                                                                                         cGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGta                          >  1:918620/1‑150 (MQ=255)
                                                                                                                          gAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtaa                         >  1:747043/1‑149 (MQ=255)
                                                                                                                           aaGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtact                        <  1:218261/150‑1 (MQ=255)
                                                                                                                           aaGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtact                        <  2:816933/150‑1 (MQ=255)
                                                                                                                           aaGACGGAAATACGGCGTGGAGCGCGGAATATGATGAAGGCGGCAACCAGCTTAATGAGGCGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtact                        <  1:313727/150‑1 (MQ=255)
                                                                                                                               cgccAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATTTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATACGGGCTGTACTATAACCGTCACCCGTACTACGa                    >  2:210266/4‑150 (MQ=17)
                                                                                                                               cGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGa                    >  2:274978/1‑150 (MQ=255)
                                                                                                                                 gCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGAta                  >  1:365304/1‑149 (MQ=255)
                                                                                                                                 gCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATc                  >  1:102743/1‑150 (MQ=255)
                                                                                                                                      aCGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTg             <  2:633543/150‑1 (MQ=255)
                                                                                                                                       cggcggGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGc            <  1:914288/150‑1 (MQ=255)
                                                                                                                                          cGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCCCCGGTACTACGCTCCGTTGCTgg         >  2:725826/1‑150 (MQ=255)
                                                                                                                                                 cgcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGta  <  1:301633/150‑1 (MQ=255)
                                                                                                                                                 cgcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGta  >  2:505431/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTA  >  NC_000913_3_bme_pgi/526649‑526943

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: