Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,296,0610CT29.0% 95.2 / 40.9 71intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (27/22);  new base T (12/8);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAT  >  NC_000913_3_bme_pgi/4295914‑4296202
                                                                                                                                                    |                                                                                                                                             
tCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACCCCAGACCGCGCGATAGCCTGATTTAGCGTGATTTTGTAGCTCG‑GATACACCGTTTACGCCGCATCCGACATCAACg                                                                                                                                              >  2:12512/1‑149 (MQ=11)
 cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTACGTCT‑GATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                             <  2:756753/149‑1 (MQ=255)
     tAAGTCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                                                                                                         <  1:961/149‑1 (MQ=17)
      tagccGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCG‑GATAGAGCGTTTACGCCGCATCCGACATCATAGCCTGag                                                                                                                                        >  2:686070/5‑148 (MQ=11)
       aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAACGTCATTTTGTAAGACA‑GAAAAGGCGTTTACGCCGCATCCGACATCAATGCCTGaag                                                                                                                                       >  1:51665/1‑147 (MQ=11)
        ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                      >  1:755101/1‑149 (MQ=17)
         gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCg                                                                                                                                     >  1:51581/1‑149 (MQ=32)
           gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGGGTTTACGCTGCATCCGCCAGCAAGGCCTGTTGCGAc                                                                                                                                   >  1:148116/1‑149 (MQ=11)
              tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGCCATCCATCCCTGATGCGACGCt                                                                                                                                >  1:334556/1‑149 (MQ=255)
               aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                                                                                               >  2:71926/1‑149 (MQ=17)
               aTGCCGCAGCCGACATCAGCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCtt                                                                                                                               <  1:226303/149‑1 (MQ=17)
                  ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCTTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCCACATCAATGCCTGATGCGACGCTTGcc                                                                                                                            >  2:529341/1‑149 (MQ=255)
                  ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGCTGCAGCGCTTGcc                                                                                                                            >  1:372077/1‑149 (MQ=255)
                   cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcg                                                                                                                           >  1:756753/1‑149 (MQ=35)
                    gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                          <  1:71926/149‑1 (MQ=21)
                     cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGTTTTTGTAGCTCG‑GCTAAGACGTTTACCCCGCATCCGCCATCAATGCCTGATGCGACGCTTGCcgcg                                                                                                                         >  2:383189/1‑149 (MQ=11)
                      aTCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                        <  2:442160/149‑1 (MQ=38)
                      aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                        <  1:561766/149‑1 (MQ=21)
                      aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                        >  1:602759/1‑149 (MQ=38)
                         cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                     <  2:51581/149‑1 (MQ=38)
                          tgaatgaaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                    <  2:320089/143‑1 (MQ=255)
                            cATCCACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                  <  1:661658/149‑1 (MQ=255)
                             aTCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCATACAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCC‑GATAAGGCGTTTACGCCGCATCCGACATCCATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                 <  2:567292/149‑1 (MQ=11)
                              tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                >  1:129084/1‑149 (MQ=25)
                              tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                >  1:605072/1‑149 (MQ=38)
                              tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCCTCTTATCAg                                                                                                                >  2:339528/1‑149 (MQ=37)
                                aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGAATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                              >  1:774787/1‑149 (MQ=255)
                                 aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                             >  2:338282/1‑149 (MQ=37)
                                     cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                         >  1:620981/1‑147 (MQ=17)
                                      tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                        >  1:203802/1‑149 (MQ=38)
                                       gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGTTTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAGGCGACGCTTGCCGCGGCTTATCAGGCCTatctt                                                                                                       >  1:619832/1‑145 (MQ=11)
                                        aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                      <  1:399177/149‑1 (MQ=39)
                                          gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                    <  1:398597/149‑1 (MQ=39)
                                                    aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATa                                                                                          >  1:214816/1‑149 (MQ=37)
                                                     cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                        >  1:663018/1‑118 (MQ=25)
                                                     cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                        <  2:663018/118‑1 (MQ=25)
                                                      gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                >  1:253890/1‑125 (MQ=40)
                                                      gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                <  2:253890/125‑1 (MQ=40)
                                                      gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                        >  1:746682/1‑149 (MQ=37)
                                                       cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGACACCGCCATAGcc                                                                                       >  2:172752/1‑149 (MQ=255)
                                                         tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                              >  2:346810/1‑124 (MQ=25)
                                                         tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                              <  1:346810/124‑1 (MQ=25)
                                                         tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGGGATTTTGTCGGTCGCGCCAATG‑GTTGACTCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCGGGCCTACGCCAGACAGCGCAATAGCCTg                                                                                     >  2:561766/1‑149 (MQ=11)
                                                          cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                              <  1:159491/123‑1 (MQ=40)
                                                          cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                              >  2:159491/1‑123 (MQ=39)
                                                                aGGCCTACGGCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                              <  1:71251/149‑1 (MQ=38)
                                                                  gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                <  1:16543/113‑1 (MQ=255)
                                                                     tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGACATCGATATAGCCTGCTTTAGCGTAAt                                                                         >  2:125149/1‑149 (MQ=11)
                                                                         ccagccagCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCTAGCCAGCGCATTAGCCTGCTTTAGCGTCATTTTg                                                                     >  2:490438/1‑149 (MQ=14)
                                                                          cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                    >  1:725891/1‑149 (MQ=38)
                                                                           aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTCCGCAGCACCGGAGATCAACGCCTGATGCGACCCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                   <  1:382548/149‑1 (MQ=11)
                                                                             aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                                                                                                <  2:68821/118‑1 (MQ=40)
                                                                             aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAgg                                                                 <  1:291049/149‑1 (MQ=39)
                                                                                gcgcAACCGCATGATTCAGCGTTGATGTGTAGGGCG‑GATCAGGCATATACGCCGCATCCGACGTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCg                                                              <  2:762941/149‑1 (MQ=11)
                                                                                      tAGCCTGATTTAGCTTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAg                                                        <  1:783984/149‑1 (MQ=38)
                                                                                         ccTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCg                                                     >  1:654977/1‑149 (MQ=38)
                                                                                                  gCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGGAATAGCCTTAAATATCGTGACTTTCGAGGTCGGATCAGGCGTTTATGCCg                                            >  1:77016/1‑149 (MQ=11)
                                                                                                  gCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                        >  1:350981/1‑105 (MQ=255)
                                                                                                  gCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                        <  2:350981/105‑1 (MQ=255)
                                                                                                   cGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                           <  1:260672/53‑1 (MQ=255)
                                                                                                      gATTTTGTGGGGCG‑GTCAGGGCGTTTACGCCCCACCAGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATc                                        <  2:340186/149‑1 (MQ=11)
                                                                                                               ggTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCaa                               >  1:476436/1‑149 (MQ=39)
                                                                                                                gTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTCATTTTGTGGGTCGGATAAGGCGTTTATGCCGAATACGACATCacc                              >  2:609971/1‑147 (MQ=255)
                                                                                                                 tCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACg                             <  1:338282/149‑1 (MQ=38)
                                                                                                                   g‑gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTACCGTGATTTTGTAGGTCGGATCATGCATTAAAGCCTCCTCCGACATctgaccc                           >  2:361783/1‑143 (MQ=14)
                                                                                                                      aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                               <  1:99849/63‑1 (MQ=255)
                                                                                                                        aaGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAt                       <  2:214816/149‑1 (MQ=38)
                                                                                                                                tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACCCt               >  2:534265/1‑149 (MQ=38)
                                                                                                                                           cGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACGCCTGATTCGCCGCTTAACGCGTCtt    >  2:768942/1‑149 (MQ=255)
                                                                                                                                            gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTa   >  1:365285/1‑149 (MQ=25)
                                                                                                                                             gcATGAACGCCCGATCCGGCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAt  <  2:322977/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCG‑GATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAT  >  NC_000913_3_bme_pgi/4295914‑4296202

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