Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,8000AG23.1% 121.3 / 27.0 65E1180E (GAA→GAG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (23/27);  new base G (6/9);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATA  >  NC_000913_3_bme_pgi/526655‑526945
                                                                                                                                                  |                                                                                                                                                 
aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAgtat                                                                                                                                                <  1:202799/149‑1 (MQ=255)
  tGCGGGCTGACCGTGGAGCAACTGGCCAGGCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatga                                                                                                                                              <  1:226576/149‑1 (MQ=11)
  ggCGGGCTGACGGTGGAGCAACTGGCCAGACACGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GAGCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                              <  1:22321/148‑1 (MQ=255)
         tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                                                       <  2:716864/149‑1 (MQ=255)
         tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                       >  1:413472/1‑149 (MQ=255)
           aCGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                                     >  1:446590/1‑149 (MQ=255)
             ggTGGAGCAAATGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCa                                                                                                                                   >  1:465848/1‑149 (MQ=255)
              gTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCGTCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCaa                                                                                                                                  >  1:475294/1‑149 (MQ=255)
               gggaggaACTGGCCAGACAGGTGAAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                                 <  1:716215/148‑1 (MQ=11)
                ggAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGCCGGCAATACGGCGTGGAGCGCGGAATATGATGACTGGGGCaaac                                                                                                                                >  2:280020/1‑147 (MQ=255)
                   gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGc                                                                                                                             <  1:387312/149‑1 (MQ=255)
                    cAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCCGCt                                                                                                                            <  2:590809/149‑1 (MQ=11)
                          gCCAGAGNGGTGGCGGCGGGATACACGCCGGNGAGAAAAGCTCTTCTTNAGCNCTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCGGCGATGAGGGCNATAAGGCGTGGAGCGCGGAATATGATGTATGTGGCAACCAGCTTAATga                                                                                                                      <  2:44589/149‑1 (MQ=255)
                          gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                      <  2:378334/149‑1 (MQ=255)
                           ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                     <  2:658771/149‑1 (MQ=255)
                           ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                     <  1:493193/149‑1 (MQ=255)
                            cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                    >  1:604746/1‑149 (MQ=255)
                            cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                    <  1:421059/149‑1 (MQ=255)
                            cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACGGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                    <  2:165111/149‑1 (MQ=255)
                             acaCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATCTTTATTCCTTC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                   <  1:353252/147‑1 (MQ=255)
                             aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                   >  2:579323/1‑149 (MQ=255)
                             aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                   <  2:297279/149‑1 (MQ=255)
                                 aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAc                                                                                                               <  2:186003/149‑1 (MQ=255)
                                 aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAATCAGCTTCACGAGGAGAAc                                                                                                               >  1:677962/1‑149 (MQ=255)
                                  ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAcc                                                                                                              >  2:134246/1‑149 (MQ=255)
                                   gTGGAGCCGGAATACACGCCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                             <  2:603953/149‑1 (MQ=255)
                                   gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                                 >  2:559058/1‑129 (MQ=255)
                                   gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                                 <  1:559058/129‑1 (MQ=255)
                                   gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                             >  2:305804/1‑149 (MQ=255)
                                    tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATCCGACGCGCTGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                            >  1:602085/1‑149 (MQ=17)
                                        gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGGGTATGATGAATGGGGCAACAGGCTTAATGAGGAGAACCCGcatc                                                                                                        >  1:94766/1‑149 (MQ=11)
                                          cGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcat                                                                                                      <  2:161852/149‑1 (MQ=255)
                                          cGGAAGACACACCGGGGTGAAAAGCTCATCTTTATCACGGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGCCAACCAGCTTAATGAGGAGAACCCGcatcac                                                                                                      <  1:346563/149‑2 (MQ=11)
                                            gAATACACACCGGCGCGAAACGCTCATCTTTATCACTGC‑GACCACCAGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                                    >  1:492338/1‑149 (MQ=255)
                                            gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGTGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                                    <  2:84902/149‑1 (MQ=255)
                                            gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                                    <  1:228305/149‑1 (MQ=255)
                                              aTACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATgtgt                                                                                                  <  1:134246/149‑1 (MQ=255)
                                                  acacCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                              >  1:198318/1‑149 (MQ=255)
                                                   cacCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                             <  1:158651/133‑1 (MQ=255)
                                                   cacCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                             >  2:158651/1‑133 (MQ=255)
                                                    acCGGCGCGAAAAGCTCATCTTTATCACAGCAGGGCTTC‑CATACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGc                                                                                            >  1:571649/1‑149 (MQ=18)
                                                     ccGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCAGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGcc                                                                                           <  2:492338/149‑1 (MQ=255)
                                                        gcgcGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTa                                                                                        >  1:106466/1‑149 (MQ=255)
                                                                   gcATCTTTATCACTGC‑GATCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                             <  1:520648/148‑1 (MQ=255)
                                                                     aTCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATGTGATGAATGGGGCAATCAGCTTAACGAGGAGAACCCGGATCATGTGTATGAGCCGTACGGTCTGCAAggg                                                                           >  2:650608/1‑149 (MQ=18)
                                                                       cTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGTTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgca                                                                         <  1:751161/149‑1 (MQ=255)
                                                                        tttATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACTGCGTTGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCGTCATGTGTATTAACCGTACCGTCTGCCAGGgcag                                                                        >  2:321040/1‑149 (MQ=17)
                                                                                tGC‑GACCACCGGGGACTGCCGCAGGCGCTTATCCGCGGAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCCACCAGCTTATTGAGGAGGACCCACATCATGGGTAACAGCCGTACCAGCTGCCAGGGCAGCGGCatga                                                                >  2:693478/1‑149 (MQ=255)
                                                                                    gACCACCGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                             <  2:114198/149‑1 (MQ=255)
                                                                                       caccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGa                                                          >  1:1296/1‑149 (MQ=255)
                                                                                            gggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAg                                                     >  1:670255/1‑149 (MQ=255)
                                                                                            gggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGCCAGCAGCATGATGAGGAATCAg                                                     >  2:430801/1‑149 (MQ=255)
                                                                                             ggACTGCCGCTGGGGCATATGAGCGAAGTAGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAGCGGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAgg                                                    <  2:262800/149‑1 (MQ=17)
                                                                                              gACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                                   <  2:198318/149‑1 (MQ=255)
                                                                                                      cTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACta                                           >  1:421005/1‑149 (MQ=255)
                                                                                                       tGGCGTTTATCACCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACtat                                          <  1:118338/149‑1 (MQ=255)
                                                                                                       tGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACtat                                          >  1:33126/1‑149 (MQ=255)
                                                                                                          cgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAc                                       <  1:531610/149‑1 (MQ=255)
                                                                                                           gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGCATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                      >  1:427783/1‑149 (MQ=255)
                                                                                                               aTCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAACGTCa                                  >  2:266422/1‑149 (MQ=255)
                                                                                                                           ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGa                      >  1:39911/1‑149 (MQ=255)
                                                                                                                                 aCGGCTTGGAGCGCGGAATATGAAGAATGGGGCAACCAGCTTACTCAGTCGTACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGtt                <  1:722874/149‑1 (MQ=17)
                                                                                                                                  cGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTg               <  1:167885/149‑1 (MQ=255)
                                                                                                                                    gCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGAGCAGCAGCAGGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCa             <  1:149224/149‑1 (MQ=255)
                                                                                                                                               ggAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGtata  <  1:280020/149‑1 (MQ=255)
                                                                                                                                               ggAACATGATGAATGGGGCAACCAGCTTAATGAGGAGACCCCGCATCATGTTTACCAGCCGTACCGTGTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGtata  <  1:398093/149‑1 (MQ=14)
                                                                                                                                                  |                                                                                                                                                 
AGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATA  >  NC_000913_3_bme_pgi/526655‑526945

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