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	29.2%	739.3 / 113.0	431	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 (98/195); new base A (92/29); total (203/228)
Fisher's exact test for biased strand distribution p-value = 1.54e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 437 total aligned reads displayed.

CATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATC  >  NC_000913/1552041‑1552341

   atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaaa                                                                                                                                                     <  8:481346/151‑1 (MQ=33)

   atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTACGCATGATAAGCGTAGCGCATAACGCAATGTTGCGTTTGTCATCAGTTTaaa                                                                                                                                                      >  1:703589/1‑150 (MQ=255)

    ttGTAGACAACGTATGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaaa                                                                                                                                                     <  5:1152944/150‑1 (MQ=18)

     tGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAAt                                                                                                                                                    <  3:1615147/150‑1 (MQ=33)

      gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATgg                                                                                                                                                  <  3:1422687/151‑1 (MQ=33)

            aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGATGt                                                                                                                                            >  5:261769/1‑151 (MQ=21)

              cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAACCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGtt                                                                                                                                           >  3:1371043/1‑149 (MQ=25)

              cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTTGCATCAGTTTCAAATGTCGCTGTa                                                                                                                                           >  4:916338/1‑150 (MQ=255)

              cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaa                                                                                                                                          <  7:45156/151‑1 (MQ=37)

                  ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTTTCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGTTGCGTTTGTCATGAGTTTCATATGGCGCGGGAGAAgg                                                                                                                                      >  3:513868/1‑151 (MQ=11)

                  ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTaaaa                                                                                                                                        >  4:644926/1‑149 (MQ=14)

                  ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTGGAAATGGCGCTGTATAAgg                                                                                                                                      >  7:94208/1‑151 (MQ=9)

                   gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGc                                                                                                                                     <  4:1497961/151‑1 (MQ=40)

                    cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAgg                                                                                                                                      >  2:771329/1‑149 (MQ=18)

                     tttGTTCATGCCGGATGCGGCGTGAGCGCCTTATCGGGCATGAAACCCCTTCAAATCGAATAGATTGCCGTGAACGTGTCTGCCTGATACGCGTCGCGCATCAGGCAATGTTGCGGTTGTCGTCAGTTTCAAATGGCGCTCTAAAAGGGGt                                                                                                                                   >  7:1202302/1‑151 (MQ=11)

                     tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGGGGTTGTCGTCAGTTTCTAATGGCGCGGTAAAAGGCGt                                                                                                                                   >  3:197259/1‑151 (MQ=255)

                      ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                                                                                                                    <  8:532323/149‑1 (MQ=255)

                       tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTATAAgggggc                                                                                                                                  >  8:1554924/1‑146 (MQ=16)

                       tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTc                                                                                                                                  <  4:311472/150‑1 (MQ=255)

                        gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCa                                                                                                                                 <  7:1531603/150‑1 (MQ=255)

                           gaTGCCGGATGCGGCGTTAACGCGTTATCCGGAATGAAAACCCTTGAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAttt                                                                                                                              <  6:817940/149‑1 (MQ=255)

                           cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAttt                                                                                                                              <  7:990514/150‑1 (MQ=255)

                           cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAAGGTTGCGTTTGTCATCAGGTTCAAATGGCGCTGTAAAAGGCGTCAgttt                                                                                                                             >  3:1182311/1‑151 (MQ=255)

                           cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCGACAAATCCAATAGATTGCAGAGAACGTGTAGGCCTGATAACGGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGGTTCAAATGGCGCTGTAAAAGGCGTCAgttt                                                                                                                             >  7:686636/1‑151 (MQ=14)

                            aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTAAAAGGGGTCAttttt                                                                                                                            >  5:1256389/1‑150 (MQ=255)

                            aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCAttttt                                                                                                                            >  8:330038/1‑150 (MQ=255)

                             tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCGGTAAAAGGCGGCATTTTc                                                                                                                            >  5:712785/1‑150 (MQ=255)

                              gCCGTATGCGTCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                                                                                                          <  3:76731/151‑1 (MQ=255)

                              gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGAGAACGTGTACGCCTGATAAGCGTAGCGAATCAGTCAATGTTGCGTTTGTCATCAGTTTCAAATGGAGAAATAAAAGGCGTCATTTTCat                                                                                                                          >  4:518290/1‑151 (MQ=14)

                               ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTAAAAGGCGTCATTTTGat                                                                                                                          >  7:1302056/1‑150 (MQ=255)

                                cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAAGAGACATTTTGatat                                                                                                                        >  4:19217/1‑151 (MQ=11)

                                cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATAGCAGTGAACGTGTAGGAATGATAAGCGTAGCGCAAAAGGAAATGTTGAGTGTGTCATCAGTTTCAAATAGAGCTGTAAAAGGCGTCAGTTTCatat                                                                                                                        >  2:49427/1‑151 (MQ=17)

                                   aTGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAACGAGTCATTTTAATATTGt                                                                                                                     >  1:1609234/1‑151 (MQ=255)

                                     gcggcgGGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAg                                                                                                                   <  3:1563424/151‑1 (MQ=255)

                                       ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTAGCGTATGACATAAGTTTAAAATATCGCTGTAAAAGGAGTCATTTTAATATTGTATAc                                                                                                                 >  6:1244169/1‑151 (MQ=14)

                                          gTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTAGAAAAc                                                                                                              >  1:1416688/1‑151 (MQ=255)

                                           tGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTATACaa                                                                                                               >  2:837550/1‑149 (MQ=255)

                                                   ttATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAACGGGACATTTTGTTATTGAATACAACGAAGGCttt                                                                                                     >  8:1294573/1‑151 (MQ=11)

                                                   ttATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTATAATGCGACATTTTGATATTGTATACAAAGAAGGCgtt                                                                                                     >  3:1502415/1‑151 (MQ=11)

                                                    tATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGAAGACAAAATAGCCttat                                                                                                    >  2:558483/1‑149 (MQ=14)

                                                     aTCCGTCATGAAAACCCTTCAAATCCAATAGATTGCACTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTg                                                                                                    <  6:244580/150‑1 (MQ=255)

                                                     aTCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGATATGGCGCTGTAACGGGCGACATAGTGAAAATCTAGACAACCTAGTATTTGt                                                                                                   >  8:704645/1‑151 (MQ=11)

                                                        cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAATGTGTAGGCCTGATAAGATTAGCGAATAAGTCAATATTGCGTTTGTAATAAGTTTAAAATGGCGCTGTAAAAGGCTTCATTTTAATATTGTAGACAACGTAGGCTTTGTTc                                                                                                 >  2:1365821/1‑150 (MQ=18)

                                                          gCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGGGGCATTTTAATATTGAATACAAAGTAGTCTTTGTTCAtt                                                                                              >  2:377248/1‑150 (MQ=14)

                                                            attaaaaCCCTTCAAATCCTATAGATTGCATTGAACGTGTAGGCCTGATAAGCTTATGTCATCAGTCAATGTTTATTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGcc                                                                                            <  1:502881/148‑1 (MQ=21)

                                                                  aCCCTTTAAAACCAATAGATTGAACTGCACGTGTAGTCCTGATAAGCGTATCGCATCAGTCAATGTTGCGTTTGTCTTAAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgc                                                                                      <  3:635485/151‑1 (MQ=25)

                                                                   cccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGGAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcg                                                                                     <  7:544714/151‑1 (MQ=255)

                                                                     cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAAGAGTCATTTTCATATTGTATAAAACGTAGGCTTAATTCAtacacaaaacgtc                                                                                   >  2:923723/1‑139 (MQ=17)

                                                                       tCAATTCCAATAGATTGCATTGAAGGTGTAGGCCTGATAAGCGTAGCGCATCATGCAATGTTGTGTTTGTCATAAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggcg                                                                                  <  8:915332/150‑1 (MQ=25)

                                                                       tCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggcg                                                                                  <  4:1587002/150‑1 (MQ=255)

                                                                       tCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggc                                                                                   <  7:24072/149‑1 (MQ=255)

                                                                         aaaaCCAATATATTTAATTGAAAGTGGAGGCCTGATAAGCGTAGCTCATCAGGCAATGTTGTGTTTGTCATCCGTTTTAATTGGCTCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGa                                                                               <  1:1337300/151‑1 (MQ=21)

                                                                           aTCCAATAGATTGCAGTTAACGTGTAGGCCTGATAAGTGTAGCGCATGAGGCAATGTTGAGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                             <  1:1042245/151‑1 (MQ=255)

                                                                            tCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGAAGGCTTTGTTCATGCCGGATCCGGCGTTAAcc                                                                            >  1:1127724/1‑150 (MQ=255)

                                                                             ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGGGTGAACGc                                                                           >  5:1633624/1‑151 (MQ=255)

                                                                             ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAAAGTAGGCTTTTTTCATCACGCATGCGGCGTGAACgg                                                                           >  3:1287588/1‑150 (MQ=18)

                                                                              cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTATAAGGGGACATTTTAATATTGTATACAACATATGCTTAGTTCATGCAGGATgggaagagaaagcc                                                                          >  2:195010/1‑138 (MQ=14)

                                                                               aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAACAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                         <  8:278659/151‑1 (MQ=255)

                                                                                 tAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGCGTCATTTTGATATTGTAGAAAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACtt                                                                        >  7:1003212/1‑150 (MQ=255)

                                                                                  aGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                       <  3:471945/150‑1 (MQ=255)

                                                                                  aGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGGCATTTTGATATTGTAGACAACGTAGGCTTTGTTCATGCCCGATGCGGCGTGAACGGCGTAt                                                                      >  5:1557541/1‑151 (MQ=21)

                                                                                    aTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                                    <  5:431021/151‑1 (MQ=255)

                                                                                      tGCAGTGAACGTGTAGGACTGATAAGCGTAGCGCATCAGTCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                   <  2:631818/150‑1 (MQ=255)

                                                                                       gCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                 <  8:1166627/151‑1 (MQ=255)

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

                                                                                          gTGAACGTGTAGGACTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATg                                                              <  1:1657510/151‑1 (MQ=255)

                                                                                            gAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGa                                                             <  2:1205127/150‑1 (MQ=255)

                                                                                            gAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaa                                                            <  6:1691117/151‑1 (MQ=255)

                                                                                            gAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGa                                                             <  5:570074/150‑1 (MQ=255)

                                                                                             aaCGTGTAGGCCTGATAAGCTTATCGCATCAGGCAATTTTGCTTTTGTCATCAGTTTCAAATTGCGCTGTAATAGGCGTCATTTTGATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaa                                                           <  2:1588462/151‑1 (MQ=21)

                                                                                               cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTATAAAACATCATATTAAAATAGTATACAACAAAAGATTTGATCATTCAGTATGAGTCGTGAAAGCCGTATACAGAATGaaaa                                                          >  2:1057881/1‑150 (MQ=2)

                                                                                                    aGGCCTTATAAGCGTATCTCATCTGGCAATGTTGTTTTTGTCATCTTTTTCATATGGAGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCtt                                                     <  4:1067810/150‑1 (MQ=17)

                                                                                                          gataattttAGCGCTTCATGCATTTTTGCTTTTTGAATAAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTATGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAAt                                                <  1:926072/141‑1 (MQ=14)

                                                                                                          gATAACCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATcc                                              <  5:1027275/151‑1 (MQ=255)

                                                                                                           aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCa                                             <  6:148277/151‑1 (MQ=255)

                                                                                                              aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATa                                          <  3:392374/151‑1 (MQ=255)

                                                                                                              aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGGGTCATTTTGATATTGTAGACAACGTAGGCTTTGTTCATTCCGGATGAGTCGTGAAAGCCTTTTCCGGCATGACACGCCTTCAACTACAAt                                           >  7:1174469/1‑150 (MQ=18)

                                                                                                                cGTAGCGCATGAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAg                                         <  1:621197/150‑1 (MQ=255)

                                                                                                                cGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCAGGATGCGTCGTGAAAGCCTTATGCGGCATGAAAACCCTTCAAATCCAAt                                           >  1:1196505/1‑148 (MQ=18)

                                                                                                                    gcccATCAGGAAATGTAGCGTTTGTCATCAGTATCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                    <  6:1340451/148‑1 (MQ=255)

                                                                                                                     cgcATCAGGCAATGTTGCGTTTGTCATCATTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                   <  1:27698/151‑1 (MQ=255)

                                                                                                                      gcATCAGGTAATGTTTCTTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAg                                  <  6:596544/151‑1 (MQ=255)

                                                                                                                       gaTCATGCAATGTTGCGTTTGTCATCAGTTTCACATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                   <  8:1208961/148‑1 (MQ=255)

                                                                                                                          gaGGCAATGTTGAGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGaa                              <  4:340212/150‑1 (MQ=255)

                                                                                                                             gCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAAGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgt                           >  1:1255746/1‑151 (MQ=255)

                                                                                                                             gCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCAGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATACATTGCAGTGAACgg                           >  5:470481/1‑150 (MQ=255)

                                                                                                                               aaTGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTGATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGGCTTATACGGCATGAAAACCCTTAAAATACAATAGATTGCAGTTAACgtgt                         >  8:165395/1‑151 (MQ=17)

                                                                                                                               aaTGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGGGTCATTTTGTTATTGTGTACAAAGTAGGCTTTGTTCATGCCGGATGCGTCGTGAACGAGTGATGCGGAATGACAACAGTTCAAATGCAATAGAGTGACGTGAACggtt                         >  7:1253360/1‑148 (MQ=9)

                                                                                                                                  gtttcttttGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAgg                      <  4:1264366/145‑1 (MQ=255)

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

                                                                                                                                     tttttttttataaGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                    <  6:398847/138‑1 (MQ=255)

                                                                                                                                     gcttttGTCATCAGTTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCt                   <  5:1448377/148‑1 (MQ=255)

                                                                                                                                     gCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCt                   >  4:1378959/1‑151 (MQ=255)

                                                                                                                                          tGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATa               >  2:278014/1‑150 (MQ=255)

                                                                                                                                            tcatcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc            >  4:1193677/1‑151 (MQ=39)

                                                                                                                                             tatcaGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc            <  7:1584705/149‑1 (MQ=25)

                                                                                                                                              atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTTTCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt          >  8:1590104/1‑151 (MQ=38)

                                                                                                                                               tcaGTTTCAATGTGCGCTTTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc            <  5:776018/148‑1 (MQ=255)

                                                                                                                                                 aGTTTCAAATGTCGATGTAAAAGGCGTCATTTTCATATTGTAGACAACTTAGGCTTTTTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGGTGTGATCGTGTAGGCCTGATAAGCGTAgg       <  4:712197/151‑2 (MQ=11)

                                                                                                                                                 aGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc       >  6:755897/1‑151 (MQ=35)

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