Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7920AC27.9% 67.5 / 42.0 71S1178R (AGC→CGC) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (19/30);  new base C (7/12);  total (28/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCG  >  NC_000913_3_bme_pgi/526644‑526940
                                                                                                                                                    |                                                                                                                                                    
gttgCTTGCGCAGGGCGGGCGGACGGTGGGGCAACGGGTCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGa                                                                                                                                                      <  2:787083/146‑1 (MQ=14)
  ggCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgc                                                                                                                                                    <  1:502228/149‑1 (MQ=255)
     gtgggcAGGGCGGGCTGACGGTGGAGCAACTGGCCAGAAAGGTGGAGCCGGAATAAACTCGGGCGTGAACAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                                 <  2:111264/145‑3 (MQ=11)
       gctagccgcgggcggACGGGTGAGCAACTGCGGACACAGGTGGAGCCGGAATAGACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGa                                                                                                                                               <  2:408606/142‑1 (MQ=11)
       gcgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGa                                                                                                                                               <  1:600667/149‑1 (MQ=255)
            gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGaggagg                                                                                                                                          >  1:272462/1‑144 (MQ=255)
             tGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                         >  1:943933/1‑149 (MQ=255)
             tGCGGGCGGACGGGGGAGGAACTGGCGAGACGGGTGGAGCGGGAATACACACCGGGGCGAAAAGCTCATCTTGATCACTGCGACCACCGGGGACTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                         <  2:285369/149‑1 (MQ=17)
              gcgggcggACGGTGGAGCAACTGGCCTGACAGGTGGAGCCGTAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGGAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                        <  2:636791/149‑1 (MQ=255)
                   gtGAGGGGGGGGCAACTGGCCAGGGAGGTGAGGCCGGAATACAGGCCGGCGGGAAGAGCGTATCTTTATCACGGGGACCAGCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGGGTGGAGCGTTTAATATGATGAATg                                                                                                                                   <  2:374070/148‑1 (MQ=9)
                   cTGAGGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATg                                                                                                                                   <  2:350196/149‑1 (MQ=255)
                    tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTCTGATGAATgg                                                                                                                                  >  2:738207/1‑149 (MQ=255)
                    tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                                                  <  2:6914/149‑1 (MQ=255)
                    tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATAAGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                  <  1:1007147/149‑1 (MQ=255)
                    tGACGGGGGAGCAACTGGCCAGACAGGTGGAGGCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCAGCTGGCGCTTATCAGCGAAGACGGAAATACGGCGTGGAGCGAGGAATATGAGGAATgg                                                                                                                                  <  2:497269/149‑1 (MQ=17)
                     gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg                                                                                                                                 <  2:819763/149‑1 (MQ=255)
                          tGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACGCCGGCGCGAAAAGCTCATCTTTATCACTGCGATCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATCGGGCAAc                                                                                                                            <  1:930596/149‑1 (MQ=11)
                          tGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                            <  2:777054/149‑1 (MQ=255)
                                aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCtt                                                                                                                      >  1:756986/1‑149 (MQ=255)
                                 aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTa                                                                                                                     >  1:862678/1‑149 (MQ=255)
                                  ccggCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCCGGCGCTAATCAACGAAGCCAGCAATCCGGCCGGGACCGCGGACACCGAAGAATGCGGCCACCAGCtgca                                                                                                                    >  2:204647/3‑146 (MQ=14)
                                   tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCCGGGACTGCCGCCGGCGCTTATCAGCGAAGACGGCAATACCGCGTGGAGCGCGGAATAAGATGAAAAGGACACCCAGCTTaac                                                                                                                   >  2:279971/1‑148 (MQ=17)
                                    ggCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATg                                                                                                                  <  1:247643/149‑1 (MQ=255)
                                     gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATga                                                                                                                 <  1:837086/149‑1 (MQ=255)
                                     gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                 >  1:448607/1‑149 (MQ=255)
                                     gCCAAACAGGAGGAGCCGGAACACCCCCCAGCGCGAAAAGATCATCTTTAACACTGCAACCACCGGGGACTGCCGCTGGCGCGTATCAGCCAAGCCGCCAACACGGCGTGTAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                 <  2:422456/149‑1 (MQ=14)
                                      ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                <  1:434156/149‑1 (MQ=255)
                                      ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                <  1:802656/149‑1 (MQ=255)
                                      ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                <  1:74071/149‑1 (MQ=255)
                                       cAGACAGGTGGAGCCGGAATACACTCCGGCGCGAAAAGCTCATCTTTATAAATACGAACACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACTAGCTTAATgagg                                                                                                               >  1:27025/1‑149 (MQ=17)
                                       cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                               <  1:39479/149‑1 (MQ=255)
                                       cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                               <  1:939188/149‑1 (MQ=255)
                                       cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAAAATGATAAATGAGGCAACCAGCTTCAtaaag                                                                                                               >  1:308203/1‑145 (MQ=17)
                                        aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                              <  2:675542/149‑1 (MQ=255)
                                        aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                              >  1:670890/1‑149 (MQ=255)
                                        aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAACACGGCAATACGGCGTGGCGCGGGGAGTCGGATGAATGGAGCAACAGGCTTAATgcgga                                                                                                              >  2:663706/1‑149 (MQ=11)
                                        aGACAGGGGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                              <  2:382255/149‑1 (MQ=255)
                                         gACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCCGGGACTGCCGCTGGCGCTTACCAGCGAAGACGGCAACACGGCCTGGAGCGAGGAATATGATCAATGAGGCGAGCAGGTGATTaaggag                                                                                                             >  2:100942/1‑149 (MQ=14)
                                              gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                        >  1:381791/1‑149 (MQ=255)
                                                 gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                                >  2:469844/1‑122 (MQ=255)
                                                 gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                                <  1:469844/122‑1 (MQ=255)
                                                   gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcacg                                                                                                   >  1:504598/1‑147 (MQ=255)
                                                   gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGACGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTGAATGAGGAGAACCCGcatc                                                                                                   >  1:956031/1‑149 (MQ=255)
                                                      ggAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtg                                                                                                <  2:984938/149‑1 (MQ=255)
                                                         aTACACACCGGCGCGAAAAGCTCATCTTTAGCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATCCGGCGTGGAGCGAGGAATACGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATgtgt                                                                                             <  2:345492/149‑1 (MQ=255)
                                                           acacacCGACGCGAAAAGCCCATCTTTATCACTGCGACCACCGGGGACTGCCGCGGGCGCATATCAGCGGATACGGCAATCCGCCGGGGGGGGCGGAATATGAGAAACGGGGAAGACACCCTACGGCGGAGAACCCGCAACACGCGTag                                                                                           >  2:204933/1‑148 (MQ=255)
                                                             acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                         >  1:424768/1‑149 (MQ=255)
                                                                   gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTa                                                                                   >  1:958380/1‑149 (MQ=255)
                                                                     gcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCCAACACGGCAATACGGCGTTGAGCCCGGACTATGATGAATGGGGCACCCAGCTTAGCGAAGCAATCCGGCACCATGAGTACCAGCCGTTcc                                                                                 >  2:668776/1‑149 (MQ=255)
                                                                       gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAAGGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                               >  1:225075/1‑149 (MQ=255)
                                                                             cgcaccTTTATCACTGCCACCACCGGGGACGCCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGGAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCa                                                                         <  2:691143/144‑1 (MQ=17)
                                                                              tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                        >  1:608845/1‑149 (MQ=255)
                                                                               cATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                            >  1:278492/1‑96 (MQ=11)
                                                                               cATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                            <  2:278492/96‑1 (MQ=11)
                                                                                     tATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAGTACGGCGTGGCGCGCGGATTATGATGAATGGGGCAACAAGCTTACTGAGGAGCACCCGCGTCAGGTGTATCAGCCGTACCGTCTGCCAGGgcagca                                                                 >  2:564059/1‑149 (MQ=14)
                                                                                       tCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                                 <  1:265710/67‑3 (MQ=255)
                                                                                       tCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                                 >  2:265710/1‑65 (MQ=255)
                                                                                            gCGGCCAACGGGGACTGCCGCTGGCGCTGATCAGCGAAGACGGCAATGCGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGTTTAACGAGGCGACCCCGCCTCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat                                                          <  2:770655/149‑1 (MQ=17)
                                                                                              gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGGACGCAGGAGATGATGAATGGGGCAACCAGCTTGATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCGGCAGCatgatta                                                        >  1:196031/1‑147 (MQ=17)
                                                                                              gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                        <  1:156386/69‑1 (MQ=255)
                                                                                              gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                        >  2:156386/1‑69 (MQ=255)
                                                                                              gACCACCGGGGACTGCCGCTAGGGCTTATCACGGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                        <  2:423000/149‑1 (MQ=255)
                                                                                               cccgccGGGGACTGCCGCGGGGGCTTTTCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGTAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                       <  1:450422/145‑1 (MQ=17)
                                                                                               agcaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                       <  1:920044/147‑1 (MQ=255)
                                                                                                     ggggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCa                                                 <  2:958380/149‑1 (MQ=255)
                                                                                                      gcggcTGCCGCTGGTGCGTATCGGCGAAGACGGCAATACGGAGTGGAGCGCGGAACATGATGAATGGGGCAACCCGCTTAAGGAGTAGAGCCCGCATCAGGTGTATCAGCCGGACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAg                                                <  2:637863/145‑1 (MQ=17)
                                                                                                        gACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                              <  1:677041/149‑1 (MQ=255)
                                                                                                        gACGGCCGCGGGCGCGTAGCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                              <  1:222679/149‑1 (MQ=255)
                                                                                                          cTGCCGCTGGCGCTTATCAGCGAATACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGGGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCt                                            <  2:399003/149‑1 (MQ=255)
                                                                                                                   gcgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATaa                                   <  1:150482/149‑1 (MQ=255)
                                                                                                                     gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                 >  2:804526/1‑149 (MQ=255)
                                                                                                                                      gCAATACGGCGTGGGGCGCGGAATATGATGAAGGGGGCAACCAGCGTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGAt                <  1:100942/149‑1 (MQ=255)
                                                                                                                                         atccGGCGTGTAGCGCGGAATATGATGAATGGGTCAAGCTGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGGACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCg             <  1:259856/146‑1 (MQ=17)
                                                                                                                                          tACGGCGTGGACCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGt            <  1:638560/149‑1 (MQ=255)
                                                                                                                                                   gAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGc   <  2:538254/149‑1 (MQ=255)
                                                                                                                                                    aGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCg  <  2:448607/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCG  >  NC_000913_3_bme_pgi/526644‑526940

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