Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,8000AG28.1% 139.8 / 72.6 98E1180E (GAA→GAG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (23/46);  new base G (8/19);  total (31/67)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 106 total aligned reads displayed.

GCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTAC  >  NC_000913_3_bme_pgi/526653‑526949
                                                                                                                                                   |                                                                                                                                                     
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)
  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)
    tGCGGGCGGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACCCCGGCGCGTAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTTTCAGCGAAGAGGGCAATACGGCGTGGCGCGGGGAGTatgat                                                                                                                                                 <  1:662954/150‑1 (MQ=11)
    tGCAGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCCCCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                                 >  2:218261/1‑150 (MQ=255)
     gCGGGCTTAAGGTGGAGCAACTGGCCAGACTGGAGGAGCCTGAATAGACACCGGCGCGAAAAGATCATCTTTATCACTGCGACCACCGGTGACTGCCGCTGGAGCTTATCAGCTAAGACTGCAATACGGCGTGGAGCGCGGAATatgatg                                                                                                                                                <  1:67956/150‑1 (MQ=14)
     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)
             aCGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGc                                                                                                                                        <  1:354173/150‑1 (MQ=255)
              cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCa                                                                                                                                       <  1:245840/150‑1 (MQ=255)
              cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCa                                                                                                                                       <  1:819933/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)
                      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:808825/150‑1 (MQ=255)
                              cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                       <  2:445242/150‑1 (MQ=255)
                              cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                       <  1:494382/150‑1 (MQ=255)
                              cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                       <  2:317771/150‑1 (MQ=255)
                              cAGACAGGTGCTGGCGGAATACACACCGGCGCGAAAAGCTCATCTATATCACAGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGTTGAATGGGGCAACCTGCTTAATgagga                                                                                                                       <  1:629537/150‑1 (MQ=11)
                               aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgaggag                                                                                                                      <  2:918505/150‑1 (MQ=255)
                                   aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACAGCGACCACCGGGGACTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAcc                                                                                                                  <  1:205489/150‑1 (MQ=255)
                                   aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCAGCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAcc                                                                                                                  <  1:208719/150‑1 (MQ=255)
                                    ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                                 >  2:158859/1‑150 (MQ=255)
                                       ggAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                              <  2:808541/150‑1 (MQ=255)
                                       ggAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                              <  2:367478/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)
                                                           cgcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTAcc                                                                                          >  1:816933/1‑150 (MQ=255)
                                                               aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCt                                                                                      <  1:479546/150‑1 (MQ=255)
                                                               aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCCTGGAGCGCGGAATATCATGAATGGGGCAACCAGCTAAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCt                                                                                      >  2:462304/1‑150 (MQ=255)
                                                                 aaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGc                                                                                    <  2:70783/150‑1 (MQ=255)
                                                                    cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGAGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                                 <  1:416674/150‑1 (MQ=255)
                                                                      cATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAggg                                                                               >  1:633543/1‑150 (MQ=255)
                                                                            tATCACAGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGAAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcag                                                                         <  2:758560/150‑1 (MQ=255)
                                                                              tCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcagca                                                                       >  2:386080/1‑150 (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)
                                                                                                            gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCATCCAGCTTTATGAGGAGAACGCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCg                                         <  2:763339/150‑1 (MQ=255)
                                                                                                            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:558388/1‑150 (MQ=255)
                                                                                                                     cGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGta                                >  1:918620/1‑150 (MQ=255)
                                                                                                                     cGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGta                                >  1:906248/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)
                                                                                                                                          gatcgtggAATATGATTAATGTGGCAACCAGCATAATGAGGAGAACCCGCATCATGTGTATCAGCCGTTCCGGCTGCCAGGGTAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCg           <  2:247578/144‑1 (MQ=17)
                                                                                                                                             cgcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGta        <  1:301633/150‑1 (MQ=255)
                                                                                                                                             cgcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGta        >  2:505431/1‑150 (MQ=255)
                                                                                                                                                   atatGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTAc  >  1:661336/1‑150 (MQ=255)
                                                                                                                                                   atatGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGAAGAAGCATGATGAGGAATCAGGGCTGTACTATACCCGTCCCCGGTACGACGCTCGGTTGTAGGGGCGGTATATTAc  >  1:561999/1‑150 (MQ=17)
                                                                                                                                                   |                                                                                                                                                     
GCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTAC  >  NC_000913_3_bme_pgi/526653‑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: ‑