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

CAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGG  >  NC_000913_3_bme_pgi/526654‑526941
                                                                                                                                                   |                                                                                                                                             
cAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAata                                                                                                                                             <  1:212661/149‑1 (MQ=255)
 agggggggCGGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAatat                                                                                                                                            <  2:291766/144‑1 (MQ=255)
  gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATatg                                                                                                                                           <  1:289832/149‑1 (MQ=255)
    gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTatgat                                                                                                                                         <  1:597874/149‑1 (MQ=255)
      gggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTAtgatga                                                                                                                                       >  1:506269/1‑149 (MQ=255)
       ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATGaa                                                                                                                                      >  1:44123/1‑149 (MQ=255)
       ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGCATatgatgga                                                                                                                                      >  1:174397/1‑147 (MQ=255)
          tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                                                   <  2:547527/149‑1 (MQ=255)
          ggagggTGGAGCAACTGGCCAGACAGGTTGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                   <  2:169075/145‑1 (MQ=255)
           ggcGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATGACTGCGACCACCGGGGACTGCAGCTGGCGCTGATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                                  <  2:461223/147‑1 (MQ=1)
           gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                                  <  2:44123/149‑1 (MQ=255)
           gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGAGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                                  <  2:461244/149‑1 (MQ=11)
            gctgTGGGGGAACAGGCCAGGCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTACCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATgggg                                                                                                                                 <  2:132325/146‑1 (MQ=14)
                tGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                             <  2:197468/149‑1 (MQ=255)
                  gAGCAACTGGCCAGACAGGTGGAGCGGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGCCCGCTGGAGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAA‑ACGCCGTGGCGCGGGGAGTATGATGAATGGGGCAGCCa                                                                                                                           >  1:95913/1‑149 (MQ=11)
                    gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGc                                                                                                                         >  1:112025/1‑149 (MQ=255)
                     cAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGCAGACGGCAAT‑ACGGCCTGGCGCCCGAACTAGGAACGATGGCGCAACCAGCt                                                                                                                        >  2:59598/1‑149 (MQ=14)
                     cAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCt                                                                                                                        >  1:493830/1‑149 (MQ=255)
                       aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCATT‑ACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTa                                                                                                                      >  1:572190/1‑149 (MQ=255)
                       aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATAAATGGGGCAACCAGc                                                                                                                         <  1:43448/146‑1 (MQ=11)
                       aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATAAATGGGGCAACCAGc                                                                                                                         >  2:43448/1‑146 (MQ=255)
                           gCCGGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACGGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                  <  2:191955/149‑1 (MQ=255)
                            ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                 >  1:444786/1‑149 (MQ=255)
                             cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTAGCACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                >  1:134987/1‑149 (MQ=11)
                                 cAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGaa                                                                                                            <  2:444786/149‑1 (MQ=255)
                                  aGGTGGAGCCGAAAGACACACCGGCGCGCAAAGCCCATCTTTATCACTGCGACCAGCGGGGACTGCCGCTGGCGCTTATCAGCGAATAGGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAc                                                                                                           <  2:621593/149‑1 (MQ=17)
                                   ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAcc                                                                                                          <  2:609156/149‑1 (MQ=255)
                                    gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCCCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATGACAGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                         >  1:347817/1‑149 (MQ=11)
                                    gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                         >  2:555850/1‑149 (MQ=255)
                                    gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                         <  1:146978/149‑1 (MQ=255)
                                     tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCACTAACGGCG‑GTGGCTGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                        >  1:20111/1‑149 (MQ=14)
                                         gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATTGTGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatc                                                                                                    >  1:537964/1‑149 (MQ=255)
                                           cGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACTGCGTGGGGCGCGGGATATGATGAATCGGGCAAACGGCTTACTGAGGGGGAACCGcatcag                                                                                                  >  2:429106/1‑148 (MQ=14)
                                               aTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATgtgt                                                                                              >  1:12405/1‑149 (MQ=255)
                                                   acacCGGCGCGAAAAGCTGATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGCCCGACGGCGTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                          <  2:149458/149‑1 (MQ=18)
                                                   acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                          <  1:357101/149‑1 (MQ=255)
                                                     acCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGc                                                                                        >  2:502890/1‑149 (MQ=255)
                                                          cgcgAAAGGCCTATCTTCATCACGGCGTCCACCGGGGACTGCCGCTGGCGCTTAGAAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGTATATGATGAATGGCGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTAc                                                                                   <  2:522039/149‑1 (MQ=18)
                                                          cgcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTAc                                                                                   <  2:12405/149‑1 (MQ=255)
                                                             gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGACAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCACACAGCTTAATGAGGAGAACCCGCATCAGATGTATCAgccgtgccgt                                                                                >  1:531521/1‑149 (MQ=18)
                                                              aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                                <  2:520016/148‑1 (MQ=255)
                                                              aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                                >  1:520016/1‑148 (MQ=255)
                                                               aaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCt                                                                              <  1:629957/149‑1 (MQ=255)
                                                                    tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                         >  1:256324/1‑149 (MQ=255)
                                                                    tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                         <  1:324106/149‑1 (MQ=255)
                                                                    tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                         >  1:329682/1‑149 (MQ=255)
                                                                      attcttATCACTGCGACCACCGGGGACTGCCGCTGGCGGTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAggg                                                                       <  1:84045/145‑1 (MQ=255)
                                                                        cTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgca                                                                     <  1:688902/149‑1 (MQ=255)
                                                                         tttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGAGGAGCGCGGAATATGACGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATAAGCGGTACCGTCTgccaggccag                                                                    >  2:518437/1‑149 (MQ=17)
                                                                                 tgggACCAGCGGGGACTGCCGCTGGCTCCTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatga                                                            <  2:568598/146‑1 (MQ=255)
                                                                                       caccGGGGACTGCCGCTGGCGCTTATCAGCGAAAACGGCACT‑ACGGGGTGGAGCGCAGAATATGAGGAATGGGGCGGCCAGCCTAATGAGGAGAACCCGCATCATGTGAATAAGCCGCACCGTCTGCCAGGCCAGAAGCAGAATCAGGa                                                      >  2:189716/1‑149 (MQ=17)
                                                                                          cGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATCTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATc                                                   >  1:475228/1‑149 (MQ=255)
                                                                                              gACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                               <  1:173660/149‑1 (MQ=255)
                                                                                                       tGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACtat                                      >  1:495183/1‑149 (MQ=255)
                                                                                                           gcTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTCTAAcc                                  >  1:22698/1‑149 (MQ=255)
                                                                                                           gcTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAAAGGGCAAATTAGCGTTATGAGAAGAACCCGCATCATGTGTATTAGCCGTACCGTCTGCCAGGGCCGCAACATGATGAGGAATCCGGGCTGTACTATaagc                                  >  2:426553/1‑147 (MQ=17)
                                                                                                                  aGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCg                           <  2:613482/149‑1 (MQ=255)
                                                                                                                    cGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGt                         >  2:370898/1‑149 (MQ=255)
                                                                                                                     gAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGta                        <  1:502890/149‑1 (MQ=255)
                                                                                                                        gACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtacta                     <  2:399478/149‑1 (MQ=255)
                                                                                                                        gACGGAAAT‑ACGGCGGGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtacta                     <  2:174397/149‑1 (MQ=255)
                                                                                                                             cAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCATGGCTGTACTATAACCGTCACCGGTACTACGATc                >  2:619418/1‑149 (MQ=255)
                                                                                                                               aT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCg              >  2:139293/1‑149 (MQ=255)
                                                                                                                                  aCGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGtt            <  1:152405/149‑1 (MQ=255)
                                                                                                                                  aCGGCGTGGAGCGCGAAACATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGtt            <  1:152374/149‑1 (MQ=255)
                                                                                                                                      catGGAGCGCGGAATATGATGAATGGGGCAACCAGCTCAATGAGGCGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAg        <  1:604113/147‑1 (MQ=255)
                                                                                                                                        tGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAggg      <  1:59598/149‑1 (MQ=255)
                                                                                                                                            gctcgGAATATGCTGAATGGGGCAACCAGCTTAATGGGGAGAACCCGGATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTATTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCgg  <  2:223798/146‑1 (MQ=17)
                                                                                                                                                   |                                                                                                                                             
CAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGG  >  NC_000913_3_bme_pgi/526654‑526941

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