Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7920AC25.4% 99.5 / 39.1 69S1178R (AGC→CGC) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base C (6/11);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCA  >  NC_000913_3_bme_pgi/526645‑526934
                                                                                                                                                    |                                                                                                                                              
tGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAg                                                                                                                                               >  1:114198/1‑149 (MQ=255)
  gCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGTGACTGCCGCTGGCGCTTCTCAGCGCAGACGGCCCAACGCCGTGGcgcg                                                                                                                                             >  2:96517/1‑149 (MQ=14)
    ttGCGCAGTGCGGGCTGACGGTGGAGCACCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCCCTGGCACTTATATGAGTAGTCGGCAATCCGGCGGAGGGcgcg                                                                                                                                           >  2:353375/1‑149 (MQ=14)
    ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                           <  2:633563/149‑3 (MQ=255)
      cggaagtgagtgCTGACGGCGGAAGCCCTGGCCAGGCCAGTTGAGCCGGACTATACCCCTGCGCCAAAAGCTCCTCTTTCTCACGGG‑GACCACCGGTGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGa                                                                                                                                         <  2:318206/145‑1 (MQ=2)
          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                                                                                                          <  1:493193/149‑1 (MQ=255)
                                     ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                          <  2:658771/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)
                                          cAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                           <  2:35071/111‑1 (MQ=255)
                                          cAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                           >  1:35071/1‑111 (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                                                                                                  >  2:158651/1‑133 (MQ=255)
                                                             cacCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                  <  1:158651/133‑1 (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)
                                                                                                                                                    |                                                                                                                                              
TGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCA  >  NC_000913_3_bme_pgi/526645‑526934

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