BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,552,191	0	C	A	30.0%	746.5 / 109.1	433	intergenic (+200/+207)	fdnI/yddM	formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand): ref base C (101/193); new base A (104/22); total (215/218)
Fisher's exact test for biased strand distribution p-value = 2.50e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 437 total aligned reads displayed.

CATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATC  >  NC_000913/1552041‑1552341

 atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                       <  6:1033744/151‑1 (MQ=31)

      gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAAt                                                                                                                                                    <  4:1221682/150‑1 (MQ=18)

       tAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATgg                                                                                                                                                  <  1:1079494/150‑1 (MQ=33)

        aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgc                                                                                                                                                 <  1:938557/150‑1 (MQ=33)

           aaaCGTAGTCTTTGTTCATGCAGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCTTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                             <  1:1632355/150‑1 (MQ=255)

             aCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTa                                                                                                                                           <  6:288999/151‑1 (MQ=35)

                 aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaaa                                                                                                                                         >  3:1059691/1‑149 (MQ=38)

                 aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAg                                                                                                                                       <  8:1085951/151‑1 (MQ=39)

                  ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAATCAAAGTAGAGTTTGTCATCAGTTTCAAATGGAGGTGTAAAAgg                                                                                                                                      >  1:545414/1‑151 (MQ=14)

                   gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGc                                                                                                                                     <  8:1702029/151‑1 (MQ=40)

                    gtttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                                                                                                                    <  2:527501/150‑1 (MQ=255)

                        gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCa                                                                                                                                 <  4:1690979/150‑1 (MQ=255)

                         ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGCGTCAtt                                                                                                                               >  1:1601939/1‑151 (MQ=255)

                              gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                                                                                                          <  1:1512408/151‑1 (MQ=255)

                                 ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATAtt                                                                                                                       <  3:182917/151‑1 (MQ=255)

                                 ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATAtt                                                                                                                       <  8:1155031/151‑1 (MQ=255)

                                 ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAACGAGTCATTTTGTTAtt                                                                                                                       >  8:67403/1‑151 (MQ=255)

                                  gATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTg                                                                                                                      <  4:1288546/151‑1 (MQ=255)

                                   aTGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGCGACATTTTGATATTGt                                                                                                                     >  1:623298/1‑151 (MQ=255)

                                     gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTgca                                                                                                                    <  4:533822/150‑3 (MQ=255)

                                     gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                                                                                                                     >  5:566883/1‑149 (MQ=255)

                                     gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa                                                                                                                    <  2:1476671/150‑1 (MQ=255)

                                     gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAg                                                                                                                   <  7:785543/151‑1 (MQ=255)

                                      cggcgTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGGGTCATTTTGATATTGTat                                                                                                                   >  3:1523242/1‑149 (MQ=255)

                                       ggcgTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGAc                                                                                                                 >  8:350012/1‑151 (MQ=255)

                                            gAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGt                                                                                                            <  5:462692/151‑1 (MQ=255)

                                             aaCTCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGt                                                                                                            <  1:321641/150‑1 (MQ=255)

                                             aaCGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTGATATTGTAGACAACGt                                                                                                            >  4:1114009/1‑150 (MQ=255)

                                               cGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAgg                                                                                                         <  8:1443859/151‑1 (MQ=255)

                                                 ccTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGACTCATTTTAATATTGAATAAAACGAAGGc                                                                                                        >  5:251212/1‑150 (MQ=14)

                                                  cTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCATGAAAAGCGTAGCCAATCAGGCAAGCTTGAGTTTGAGATGAGTTGCAAATGGCCATGTAAAACGCGTCATTTTCAGATTGTAGCCACCGTAGGCtt                                                                                                      >  4:999137/1‑151 (MQ=17)

                                                    tATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAACGAGGCATTTTTATATTGAATAAAACGTAggaat                                                                                                      >  4:353024/1‑146 (MQ=11)

                                                     aTCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGt                                                                                                   <  3:496990/151‑1 (MQ=255)

                                                      tCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAATGAGACATTTTAATATTGTATACAAAGAAGGCTTTAtt                                                                                                  >  7:1405399/1‑151 (MQ=14)

                                                       cccgCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGAAATGTTGCGTTTGTAATAAGTTAAAAATGGCGCTGTATAAAGAGACATTTTAAAATAAAAGAAAAGCTACGCTTTGtt                                                                                                  >  1:1225970/4‑150 (MQ=11)

                                                       ccGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGAGCTCTAAAAAGCGTCATTTTCATATTGTAGAAAACGTATGCTTTAtt                                                                                                  >  6:1309579/1‑150 (MQ=17)

                                                        cGGAATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCa                                                                                                <  7:1724/151‑1 (MQ=255)

                                                          gCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACgaacgatttcttaatc                                                                                              >  2:1666206/1‑136 (MQ=17)

                                                           cATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGc                                                                                             <  2:387472/151‑1 (MQ=255)

                                                             tGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGGGTCATTTTATTATTGTAGAAAACGTAGGATTTGTTCATga                                                                                             >  7:899968/1‑148 (MQ=17)

                                                                aaaCCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTAAAAGGCGTCATTTTGATATTGTACACAACGTAGTCTTTGTTCATGCCGCAt                                                                                        >  3:384994/1‑151 (MQ=255)

                                                                 aaCCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTCGGCTTTGTTCATGCCGGAtc                                                                                       >  7:1500333/1‑150 (MQ=255)

                                                                  aCCCTTC‑TATTCCTATAGATTGCAGTGAACGTTTAGGCCTGATTAGCGTAGCGCATCAGGCAATGTTTCGTTTGTGATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgc                                                                                      <  6:888897/151‑1 (MQ=25)

                                                                  aCCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAACGCGTCATTTTAATATTGTATACAACGTAGGCTTTGTTCATCCCCGAAgc                                                                                      >  7:399249/1‑151 (MQ=18)

                                                                   cccTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcg                                                                                     <  2:691011/151‑1 (MQ=255)

                                                                   cccTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTATAAGGCGGCATTTTCTTATTGTATACAAAGTAGGCTTTGTTCATGCAGGATgcg                                                                                     >  8:1007892/1‑151 (MQ=17)

                                                                   cccTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAAGGGACATTTTAATATTGTTTAAAAAGAAGGCTTTGTGCATGCCGGAAgcg                                                                                     >  6:979313/1‑151 (MQ=14)

                                                                    ccTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGAAGGCTTTGTTCATTCCGGATgcgg                                                                                    >  7:1268856/1‑151 (MQ=255)

                                                                    ccTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAACGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAAgcgg                                                                                    >  4:1591548/1‑151 (MQ=255)

                                                                        c‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGAGATCTAAAAAGCGTCATTTACATATTGTAGAAAACGTATGATTTGTTCATGCAGGATgcggaatg                                                                                >  6:744622/1‑147 (MQ=17)

                                                                        c‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCATATCGCATCAGGCTATTTTTCGTTTGTCATCAGTTTCAAATGGAGCTGTAAAATGCTTCATTTTAATTTTGTAGAAAACGTAGGCTTTGTTCATGCCGGATGCGGCGTg                                                                                >  6:323023/1‑151 (MQ=21)

                                                                           aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                              <  7:307154/151‑1 (MQ=255)

                                                                            aTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                             <  5:1516675/151‑1 (MQ=255)

                                                                             tCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTATACAAAGTAGGCTTTGTTCATCCAGAATGCGGCGTTaaac                                                                            >  1:1133513/1‑149 (MQ=18)

                                                                              ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                           <  3:977796/151‑1 (MQ=255)

                                                                                aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCTTATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                         <  3:1624983/151‑1 (MQ=255)

                                                                                aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                         <  4:1425450/151‑1 (MQ=255)

                                                                                 aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGGGTCATTTTGTTATTGTATACAACGTAGGCTTTGTTCATGCCGGAGGCGGCGTGATCgggtt                                                                        >  2:1269093/1‑147 (MQ=18)

                                                                                 aTAGATTGCAGTGAACGTGTAGGACTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                        <  7:1604761/151‑1 (MQ=255)

                                                                                     aTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAATTGGCGCTGTAATAGGCGTGATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCAGGATGAGGAGTGAACGCCTTATcc                                                                    >  3:1170195/1‑151 (MQ=25)

                                                                                     aTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                                    <  7:1057244/151‑1 (MQ=255)

                                                                                       tGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                   <  2:230867/150‑1 (MQ=255)

                                                                                        gCAGTGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                  <  6:606913/150‑1 (MQ=255)

                                                                                          aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGGGTTTGTCACAAGTTTCAAATGGCGCTGTAAAATGAGGCATTTTATTATTGTAGACAACGAATGCGTTGTTCATGCCGGAAGCGGCGTGAACGcgaaattcggc                                                                 >  3:872419/1‑139 (MQ=18)

                                                                                             gAACGGGTAGGCCTGATAAGCGTAGCGCATCAGTCAATGTTGCGTTTGTCATCAGGTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGa                                                             <  4:841912/150‑1 (MQ=255)

                                                                                                   gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCt                                                      <  6:185058/151‑1 (MQ=255)

                                                                                                    tAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTGATATTGTAGACAACGTAGGCTTTGTTCATCCCGGATGAGGCGTGAACGCCGTATACGGCATGAAAACCCtt                                                     >  4:802107/1‑151 (MQ=255)

                                                                                                    tAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTATACAAAGTAGGCTTTATTCATGCAGGATGCGGCGTCAACTACTTAAACGGAATGaaaaaaatt                                                     >  7:115469/1‑146 (MQ=21)

                                                                                                      ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCa                                                   <  2:1428199/151‑1 (MQ=255)

                                                                                                         cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAAt                                                >  4:447185/1‑151 (MQ=255)

                                                                                                          tGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTATACAACGTAGGCTTTGTTCATGCCGGATGCCGCGTTAAAGCATTATCCGGCATGAAAACCCTTAAAATc                                               >  4:678838/1‑151 (MQ=18)

                                                                                                           gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCATTTTAAAATGGAGCTTTAAAATGCGTCAGTTTAATATTTTAGACAACGTAGGCTTTGTTCATGAAGTATGAGGCGTTACCGAAGTATCCGGCATGAAAAACATGCAAAtac                                              >  3:785/1‑149 (MQ=16)

                                                                                                           gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAAt                                                <  8:282641/149‑1 (MQ=255)

                                                                                                           gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTCTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATc                                               >  8:1365302/1‑150 (MQ=255)

                                                                                                                 cGTAGCGCATCAGGCAATGTTGCGTTTGTCATCTGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAAt                                           >  2:15765/1‑148 (MQ=255)

                                                                                                                  gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAGGCCTTATCCGGCATGAAAACCCTTCAAATCCaaa                                           >  8:441819/1‑146 (MQ=255)

                                                                                                                  gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAGGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCATTAgct                                       >  6:1285389/1‑149 (MQ=255)

                                                                                                                    aGCGCATCAGGGAATGTTGCGTTTGTCATTAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTg                                     <  1:71461/151‑1 (MQ=255)

                                                                                                                    aGCGCATCAGGCAATGTTGCGTTTGTAATCAGTTTCAAATTGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGAAGGGTGAGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATGCAATAGATTg                                     >  5:918339/1‑151 (MQ=18)

                                                                                                                     tcgcATCAGGCAATGTTGCGTTTGTCATCAGTGTCAAATGGCGCTTTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTg                                     <  8:522933/149‑1 (MQ=255)

                                                                                                                     gcgcATGAGTCAATGTTGCTTTTGTGATCAGGTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                    <  2:965766/151‑1 (MQ=255)

                                                                                                                     gcgcATCAGGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                    <  2:1050725/151‑1 (MQ=255)

                                                                                                                     ccgcATCAGGCAATGTTGCGTGTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                    <  7:704351/150‑1 (MQ=255)

                                                                                                                        cATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCAGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGt                                 >  3:609581/1‑151 (MQ=255)

                                                                                                                          tgatgCAATGTTGTGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTg                                <  6:559105/146‑1 (MQ=255)

                                                                                                                            agtgaaTGTTGGGTTTTTCTTGTGTGTCAAATGGTGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAAc                             <  2:600270/147‑1 (MQ=11)

                                                                                                                            aGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGACATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCCGCGTGAAAGACAAAAACGGAAAGAAAACCCTTCAAAACCaaaaaaaaacaatgaa                              >  5:845202/1‑136 (MQ=17)

                                                                                                                             ggCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGACATTTTCATATTGTATACAACGTAGGCTTTGTACATGCCGGATGCGGCGTGAACGCCGTAAACGGCAAGAAAACCCTTCAAATCCATTAGCTAGCAGTGAACg                            >  5:1342454/1‑151 (MQ=14)

                                                                                                                               taatgttttttttGTCTTCAGTTTAAAATGGCGCTGTAAAAGGCTTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgt                           <  5:126209/140‑1 (MQ=11)

                                                                                                                               cAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATCCCGGATGCGGCGTGAACGCCAGATGCGGCATGAAAACCCGTAAAATCCAATAGAATGCAGTTCACgtg                          >  6:1533642/1‑151 (MQ=17)

                                                                                                                               cAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTTAAAGACTTATCCGGAATGAAAACCCTTCAAAACCAATAGATTGCAGTGAACgtg                          >  1:435862/1‑151 (MQ=17)

                                                                                                                                aaTGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgtgt                         >  4:1456475/1‑151 (MQ=255)

                                                                                                                                 aTGTTGCGTTTGGCATCAGGTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTa                        <  2:838500/151‑1 (MQ=255)

                                                                                                                                        gTTTGACATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTg                  >  8:888209/1‑150 (MQ=255)

                                                                                                                                              tattactttCTAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc            <  7:611611/144‑1 (MQ=16)

                                                                                                                                              catcaGGTTGTAATTGCTCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc            <  8:221297/150‑1 (MQ=255)

                                                                                                                                              catcaGCTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa              <  2:1081356/148‑1 (MQ=37)

                                                                                                                                               atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt          >  5:647768/1‑151 (MQ=38)

                                                                                                                                                 caGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa         <  2:125434/150‑1 (MQ=37)

                                                                                                                                                  aGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc       >  3:1326140/1‑151 (MQ=35)

CATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTC‑AAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATC  >  NC_000913/1552041‑1552341

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