Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7920AC27.1% 159.8 / 64.6 107S1178R (AGC→CGC) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (29/49);  new base C (10/19);  total (39/68)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 107 total aligned reads displayed.

GTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGG  >  NC_000913_3_bme_pgi/526644‑526936
                                                                                                                                                    |                                                                                                                                                
gTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAg                                                                                                                                                 >  1:144715/1‑150 (MQ=255)
gTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAg                                                                                                                                                 >  1:148332/1‑150 (MQ=255)
     ttGCGCAGGGCGGGCTGACGGTGGAGCAACTGGCCAGCCAGGTGGCGCCGGAATACACACCGGCGCGCAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCgg                                                                                                                                            <  2:745357/150‑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:769753/1‑150 (MQ=255)
          cAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                       >  1:763339/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)
                      tcGGTGGAGCAACTGGCCAGACAGGTTGAGCCGGAATACACACCGGCTCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGc                                                                                                                           <  1:595271/149‑1 (MQ=11)
                      aCGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGc                                                                                                                           <  1:354173/150‑1 (MQ=255)
                       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                                                                                                          <  2:445242/150‑1 (MQ=255)
                                       cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                          <  2:317771/150‑1 (MQ=255)
                                       cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                          <  1:494382/150‑1 (MQ=255)
                                        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)
                                             ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                    >  2:158859/1‑150 (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)
                                                                    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)
                                                                            tcgcATCTTTATCACTGCGACCAGCGGGGACTGCCGCTGGCGCTTATAAGCGAAGACGGCTATACGGCGTGGAGCGCGGATTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCa                                                                     <  1:856737/147‑1 (MQ=17)
                                                                             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                                                  <  1:5985/150‑1 (MQ=255)
                                                                                               accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAgg                                                  <  2:239176/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: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)
                                                                                                                                               cGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCCCCGGTACTACGCTCCGTTGCTgg  >  2:725826/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGG  >  NC_000913_3_bme_pgi/526644‑526936

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: ‑