Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7960CG35.3% 97.9 / 91.4 102A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (29/37);  new base G (12/24);  total (41/61)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 106 total aligned reads displayed.

CTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTAT  >  NC_000913_3_bme_pgi/526648‑526944
                                                                                                                                                    |                                                                                                                                                    
cTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgcgc                                                                                                                                                      <  2:845763/149‑1 (MQ=255)
cTTGCGCAGTGCGGGCTGACGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgcgc                                                                                                                                                      <  2:928158/149‑1 (MQ=255)
 gtGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                                     <  1:891427/148‑1 (MQ=255)
  tGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCgg                                                                                                                                                    >  1:635480/1‑149 (MQ=255)
   gcgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGTGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGa                                                                                                                                                   <  2:300084/149‑1 (MQ=255)
    cgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACCGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGag                                                                                                                                                  <  1:1161600/149‑2 (MQ=255)
    cgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGCGCCGGAATACACACCGGCGCGAAAAGCTGATCTTTATCACTGCGACCACCGGGGGCTGCCGCGGGCGCTTAGCAGCGCTGACGGCAATTCGGCGGGGAGCGCGGaa                                                                                                                                                  >  2:1312664/1‑149 (MQ=14)
     gcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAat                                                                                                                                                 <  1:1218874/149‑1 (MQ=255)
     gcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAat                                                                                                                                                 >  1:314283/1‑149 (MQ=255)
      cAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatt                                                                                                                                                >  2:1459056/1‑148 (MQ=255)
      cAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAata                                                                                                                                                <  2:701368/149‑1 (MQ=255)
       aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                               <  1:319025/149‑1 (MQ=255)
        gggCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatg                                                                                                                                              <  2:203478/147‑1 (MQ=255)
        gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatg                                                                                                                                              <  1:918377/149‑1 (MQ=255)
        gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatg                                                                                                                                              <  2:658406/149‑1 (MQ=255)
          gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                            <  2:342149/149‑1 (MQ=255)
           cGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgatg                                                                                                                                           >  1:196837/1‑149 (MQ=255)
              gCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                                                                        <  2:930503/149‑1 (MQ=255)
              gCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAAt                                                                                                                                        <  2:1153575/149‑1 (MQ=255)
                tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGCCCACCTGGGACTGCCGCTGGCGCTTTTCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                      >  2:1382550/1‑149 (MQ=255)
                tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                                                      <  1:1163698/149‑1 (MQ=255)
                 gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCCGCGCGAAAAGCTCATCTTTAGCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                                     <  1:1154004/149‑1 (MQ=11)
                      tGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                                <  2:244330/149‑1 (MQ=255)
                      gggAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                                <  2:1409991/148‑1 (MQ=255)
                       ggAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGTAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAcc                                                                                                                               <  1:1304750/149‑1 (MQ=255)
                          gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGAACACCCGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGAAATACGGCGAGGCGCGGGGAGTATCATGAATGGGGCAACCAGc                                                                                                                            >  1:365972/1‑149 (MQ=11)
                           cAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCt                                                                                                                           <  1:426792/149‑1 (MQ=255)
                            aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACGGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCtt                                                                                                                          <  2:1381655/149‑1 (MQ=255)
                             aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTa                                                                                                                         <  1:977433/149‑1 (MQ=255)
                             aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTa                                                                                                                         <  1:1093605/149‑1 (MQ=255)
                             aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTa                                                                                                                         <  1:479102/149‑1 (MQ=255)
                              cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGAGGAATGGGGCAACCAGCTTaa                                                                                                                        >  1:218282/1‑149 (MQ=255)
                              cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGGGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTaa                                                                                                                        >  1:327586/1‑149 (MQ=255)
                               tgtccagccagGTGGAGACGGAATACACCCCGGCGCGAAAATCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAAt                                                                                                                       <  2:109714/146‑1 (MQ=17)
                               tGGCCAGACAGGTGGAGGCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAAt                                                                                                                       <  1:16551/149‑1 (MQ=255)
                                ggCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATg                                                                                                                      <  1:288141/149‑1 (MQ=255)
                                 gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATga                                                                                                                     >  1:25743/1‑149 (MQ=255)
                                 gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                     <  1:448180/149‑1 (MQ=255)
                                  ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCCCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                    <  2:1281516/149‑1 (MQ=255)
                                  ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                    <  1:437383/149‑1 (MQ=255)
                                  ccAGACAGGTGGAGCCGCAATACACACCGGCGCGAAAAGCTCATCTTTAACACTGCCACCACCGGGGACTGCCGCAGGCGCATATCACCGAAGACGGCTATACGGCGTGTAGCGCGTAATATGATGAATGGGGAAACCAGCTTAATtag                                                                                                                    >  2:683236/1‑149 (MQ=14)
                                    aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                  >  2:1360445/1‑149 (MQ=255)
                                       cAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGaa                                                                                                               <  1:273780/149‑1 (MQ=255)
                                        aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTGAATGAGGAGAAc                                                                                                              >  2:951758/1‑149 (MQ=255)
                                         gggggAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAcc                                                                                                             <  2:138086/146‑1 (MQ=255)
                                          gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                            >  1:605752/1‑149 (MQ=255)
                                          gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                            <  2:1405251/149‑1 (MQ=255)
                                          gTGGAGCCGGAAGACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCGTAATGAGGAGAAccc                                                                                                            >  2:1014949/1‑149 (MQ=11)
                                             gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                         >  2:854032/1‑149 (MQ=255)
                                             gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATAAGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGct                                                                                                         >  2:291756/1‑148 (MQ=255)
                                              aGCCGGAATTCACACCGGCGCTAAAAGCTCATCTTTATCACTGCGACGACCGGGGACTGGCGCTGGCGGTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcat                                                                                                        <  2:500704/149‑1 (MQ=11)
                                              aGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcat                                                                                                        <  1:763670/149‑1 (MQ=255)
                                              aGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcat                                                                                                        <  2:1356297/149‑1 (MQ=255)
                                               gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGCCCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatc                                                                                                       >  2:1378302/1‑149 (MQ=255)
                                               gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatc                                                                                                       >  1:99743/1‑149 (MQ=255)
                                               gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatc                                                                                                       <  2:432782/149‑1 (MQ=255)
                                                ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATAAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatca                                                                                                      <  1:452321/149‑1 (MQ=255)
                                                 cGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcac                                                                                                     <  2:604194/149‑2 (MQ=255)
                                                   gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGCCGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                                   >  2:1276912/1‑149 (MQ=255)
                                                    aaTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctg                                                                                                  <  2:196837/149‑5 (MQ=255)
                                                       acacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGTGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAt                                                                                               <  1:1474620/149‑1 (MQ=255)
                                                         acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAGGTGTATcc                                                                                             >  1:301096/1‑148 (MQ=255)
                                                         acacCGGCGCGAAAAGCTCATCTTTAACACTGCGACCACCGGGGACAGCCCCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATTATGTGTATCa                                                                                             >  1:907541/1‑149 (MQ=255)
                                                           acCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGc                                                                                           <  2:1218433/149‑1 (MQ=255)
                                                           acCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGATTATCAGCGAAGACGGCAATACGGCGTTGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGc                                                                                           <  2:260663/149‑1 (MQ=255)
                                                               gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTa                                                                                       >  1:1290155/1‑149 (MQ=255)
                                                                        gCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGTCGGCAATACGGCGTGGAGGGCGGAATATGATGAATGGGGCAACCAGCTTACTGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGGCTGcc                                                                              >  1:129385/1‑149 (MQ=255)
                                                                         cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCa                                                                             <  1:1375019/149‑1 (MQ=255)
                                                                         cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCa                                                                             <  2:426601/149‑1 (MQ=255)
                                                                            aTCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGTGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAggg                                                                          <  2:750881/149‑1 (MQ=255)
                                                                               tttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGAGGAGAGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGACAACCCGCATCATGTGTATCAGCCGTACAGTCTGCCAGCgcac                                                                       >  2:1268987/1‑148 (MQ=17)
                                                                                   tCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATCAATGGGGCAACCAGCAAAATGAGGAGAACCCGCATCATGTCTATCAGACGTACCGTCTGCCAGGgcagcagc                                                                   >  2:254393/1‑149 (MQ=17)
                                                                                      cTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatg                                                                >  1:672586/1‑149 (MQ=255)
                                                                                          gACCACCGGGTACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATATGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                            <  2:251205/149‑1 (MQ=255)
                                                                                          gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCGTAATGAGGAGAACCCGCATCCTGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAggatga                                                            >  1:342299/1‑149 (MQ=255)
                                                                                          gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGGGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                            <  2:948734/149‑1 (MQ=255)
                                                                                           accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGATCCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                           >  2:754973/1‑149 (MQ=255)
                                                                                                cGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATc                                                      >  2:46484/1‑149 (MQ=255)
                                                                                                     aCTGACGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGc                                                 <  2:635480/149‑1 (MQ=255)
                                                                                                      cTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATTATGAAGGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCt                                                <  1:1295719/149‑1 (MQ=255)
                                                                                                      atGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGTCAGGGCAGCAGCATGATGAGGAATCAGGGCt                                                <  2:1267634/148‑1 (MQ=255)
                                                                                                       tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                                  <  2:320676/146‑1 (MQ=255)
                                                                                                       tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                                  >  1:320676/1‑146 (MQ=255)
                                                                                                       tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTg                                               >  1:678841/1‑149 (MQ=255)
                                                                                                       tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGACAGGGCACCAGCATGATGAGGAATCAGGGCTg                                               >  1:656169/1‑149 (MQ=255)
                                                                                                         ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGATTATGATGAATGGGGCAACCAGCTTAATGAGTAGAACCCGCATTATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                             <  2:907541/149‑1 (MQ=255)
                                                                                                         ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGCGGAATCAGGGCTGtc                                             >  2:307197/1‑148 (MQ=255)
                                                                                                         ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                             >  2:302171/1‑149 (MQ=255)
                                                                                                           gCTGGCGGTTAACAGCGAAGACGGCAATACGGCGTGGAGAGCGGAATATGATGAATGGGGCAAACAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACt                                           <  2:972810/149‑1 (MQ=255)
                                                                                                            gtGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCTGAATATGATGATTGGGGCAACCAGGTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACta                                          <  1:131903/148‑1 (MQ=255)
                                                                                                             tGGTGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGCACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACtat                                         >  2:1288919/1‑149 (MQ=255)
                                                                                                                       cAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCACGGCAGCAGCatgat                                                              >  1:164691/1‑118 (MQ=255)
                                                                                                                         gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGCGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAGGATAAGGAATCAGGGCTGTACTATAACCGTCACCgg                             >  1:948401/1‑149 (MQ=255)
                                                                                                                            aaGACGGCAATACGGCGTGGCGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGACCCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCGGTACTATAAACGTCACCGGtac                          >  2:1428161/1‑149 (MQ=255)
                                                                                                                                      tACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGt                <  2:314283/149‑1 (MQ=255)
                                                                                                                                       aCGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGtt               <  2:605752/149‑1 (MQ=255)
                                                                                                                                        cGGCGTTGAGCGCGGAATATGATGAATGGTGCAACCAGCTTAATGAGGAGAACCAGCATCATGTGTATCAGCCGTACCTTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTg              <  1:1268987/149‑1 (MQ=255)
                                                                                                                                           cGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAg           <  1:564927/149‑1 (MQ=255)
                                                                                                                                              ggAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCTTGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAgggg        <  1:291756/149‑1 (MQ=255)
                                                                                                                                                    cgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGtat  >  1:786206/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTAT  >  NC_000913_3_bme_pgi/526648‑526944

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: