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	31.3%	680.9 / 115.4	383	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 (79/178); new base A (97/20); total (183/200)
Fisher's exact test for biased strand distribution p-value = 9.66e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 389 total aligned reads displayed.

CATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATC  >  NC_000913/1552041‑1552341

cATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                        <  4:1606060/151‑1 (MQ=31)

 atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                        >  6:503986/1‑150 (MQ=31)

       tAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGgc                                                                                                                                                 >  2:897009/1‑151 (MQ=14)

        aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcg                                                                                                                                                <  7:1010229/151‑1 (MQ=33)

        aGAAAACGTTTGGTTTTTTCATGCCGTTTGCGGGGTTAACGCCTTATCCGGCATTAAAACCCTTAAAATCAAATAGATTGCAGTTAACGTGTAGGCCTTATAAGCTTAGCGCATCAGGAAATGTTGCGTTTGTCATCAGTTTAAAATGgc                                                                                                                                                 <  4:1517521/150‑1 (MQ=255)

           cAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCtt                                                                                                                                             >  1:730319/1‑150 (MQ=14)

             aCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                                                                                                            >  3:1154333/1‑150 (MQ=35)

             aCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCACATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                                                                                                            >  6:739694/1‑150 (MQ=25)

              cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaa                                                                                                                                          <  5:38888/151‑1 (MQ=37)

               gTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCGAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCGGTaa                                                                                                                                          >  3:460588/1‑150 (MQ=25)

                  ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGAGAAGGCCTGATAAGCGTAGCGCATAAGGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTaaaat                                                                                                                                       >  8:656626/1‑149 (MQ=255)

                     tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCgg                                                                                                                                   >  5:886866/1‑150 (MQ=255)

                      ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGt                                                                                                                                   <  1:1539135/150‑1 (MQ=255)

                       tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAATGGCGTc                                                                                                                                  >  5:624309/1‑150 (MQ=255)

                       tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTc                                                                                                                                  <  6:389143/150‑1 (MQ=255)

                          tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCAttt                                                                                                                              >  1:585834/1‑151 (MQ=255)

                            aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                            <  7:1343822/151‑1 (MQ=255)

                               ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCata                                                                                                                         <  3:900882/151‑1 (MQ=255)

                                 ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGGGGCATTTTCTTAtt                                                                                                                       >  2:668017/1‑151 (MQ=255)

                                  gATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTg                                                                                                                      <  6:1150911/151‑1 (MQ=255)

                                       ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTATTCAGTTTCAAATGGCGATGTAAAAGGCGTCATTTTCATATTGTAGa                                                                                                                  >  2:1423349/1‑150 (MQ=255)

                                        gcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACa                                                                                                                <  2:250023/151‑1 (MQ=255)

                                         cgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACaa                                                                                                               <  4:1321443/151‑1 (MQ=255)

                                             aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGTAATTATAAGCGTAGCTCATCAGGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTa                                                                                                           <  1:1455971/151‑1 (MQ=18)

                                             aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGGTAAGCGTAGCGAATCATGGGATGTGGCGTGTGTCTGACCTTTCAAATGGCGCTGTAAAAGGCGTAAATTTCATATTGTAGACAACGtc                                                                                                           >  2:1502015/1‑150 (MQ=17)

                                             aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTa                                                                                                           <  5:705403/151‑1 (MQ=255)

                                             aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTa                                                                                                           <  5:1166630/151‑1 (MQ=255)

                                              aCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGt                                                                                                            >  3:1422510/1‑149 (MQ=255)

                                               cGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAATTCGTAATATACAAATTTTAGCAAACTTAgt                                                                                                         >  8:1043815/1‑150 (MQ=1)

                                                gCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCTTAGCGCATCAGGCAATGTTGAGTTTGTCATCAGTTAAAAATGGCGCTTTATAGATCGTAATATACAAATTGTAGAAAACGttgt                                                                                                         >  4:948947/1‑147 (MQ=11)

                                                       ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGTGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTc                                                                                                 <  1:120306/151‑1 (MQ=255)

                                                       ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACCAAGTAGGCTTTGTTc                                                                                                 >  2:1451101/1‑151 (MQ=255)

                                                        cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAAGAGTCATTTTAATATTGTAGAAAACGTAGACTTTGtta                                                                                                 >  1:883064/1‑149 (MQ=17)

                                                         ggCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGGGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCAt                                                                                               >  8:1281823/1‑151 (MQ=255)

                                                         ggCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAGAAATGGCGCTGTATAAGGCGTCGTATTAAGATTGTAGACAACGTATGGTTTATTCAt                                                                                               >  2:323105/1‑151 (MQ=14)

                                                              gAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTATAAAACGTAGGCTTTTTTCATGCCgg                                                                                          >  7:1278927/1‑151 (MQ=18)

                                                               aaaaCCCTTCAAATCCAATAGATTGTAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAAGAGACATTTTAATATTGAAGAAAACGTAGACATAGTACATGCCgc                                                                                          >  7:1152219/1‑149 (MQ=11)

                                                               aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCgg                                                                                          <  4:267064/150‑1 (MQ=255)

                                                                aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAt                                                                                        <  4:785412/151‑1 (MQ=255)

                                                                aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAt                                                                                        <  8:891745/151‑1 (MQ=255)

                                                                 aaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATg                                                                                       <  5:1107730/151‑1 (MQ=255)

                                                                 aaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGAAGACAACGAAGGCTTTGTTCATGCAGCATg                                                                                       >  6:567760/1‑151 (MQ=255)

                                                                    ccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTATAAAAGTCGTCATTATGATATTATAGAAAACTTATGCTTTCTTCATTCCggctaacc                                                                                    >  4:1085414/1‑145 (MQ=17)

                                                                       tCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAAGAGACATTTTAATATTGAAGAAAACCTCGGCTTTGTTCATTCCGAAAGAGGCGt                                                                                 >  8:249687/1‑151 (MQ=14)

                                                                        cAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGt                                                                                 <  2:98032/150‑1 (MQ=255)

                                                                        cAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTg                                                                                <  5:469499/151‑1 (MQ=255)

                                                                         aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGACATTTTCATATTGAATACAACGTAGGCTTTGTTCATCCCCGATCCCGCGtt                                                                                >  8:409152/1‑149 (MQ=18)

                                                                          aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGGGTCATTTTAATATTGTATACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGat                                                                              >  3:1564281/1‑150 (MQ=255)

                                                                          aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGGCATTTTATTATTGTAGACAAAGTAGGCTTAGTTCATGCCGGAAGCGGCGtta                                                                               >  4:571989/1‑148 (MQ=17)

                                                                           aTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                             <  6:225370/151‑1 (MQ=255)

                                                                           aTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                             <  3:1033033/151‑1 (MQ=255)

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

                                                                             ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAAGc                                                                           >  5:187695/1‑151 (MQ=255)

                                                                             ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAAGAGTCATTTTAATATTGTATAAAAAGTAGGCTTTGTTCATGCCGAATACCGCCTGatagc                                                                           >  8:1600895/1‑147 (MQ=17)

                                                                                  aGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                      <  4:1011897/151‑1 (MQ=255)

                                                                                  aGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                      <  6:1406405/151‑1 (MQ=255)

                                                                                  aGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                      <  5:1257426/151‑1 (MQ=255)

                                                                                  aGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAGGAGTCATTTTAATATTGTAGAAAACGTAGGCTTTGTTAAAGCAGGATGCTCCGTGAACGAATTAt                                                                      >  1:1011935/1‑151 (MQ=18)

                                                                                   gATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGGAGACAACGTAGGCTTTGTTCATACAGGATCCGGCGTGAAAGCCTTATc                                                                     >  1:422648/1‑151 (MQ=25)

                                                                                   gATTGCAGTGAACGAGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAAGCGTCATTTTCATATTGTAGACAACGTAGGCTTTTTTCATTCAGGAAGCTCCGTGAAAGACTTAt                                                                      >  7:509532/1‑150 (MQ=18)

                                                                                    aTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                                    <  2:1602519/151‑1 (MQ=255)

                                                                                    aTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                                    <  3:509781/151‑1 (MQ=255)

                                                                                       gCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAGCCGGCGTGAACGCCATAACCGGc                                                                 >  5:1002300/1‑151 (MQ=255)

                                                                                          gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATg                                                              <  3:163630/151‑1 (MQ=255)

                                                                                          gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAGGAGACATTTTAATATTGTATACAACGTAGGCTTTCTTCATGACGTATGCGACGTGAACGCCTTATCCgtaatt                                                              >  8:1444166/1‑146 (MQ=18)

                                                                                                  gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCt                                                      <  3:664156/151‑1 (MQ=255)

                                                                                                   tAGGCCTGATAAGCGTAGCGCAACAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGAAAAAGGCGTCATTTTCATATTGTAGACAACGAACGCTTTTTACATGCCGGATCCAGCATGAACCCCTAATCCTCCAACAAAACCCt                                                      >  8:1164244/1‑150 (MQ=21)

                                                                                                    aGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                    <  5:733807/151‑1 (MQ=255)

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

                                                                                                     ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGCATGAAAACCCTTc                                                    >  4:1221414/1‑150 (MQ=255)

                                                                                                         tGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTGATCCGGCATGAAAACCCTTCaaaac                                               >  4:1426168/1‑149 (MQ=255)

                                                                                                            tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAGGAGGCATTTTCATATTGTAGACAAAGTAGGCTTTGTTCATTCAGGATGCGGAGTGAACGCCTTATCCGCCATTAAAACCCTTCAAATCCa                                             >  2:763942/1‑150 (MQ=18)

                                                                                                             aatcGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAAt                                           <  7:373826/148‑1 (MQ=255)

                                                                                                              aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATa                                          <  2:1000243/151‑1 (MQ=255)

                                                                                                                cGTAGCGCATCAGTCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGa                                        <  1:56747/151‑1 (MQ=255)

                                                                                                                   aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTg                                     <  3:1270473/151‑1 (MQ=255)

                                                                                                                      gcATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTATACAACGTAGGCTTTATACATGCCCGATCCGTCTTGAAAAACATAACCGGCAAGAAAACCCTTAAAATCCAATAGATTGc                                    >  1:827607/1‑149 (MQ=17)

                                                                                                                       cATCAGTGAATGTTGTTTTTCTCTTAAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAg                                  <  7:139497/150‑1 (MQ=11)

                                                                                                                       cATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGt                                 <  3:1350314/151‑1 (MQ=255)

                                                                                                                         tCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGACGCGTGATCCCCTTATCCGGCATGAACACCCTTGAAATACAATAGATTGCAGTg                                >  3:827316/1‑150 (MQ=17)

                                                                                                                            ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                            >  6:812346/1‑151 (MQ=255)

                                                                                                                            ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                            <  6:927871/151‑1 (MQ=255)

                                                                                                                             gCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                            <  2:58766/150‑1 (MQ=255)

                                                                                                                                aTGTTGCGTTTGTCATCAGGTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTa                        <  2:46204/151‑1 (MQ=255)

                                                                                                                                 tGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgtgt                         >  2:214218/1‑149 (MQ=255)

                                                                                                                                  gtttttttttttttaaGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                       <  4:892723/135‑1 (MQ=11)

                                                                                                                                  gTTGCGTTTGTAATGAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAgg                      <  2:133470/151‑1 (MQ=255)

                                                                                                                                   ttGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATGCGGCATGAAAACCCTTCAAATACAATAGATTGCAGTGAACGTGTAgtc                     >  7:858902/1‑149 (MQ=255)

                                                                                                                                   ttGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGc                     >  8:72591/1‑151 (MQ=255)

                                                                                                                                        tttGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt                >  3:296443/1‑151 (MQ=255)

                                                                                                                                         ttgtaatgtctttCAAATTGAGCTTTAAAATTCGTTATTTTCTTTTTCTAGACTACGTTCTTTCTTTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt                <  4:693788/140‑1 (MQ=2)

                                                                                                                                         ttGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt                <  5:498591/150‑1 (MQ=255)

                                                                                                                                          ttgcatcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATa               <  4:76673/147‑1 (MQ=255)

                                                                                                                                           gttttaaGTTTAAAATGTCGCTGTAAATGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa              <  7:696353/144‑1 (MQ=255)

                                                                                                                                            ttttaaGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg             <  8:505419/145‑1 (MQ=9)

                                                                                                                                            ttttaaGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg             <  7:1544156/145‑1 (MQ=9)

                                                                                                                                            tcataaGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg             <  1:1384931/150‑1 (MQ=18)

                                                                                                                                             catcaGTTTCTCATTGCGCGTTTATTGGAGTCTTTTTCTTATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa              <  2:1199396/148‑1 (MQ=11)

                                                                                                                                             catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg           >  6:1237492/1‑151 (MQ=39)

                                                                                                                                              atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt          >  6:415905/1‑151 (MQ=38)

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