Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7960CG26.2% 112.9 / 44.7 65A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (23/25);  new base G (6/11);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCA  >  NC_000913_3_bme_pgi/526649‑526934
                                                                                                                                                    |                                                                                                                                          
ttGCGCAGTGCGGGCTGACGGTGGAGCACCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCCCTGGCACTTATATGAGTAGTCGGCAATCCGGCGGAGGGcgcg                                                                                                                                           >  2:353375/1‑149 (MQ=14)
  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                                                                                                          <  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)
                                      cAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                           >  1:35071/1‑111 (MQ=255)
                                      cAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                           <  2:35071/111‑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                                                                                                                      <  1:559058/129‑1 (MQ=255)
                                         gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                      >  2:559058/1‑129 (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)
                                                                                                                                                    |                                                                                                                                          
TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGC‑GACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCA  >  NC_000913_3_bme_pgi/526649‑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: ‑