Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7960CG27.9% 111.3 / 47.1 68A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (21/28);  new base G (10/9);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGG  >  NC_000913_3_bme_pgi/526648‑526941
                                                                                                                                                     |                                                                                                                                                 
cTTGCGCAGGGCGGGCTGAGGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAgcgc                                                                                                                                                   <  2:112969/149‑1 (MQ=255)
 ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTCCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGcgcg                                                                                                                                                  <  1:548869/149‑1 (MQ=255)
   gcgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGTGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGa                                                                                                                                                <  1:323990/149‑1 (MQ=255)
      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                                                                                                                         >  2:43448/1‑146 (MQ=255)
                             aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGCGCGGGGAGTATGATAAATGGGGCAACCAGc                                                                                                                         <  1:43448/146‑1 (MQ=11)
                                 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)
                                                                                                                                                     |                                                                                                                                                 
CTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAT‑ACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGG  >  NC_000913_3_bme_pgi/526648‑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: ‑