Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0390TG44.1% 85.5 / 15.0 62N108T (AAC→ACC) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (7/26);  new base G (26/0);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 1.12e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677890‑2678187
                                                                                                                                                      |                                                                                                                                                     
ccaGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGt                                                                                                                                                       <  2:609840/150‑1 (MQ=255)
 caGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGC‑GGGCCACCGACCGAGGCAGA‑GGgttt                                                                                                                                                    >  2:559568/1‑151 (MQ=255)
      gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGACCGAGGCGGA‑GGgttt                                                                                                                                                    >  1:1134562/1‑147 (MQ=255)
      gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTcc                                                                                                                                                  <  2:395237/149‑1 (MQ=255)
      gctATTGCGTCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTATCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGGATCCTGCGCCATTGGCAGGAAATAAAGCATGCCGTTGCGGCGGGGCCACTGCCCGAGGCAGA‑GGTTTTcc                                                                                                                                                  <  1:31920/149‑1 (MQ=255)
       ctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGACCCAGGCAGAGGGTTTTcccc                                                                                                                                                >  2:30126/1‑149 (MQ=255)
       ctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGc                                                                                                                                                <  1:819871/150‑1 (MQ=255)
        ttttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGA‑GGGTTTCCGcc                                                                                                                                               >  1:55997/1‑150 (MQ=255)
            gcgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCCCGGCAGAGGGTTTTCccccccg                                                                                                                                            >  1:308903/1‑147 (MQ=255)
              gcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCGCCGACCGCGGGAGG‑GGGTTTTCGCCAccctc                                                                                                                                         >  1:133077/1‑146 (MQ=255)
                cGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGAGCCGCCTAACACAGCAGA‑AGATTTCCGCCaacgcttac                                                                                                                                      >  1:321171/1‑143 (MQ=255)
                  ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGGGGCAGA‑GGTTTTCCGCCACTGTC‑AGccc                                                                                                                                   >  2:306675/1‑151 (MQ=255)
                     ttAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGAGGCAGA‑GGGTTTTCGCCACCGTC‑AGCCCtt                                                                                                                                 >  1:440643/1‑150 (MQ=255)
                        gCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTCCCCCACTGTC‑AGCCCTTccac                                                                                                                             >  2:523759/1‑151 (MQ=255)
                        gCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTcca                                                                                                                              <  1:497597/150‑1 (MQ=255)
                        gCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGA‑GGGTTTTCGCCCCCGTC‑AGCCCCTccac                                                                                                                             >  2:166662/1‑151 (MQ=255)
                              aaaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTCCGCCCCTGTC‑AGCCCTTCCACCAGAc                                                                                                                        >  2:1028533/1‑150 (MQ=255)
                               aaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCAACCGAAGGAGA‑AGAATTTTTCCACTGTC‑GGGCCCTCCACCCCACtt                                                                                                                      >  1:372169/1‑151 (MQ=255)
                                 cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTaa                                                                                                                    <  2:535474/151‑1 (MQ=255)
                                 cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGGCCCAGGCAGA‑GGGTTTTCCCCCCTGTC‑AGCCCGCCCACCAGACt                                                                                                                       >  2:452936/1‑148 (MQ=255)
                                  aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAAc                                                                                                                   <  2:751071/151‑1 (MQ=255)
                                   aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGGGGGAGA‑GGGTTTTCCCCGCCGTC‑AGCCCCTCCACCCGACt                                                                                                                       >  2:69027/1‑146 (MQ=255)
                                    gttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGCACCAGGGAGA‑GGGTTTTCGCCCCTGTC‑AGCCCCTCCCCCCAACttca                                                                                                                    >  2:847730/1‑146 (MQ=255)
                                      tgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGAGGAAGG‑GGGTTTCCGCCCCCGTC‑AGCCCCTCCACCAAAATTAAcctc                                                                                                                >  2:1239694/1‑147 (MQ=255)
                                       gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTg                                                                                                              <  1:390248/151‑1 (MQ=255)
                                       gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGAGCCCCCGACCCAGGCCGA‑TGGTTTTCGCCGCTGTC‑CGCCCCTccacc                                                                                                                            >  1:1206162/1‑137 (MQ=255)
                                         tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGAGGGAGA‑GGGTTTTCGCCACCGTC‑AGCCCCTCCACCAGACTTTACTCCTGAt                                                                                                            >  2:460119/1‑151 (MQ=255)
                                         tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCAGA‑GGGTTTTCGCCACCGTC‑AGCCCTTCCCTCAGAATTAACTcc                                                                                                                >  2:634168/1‑147 (MQ=255)
                                                 aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTTCGCCCCCGTC‑AGCCCCTCCCCCCGACCTAACTCCCGGTAAAAcg                                                                                                      >  1:577291/1‑149 (MQ=255)
                                                      cAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTTACCCAGGCAGA‑GGGTTTCCGCCACTGTC‑CGCCCTTCCCCCAGACTTAAATCCTGATAAAAcccgc                                                                                                   >  2:603547/1‑147 (MQ=255)
                                                          gCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCCAGGGAGA‑GGGTTTTCGCCCCCGTC‑CGCCCCTCCACCCCACCTAACTCCTGATAAAACCCGCGCAccccc                                                                                           >  1:871870/1‑151 (MQ=255)
                                                          gCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCAccccc                                                                                           <  2:704038/151‑1 (MQ=255)
                                                           cccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCa                                                                                          <  1:987387/151‑1 (MQ=255)
                                                           cccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCCCGGGAGA‑GGGTTTCCGCCACTGTC‑AGCCCCTCCCCCCGCCCTAAATCCCCATAAACCCCGCGCccccccc                                                                                          >  1:234319/1‑150 (MQ=255)
                                                            ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGAATGGGGGAGA‑GGGTTTTCGCCACTGTC‑GGGTCCTCCACCCAACTTAACTCCTGTGAAAACCCGCGCAacccc                                                                                           >  1:881717/1‑149 (MQ=255)
                                                            ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAg                                                                                         <  1:260785/151‑1 (MQ=255)
                                                                 gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGCGGCAGA‑GGGTTTTCCCCACTGTC‑AGCCCCTCCACCCAACTTAACTCCCGATAAAACGCGCGCACCCCCacgtc                                                                                      >  2:778817/1‑145 (MQ=255)
                                                                  ccTGTAAGCAACGTTGGGCATCCTGCGCCATGTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACGGACCGAGCCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGc                                                                                   <  2:1105573/151‑1 (MQ=255)
                                                                   cTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACCGAGGGAGA‑GGGTTTTCGCCCCTGTC‑AGCCCTTCCCCCCGACTTAACTCCCTATAAAACGCGCGCccccccccctc                                                                                      >  1:749805/1‑142 (MQ=255)
                                                                      tAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCCCGGCCGA‑GGGTTTCCGCCCCTCTC‑ACCCCTTCCCCCAGACCTTACTCCCGATAAACCCCCCCCCCCCCCCATTCACCg                                                                                  >  1:158729/1‑148 (MQ=255)
                                                                      tAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACCGCGGCGGA‑GGGTTTTCGCCCCTTTC‑CGCCCCTCCCCCCGACCTTACTCCTTATAAAACCCGCGCACCCCCCCTTCAGccg                                                                                 >  1:2361/1‑147 (MQ=255)
                                                                           aaCGTTGGGCATCCTGCGCCATTTGTAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCg                                                                          <  1:1152468/151‑1 (MQ=255)
                                                                              gTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTGGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGa                                                                         <  2:787781/149‑1 (MQ=255)
                                                                                  ggCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTTACCCGAGCAGA‑GGGTTTCCCCCACTCTC‑AGCCCTTCCCCCAGACCTAAATCCTGATAAAACCCGCGCCCCCCCCCTTCAACGGGACAGCGAGTCa                                                                     >  2:1228995/1‑149 (MQ=255)
                                                                                   gCATCCTGCCCCATTTGCAGGTAATCAATCATGCCGGTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCg                                                                  <  2:144025/151‑1 (MQ=255)
                                                                                        cTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGGCCGAGGCAGA‑GGGTTTTCTCCACCTTC‑CGCCCCTCCACCAGACTTAAATCCTGATAAACCCCCCGCACCCCCCGGTCCGCGGCACCGCGGGGCAGCGGAg                                                               >  2:899991/1‑149 (MQ=255)
                                                                                              cATTTGCAGGTAATCAATCATGCCGTTGCGGCGGCGCCACTCATTAAGGCAGA‑GGCTTTCCCCGACTCTCTGGCCCTTCC‑CCAGCGTTCACACCTGATTAATCCCGTGCACCCCCAACTCAACGGCCAGGCGAGTCCGcagatatttg                                                          >  2:1158529/1‑139 (MQ=255)
                                                                                                tttGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGGGGA‑GGGTTTTCGCCCCTCTC‑AGCCCCTCCCCCAGACCTAACTCCTGATAAAACCCGCGCCCCCCCCCTTCAGCGGCACGGCGGGGCAGCGGAGCCTCGTcccc                                                     >  1:1060357/1‑151 (MQ=255)
                                                                                                  tGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCt                                                    <  2:1163329/150‑1 (MQ=255)
                                                                                                     aGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCa                                                <  2:317966/151‑1 (MQ=255)
                                                                                                        tAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGCAGA‑GGGTTTTCGCCACTGTC‑AGCCCCTTCACCCGACTTAAAACCCGATAAACCGCGCGCACCCCCCCTTCAGCCGCACGGCGAGGCAGGGGAGATTCGTTCCCTTcccctc                                             >  2:812239/1‑147 (MQ=255)
                                                                                                          atcaatcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcgc                                           <  1:325493/151‑1 (MQ=255)
                                                                                                          atcaatcaTGCCGTTGCGGCGGGGCCACCGACCGCGGCAGA‑GGTTTTCCCCCCCTGTC‑AGCCCCTCCACCAGACTTTACTCCTGGTAAACCCCCCGCACCCCCAGTTCAGCGGCACGGCGCGTCAGCGTAGCTTCGTCACCTTCCAGTcgc                                           >  2:221012/1‑151 (MQ=255)
                                                                                                                 aTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGc                                    <  2:829675/151‑1 (MQ=255)
                                                                                                                   gcggTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtct                                       <  1:1124109/143‑1 (MQ=255)
                                                                                                                           ggcggGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                              >  1:1098154/1‑147 (MQ=255)
                                                                                                                             cggGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGGC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                        >  1:378660/1‑151 (MQ=255)
                                                                                                                               gggCCACTGACCTAGGCGGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGCCTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCATCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCa                      <  2:1082364/151‑1 (MQ=255)
                                                                                                                                  ccACTGACCGAGGCAGA‑GGTTTTCCGCCTCTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCAcct                   <  1:980361/151‑2 (MQ=255)
                                                                                                                                     cTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt                 <  2:1150481/150‑1 (MQ=255)
                                                                                                                                      tGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTAt               >  2:816105/1‑151 (MQ=255)
                                                                                                                                              gCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGt       <  2:760283/151‑1 (MQ=255)
                                                                                                                                                aga‑gGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGtt      <  1:420507/150‑1 (MQ=255)
                                                                                                                                                aga‑gGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc     <  1:30125/151‑1 (MQ=255)
                                                                                                                                                    ggTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  <  1:1124839/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTC‑AGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677890‑2678187

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑