Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7920AC34.9% 106.7 / 86.6 107S1178R (AGC→CGC) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (29/40);  new base C (19/18);  total (48/59)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-02
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 110 total aligned reads displayed.

TGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAG  >  NC_000913_3_bme_pgi/526645‑526935
                                                                                                                                                   |                                                                                                                                                
tGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAg                                                                                                                                                 <  2:850328/149‑1 (MQ=255)
  gCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGACAGACAGGTGGAGCCGGAATACAAACCGGCGCGAAATGCCCATCTTTATCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcg                                                                                                                                               >  1:917491/1‑149 (MQ=11)
    ttGCGCAGTGCGGGCTGTCGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                             <  1:78331/149‑3 (MQ=255)
    ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCAAGAAAGGTGGACCCGGAATACACACCGGCGCGAAAAGCTCATCTTTCTCACTGTGTGCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                             <  1:448822/149‑3 (MQ=14)
    ttGCGCAGGGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                             <  2:1293534/149‑1 (MQ=255)
       cgcAGCGCGGGCGGACGGTGGAGCAACGGGCCAGCCAGGTGGAGCCGGAATACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGGCCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGag                                                                                                                                          <  2:797832/149‑2 (MQ=11)
        gcAGTGCGGGCTGACGGTGGAGCAACTGGCGAGACAGGGGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCAACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAGACGGCGTGGAGCGCGGAat                                                                                                                                         >  2:416537/1‑149 (MQ=255)
        gcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAat                                                                                                                                         >  1:958815/1‑149 (MQ=255)
          aGTGCGGGGTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAgtat                                                                                                                                       <  2:252798/149‑1 (MQ=255)
          aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                       >  1:1014435/1‑149 (MQ=255)
          aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                       >  1:1246750/1‑149 (MQ=255)
            tGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                     <  1:1210178/149‑1 (MQ=255)
            tGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGACGGCAGCTCTTCTCAGCGAAGACGGCACTACGGCGTGGAGCGCGGAATgtga                                                                                                                                     >  2:139826/1‑149 (MQ=17)
            ggCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatga                                                                                                                                     <  1:60418/148‑1 (MQ=255)
                 gCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAtgatgagt                                                                                                                                >  1:571993/1‑147 (MQ=255)
                   tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                                              <  1:518451/149‑1 (MQ=255)
                   tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                              <  2:344600/149‑1 (MQ=255)
                    taCGGTTGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                             <  1:745023/148‑1 (MQ=11)
                    ggcGGGGGAGCGAGTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                             <  2:253363/147‑1 (MQ=11)
                    gagggTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                             <  2:233685/146‑1 (MQ=255)
                    gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGTGGCTTCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                             <  1:705629/149‑1 (MQ=11)
                     aCGGTGGAGCAACTGGCCAGAAAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgggg                                                                                                                            <  1:1163754/149‑1 (MQ=255)
                         tGGAGCAACTGGCCAGACAGGTGGAGCCGTAATACACACCGGCGTGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                        <  1:747125/149‑1 (MQ=255)
                         gcgAGCAACTGGCCAGACAGCTGGAGCCGTACTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACTGGGGTCTGCCGCTGTCGCTTATCTGCGAAGACGGCAATACGGCGTGGTGCGCGGAATATGCTGAATGGGGCAAc                                                                                                                        <  1:1087601/147‑1 (MQ=14)
                           gAGCAACTGGCCAGACAGGTGGAGCCGGAATACACGCCGGCGCGAAAAGCTCATCTTTATCACTACGACCACCGGGGACTGCCGCTGGCGCCTATCCGCGAAGACGGCAATACGGAGGGGCGAGGGGAGCATGATGACTGGGGCAAACa                                                                                                                      >  2:849206/1‑149 (MQ=11)
                           gAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCa                                                                                                                      >  1:571595/1‑149 (MQ=255)
                             gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGc                                                                                                                    >  1:398939/1‑149 (MQ=255)
                             gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCACCGGATCCGGGAATACGGCGTGGCGGGGGGGGTATGATGAATGGGGCAGCCAga                                                                                                                    >  1:251438/1‑148 (MQ=11)
                               aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCtt                                                                                                                  >  1:970449/1‑149 (MQ=255)
                                 cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGGGAAAAGATAATATTTATGACTGCGACCAAAGGGGACTGCAGCAGGAGCTTATCAGCGAAGACGGCAATAAGGCGTGGAGGGCGGAATATGATGAATGGGGGAAGCAgactaa                                                                                                                >  2:1284789/1‑144 (MQ=14)
                                 cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTaa                                                                                                                <  1:479862/149‑1 (MQ=255)
                                 cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTaa                                                                                                                >  1:1247181/1‑149 (MQ=255)
                                 cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTaa                                                                                                                <  1:823654/149‑1 (MQ=255)
                                    gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCGGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGTCAACCAGCTTAATga                                                                                                             >  1:113616/1‑149 (MQ=11)
                                    gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATga                                                                                                             >  1:105077/1‑149 (MQ=255)
                                     ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgag                                                                                                            <  1:1116028/149‑1 (MQ=255)
                                     ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgag                                                                                                            <  1:1062377/149‑1 (MQ=255)
                                     ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                            <  2:989505/149‑1 (MQ=255)
                                      cAGATAGGTGGAGCCGGAATACACACCGGCGGGAAGGGCGCATGTGTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                           <  1:33137/149‑1 (MQ=1)
                                      cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                           >  1:540311/1‑149 (MQ=255)
                                      cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                           >  2:715163/1‑149 (MQ=255)
                                       aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCCACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGAGAAGCGAGGGGAGGATGAGGAATGGGGCAACCAGCTTAATgaggc                                                                                                          >  1:245533/1‑148 (MQ=11)
                                         aCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATCCGGCGTGTCGCGGGGAGTATGATGAATGGGGCAACCAGCTTCAGGAGggga                                                                                                        >  1:512514/1‑149 (MQ=11)
                                          cAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGaa                                                                                                       <  1:619554/149‑1 (MQ=255)
                                           aGGTGGAGCCGGAATACACACCGGCGCGCAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCACACCCTTCCATACGGCGTGGCGCGGGGAGGATGATGAATGGGGCAACCAGCTTAATGAGGAGAAc                                                                                                      >  1:1300868/1‑149 (MQ=11)
                                             gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGTACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                    >  2:702906/1‑149 (MQ=11)
                                             gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                    <  1:441950/149‑1 (MQ=255)
                                              tGGAGCCGGAATAGCCAGCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTCGCGCTTATCAGCGAAGACGTCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                   <  1:86565/149‑1 (MQ=17)
                                              tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                   <  2:571595/149‑1 (MQ=255)
                                              tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                   <  1:1204326/149‑1 (MQ=255)
                                                gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAACAGCTTAATGAGGAGAACCCGca                                                                                                 >  1:333330/1‑149 (MQ=255)
                                                  gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGGGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatc                                                                                               >  1:147574/1‑149 (MQ=255)
                                                   ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                        <  2:214073/123‑1 (MQ=255)
                                                   ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                        >  1:214073/1‑123 (MQ=255)
                                                      gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                             <  2:9492/115‑1 (MQ=11)
                                                      gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                             >  1:9492/1‑115 (MQ=255)
                                                      gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                           <  2:970449/149‑1 (MQ=255)
                                                         tACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTa                                                                                        <  1:918581/149‑1 (MQ=255)
                                                          acacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAt                                                                                       <  2:958815/149‑1 (MQ=255)
                                                           cacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGTAGGCGGCAATACGGCGTGGAGCGGGAGATATGATGTATGGGGCAGCCAGCGTAATGGAGAGAACCCGCACCATGTTTAta                                                                                      >  2:1035034/1‑148 (MQ=255)
                                                            ccacCGGCCCGAAAAGCTCATCTTTATCACTGCGACAACCGGGGACTGCCGCCGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                     <  2:249065/148‑1 (MQ=255)
                                                                 ggCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGt                                                                                >  1:805541/1‑149 (MQ=255)
                                                                  gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                             >  2:996219/1‑119 (MQ=255)
                                                                  gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                             <  1:996219/119‑1 (MQ=255)
                                                                      caaaaGGTCGTCTTTAGCCCGGCGAACACCGGGACCTGCCGCTGGCGCCGATCAGCGAAGACGGCAATAAGGCGTGGAGGTCTGGGTGTGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                           <  2:1208434/148‑1 (MQ=14)
                                                                       aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGCCGGCAGTACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGGCGTACCGTc                                                                          >  1:780696/1‑149 (MQ=255)
                                                                        aaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCt                                                                         <  1:211953/149‑1 (MQ=255)
                                                                         aaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTg                                                                        >  1:556896/1‑149 (MQ=255)
                                                                          tgCTCATCTTTATCACGGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCTTCATGTGTATCAGCCGTACCGTCTGc                                                                       <  2:225706/148‑1 (MQ=255)
                                                                            cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcat                                                                                             >  2:57265/1‑125 (MQ=255)
                                                                             tCATGTTTATCACTGCGCCCGCCGGGGACTGCCGCTGGCGATTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                    <  2:562994/149‑1 (MQ=255)
                                                                                 cTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgat                                                                                                                                    >  1:56032/1‑81 (MQ=11)
                                                                                  tttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                       <  1:866339/77‑1 (MQ=255)
                                                                                  tttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                       >  2:866339/1‑77 (MQ=255)
                                                                                   ttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGag                                                                                                                                          >  2:644394/1‑72 (MQ=11)
                                                                                   ttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGag                                                                                                                                          <  1:644394/73‑2 (MQ=11)
                                                                                        aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑Cat                                                        <  1:1309571/149‑1 (MQ=255)
                                                                                         cTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCAGTACCGTCTGCCAGGGCAGCAG‑Catg                                                       >  1:1115360/1‑149 (MQ=255)
                                                                                             gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CAtgatga                                                   <  1:1134810/149‑1 (MQ=255)
                                                                                             gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CAtgatga                                                   <  2:1014435/149‑1 (MQ=255)
                                                                                                  ccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAAt                                              <  1:139826/149‑1 (MQ=255)
                                                                                                    tgggTCTGCCGGTGGCGGTTATCAGGGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCa                                            <  2:589119/148‑1 (MQ=255)
                                                                                                     gggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAg                                           >  2:1015361/1‑149 (MQ=255)
                                                                                                      ggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAga                                          >  2:856886/1‑148 (MQ=255)
                                                                                                       gACGGCCGCTGGCGGTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAggg                                         <  1:1231058/149‑1 (MQ=255)
                                                                                                          tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGCACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAGTCAGGGCTg                                      >  2:124109/1‑149 (MQ=255)
                                                                                                          tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGTGCTg                                      >  2:672681/1‑149 (MQ=255)
                                                                                                                tGGCGCGTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACtat                                <  1:588906/149‑1 (MQ=255)
                                                                                                                 ggCGCTTAGCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACtata                               <  1:1093980/149‑1 (MQ=255)
                                                                                                                  gcgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTCATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACTATaa                              >  1:342308/1‑149 (MQ=255)
                                                                                                                  gcgcTTATCAGCGAAGACGGCAAGACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACTATaa                              <  2:344076/149‑1 (MQ=255)
                                                                                                                    gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACGCGCATCATTTTTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACTATAAcc                            <  1:1246822/149‑1 (MQ=255)
                                                                                                                    gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCCGGAG‑CATGATGAGGCAGCAGGGCTGTACTATAAcc                            >  1:2070/1‑149 (MQ=255)
                                                                                                                    gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAC‑CATGATGAGGAATCAGGGCTGTACTATAAcc                            >  2:954286/1‑149 (MQ=255)
                                                                                                                          cAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGTAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGGCGCCATCCT‑TGGAGTCAGGGCTGTACTATAACCGTCAcc                      <  2:1280940/149‑1 (MQ=17)
                                                                                                                           aGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACTATAACCGTCACCg                     <  2:1161093/149‑1 (MQ=255)
                                                                                                                            gCGAGGACGGAAATAGTGCGTGGATCGCGTAATAGGAAGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAGGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACTATAACCGCCACCgg                    <  2:541793/149‑1 (MQ=17)
                                                                                                                                      cAATATGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATc          >  2:741317/1‑149 (MQ=255)
                                                                                                                                             gcttGGAGCGCGGAATATGATGATTGGTGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCa   <  1:336954/146‑1 (MQ=255)
                                                                                                                                              cGTGGAGCGCGGAATATGATGAAGTGGGCAACCAGCTTAATGAGGAGACCCTGCTCCCCCGCTATTAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAg  <  1:640808/149‑1 (MQ=14)
                                                                                                                                                   |                                                                                                                                                
TGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAG‑CATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAG  >  NC_000913_3_bme_pgi/526645‑526935

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