Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,8000AG21.4% 110.9 / 28.0 73E1180E (GAA→GAG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (20/35);  new base G (5/10);  total (28/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATAT  >  NC_000913_3_bme_pgi/526657‑526946
                                                                                                                                               |                                                                                                                                                  
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:802656/149‑1 (MQ=255)
                         ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                      <  1:434156/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                                                                                                                                      <  1:469844/122‑1 (MQ=255)
                                    gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                                      >  2:469844/1‑122 (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)
                                                                               gCGGCCAACGGGGACTGCCGCTGGCGCTGATCAGCGAAGACGGCAATGCGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGTTTAACGAGGCGACCCCGCCTCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat                                                                <  2:770655/149‑1 (MQ=17)
                                                                                 gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGGACGCAGGAGATGATGAATGGGGCAACCAGCTTGATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCGGCAGCatgatta                                                              >  1:196031/1‑147 (MQ=17)
                                                                                 gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                              >  2:156386/1‑69 (MQ=255)
                                                                                 gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                              <  1:156386/69‑1 (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)
                                                                                                                                        gcgcGGAATATGAGGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCgg       <  2:424768/149‑1 (MQ=255)
                                                                                                                                           cgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATATCCTTCACCGGTACTACGATCCGTTGCAGGGGCGGtat    <  1:995298/149‑1 (MQ=255)
                                                                                                                                           cgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGtat    <  2:670890/149‑1 (MQ=255)
                                                                                                                                            ggAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGtata   >  2:652868/1‑149 (MQ=255)
                                                                                                                                             gAATATGATGAATGGTGCAACCACACTACTGAGTGCAACGCGCATAATGTGTATCAGCCGTACCGTCTGCCACGGCAGCAGTATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTatat  <  2:272433/149‑1 (MQ=17)
                                                                                                                                               |                                                                                                                                                  
TGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATAT  >  NC_000913_3_bme_pgi/526657‑526946

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