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	25.4%	924.9 / 86.9	456	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 (117/215); new base A (95/18); total (222/234)
Fisher's exact test for biased strand distribution p-value = 4.08e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 460 total aligned reads displayed.

CATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATC  >  NC_000913/1552041‑1552341

 atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                       <  3:96547/151‑1 (MQ=31)

    ttGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAAt                                                                                                                                                    >  2:1402225/1‑151 (MQ=14)

    ttGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCATGATAAGCTTAACGCATCAGGAAATGTTGCGTTTGTCATAAGTTTAAAAt                                                                                                                                                    >  6:1196474/1‑151 (MQ=11)

     tGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTGATCAGTTTCAAAt                                                                                                                                                    <  5:245760/150‑1 (MQ=18)

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

         caCAACGTAGGCTTTCTTCATGCCGGATGCGGCGTGAACGCCTTATCCGTCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcgc                                                                                                                                               <  6:211629/150‑1 (MQ=18)

           cAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                             <  1:100573/151‑1 (MQ=34)

           cAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                             <  5:872929/151‑1 (MQ=34)

              cctAGGCTTTGGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaa                                                                                                                                          <  6:1063089/149‑1 (MQ=37)

              cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCAGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAGAAGGCAATGGTGCGTTTGTCATCAGTTTACAATGGCGGTGTa                                                                                                                                           >  4:1043238/1‑150 (MQ=11)

               gTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaaa                                                                                                                                         <  1:1430801/151‑1 (MQ=37)

                tAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaact                                                                                                                                        >  2:1097774/1‑149 (MQ=37)

                 aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaaaa                                                                                                                                        <  7:1106529/150‑1 (MQ=38)

                    cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGACACCCCTTCAAATCCAATAGATTGCTGTGCCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGTTGCGTTTGTCCTCAGTTGGAAATGGCGCGGTAAAAgggg                                                                                                                                    >  1:1378698/1‑149 (MQ=11)

                     tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGAATAAGGAAAAGTTTCGTTTTTCATCAGTTTCAAATGGATATGTaaaatgagt                                                                                                                                   >  8:1010166/1‑146 (MQ=11)

                     tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCATTCAAATCCAATAGATTGCAGTGAACGTGTAGCCCTGATAAGCGAAGCGCATCAGGCAATGTTGCGTTTTTCATCAGTTAAAAATGGACCTGTAAAAAGcc                                                                                                                                    >  8:1142510/1‑149 (MQ=11)

                      ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTc                                                                                                                                  >  5:730575/1‑151 (MQ=255)

                      ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCAGGCATGAAAACACTTCAAATCCAATAGATTGCAGTGAACGTGTATGCCTGATAAGCGTAGCGCATCAGGCAATCTTGCGTTTGTCATCAGTTTCAAATGTCGCTGTAAAAGGCGtg                                                                                                                                  >  8:1374949/1‑150 (MQ=255)

                       tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTc                                                                                                                                  >  1:1168571/1‑150 (MQ=255)

                          tCATGCCGGATGCGGCGTGAACGCCTTATCCGTCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGCCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtt                                                                                                                               <  3:359025/150‑1 (MQ=255)

                          tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtt                                                                                                                               <  1:875672/150‑1 (MQ=255)

                          tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGAATAGCGAATCAGTCAATAATACGTTTGTAATAAGTTTAAAATGGCGCTGTAAAAGTAGTAAttt                                                                                                                              >  6:11422/1‑151 (MQ=11)

                           cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtttt                                                                                                                             <  1:1259756/151‑1 (MQ=255)

                                cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                                                                                                        <  7:792358/151‑1 (MQ=255)

                                cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                                                                                                        <  2:488692/151‑1 (MQ=255)

                                 ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTCATAtt                                                                                                                       >  4:434870/1‑151 (MQ=255)

                                 ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGGCATTTTGTTAtt                                                                                                                       >  7:405437/1‑151 (MQ=255)

                                 ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTATAAGGAGTCATTTTAATAtt                                                                                                                       >  8:54853/1‑151 (MQ=255)

                                    tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGtt                                                                                                                    >  4:588510/1‑150 (MQ=255)

                                       ggcgTGAACGCCTTATCCCGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGAc                                                                                                                 <  2:701112/151‑1 (MQ=255)

                                          gTGAACGCGTTATCCGGCATGAAAACCCTTCAAATCCTTTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAc                                                                                                              <  7:13879/151‑1 (MQ=255)

                                          gTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTataaaa                                                                                                               >  7:1542327/1‑145 (MQ=255)

                                          gTGAACGCCTTATACGTAATGAAAACCCTTCAAATCCAATAGATTGCAGTGAAAGTGTAGTTCTGATAATCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAc                                                                                                              <  8:73761/151‑1 (MQ=18)

                                             aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTa                                                                                                           <  4:821831/151‑1 (MQ=255)

                                             aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTa                                                                                                           <  3:1068357/151‑1 (MQ=255)

                                                 ccTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCt                                                                                                       <  2:1014220/151‑1 (MQ=255)

                                                     aTCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTGATATTGTATAACAAGTAGGCTTAGt                                                                                                   >  6:1537091/1‑151 (MQ=17)

                                                     aTCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGt                                                                                                   <  2:1133969/151‑1 (MQ=255)

                                                           cATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAACGAGACATTTTCATATTGTATACAAAGTAGGCTTAGTTCATGc                                                                                             >  7:1016940/1‑151 (MQ=17)

                                                            aTGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTAAAAGGGGTCATTTTTATATTGTGTACCAAGTAGGCTTTGTGCATGcc                                                                                            >  3:83657/1‑151 (MQ=17)

                                                               aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTATAAGGAGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATTCCgt                                                                                          >  3:1550931/1‑149 (MQ=17)

                                                                aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGTACTGATAAGCGTAGCTCATCAGGCAATGTTTCGTTTCTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAt                                                                                        <  8:154550/151‑1 (MQ=255)

                                                                aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAt                                                                                        >  5:572199/1‑151 (MQ=255)

                                                                 aaCGCTTCAAATCCAATAGATTGCAGTGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATg                                                                                       <  7:1216346/151‑1 (MQ=255)

                                                                 aaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCATATCTCAACAGATAAAAATTGCGCTCAAATAGGTGTCATTTTTATATTGTAGACAAAGTAGGATTTTTTAATgaatgatg                                                                                       >  2:369876/1‑144 (MQ=9)

                                                                  aCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgc                                                                                      <  1:334679/151‑1 (MQ=255)

                                                                  aCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGGTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCAGGATgc                                                                                      >  5:1576638/1‑151 (MQ=255)

                                                                    ccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTATAAAAAGAAGGCTTTCTTCAtacagcatcagt                                                                                    >  1:1658059/1‑140 (MQ=17)

                                                                     cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGAGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTAAAAGGGGTCATTTTGATATTGTATACAAAGAAGGCTTTGTTCATGCCGGATgcggc                                                                                   >  2:1218389/1‑151 (MQ=18)

                                                                       tcttaTCCAATAGATTGAAGTGAACGTGTAGGCCTGATAAGCGTAGCCCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggc                                                                                   <  6:1615151/145‑1 (MQ=255)

                                                                       tCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggc                                                                                   <  2:215119/149‑1 (MQ=255)

                                                                        cAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTg                                                                                <  3:1486760/151‑1 (MQ=255)

                                                                           aTCCCATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                             <  7:246627/151‑1 (MQ=255)

                                                                           aTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                             >  5:628517/1‑151 (MQ=255)

                                                                             ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                            <  3:125359/150‑1 (MQ=255)

                                                                             ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGACCTAAAAAAGGCGTCATTTTCATATTGTAGACAAAGTAGGAATTGATCATGCCGGAAgcgaagagaaca                                                                            >  6:845156/1‑141 (MQ=17)

                                                                              cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                          <  5:963562/151‑1 (MQ=255)

                                                                                  aGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                      >  3:1651503/1‑151 (MQ=255)

                                                                                   gATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                     >  3:523274/1‑151 (MQ=255)

                                                                                   gATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGACATTTTCATATTGTATAAAACGTACGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATc                                                                     >  1:72719/1‑151 (MQ=25)

                                                                                   gATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGAATCAGTCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTATGCTTTGTTCATGCCGGATGCGGCGTGAACTCCTTATc                                                                     <  3:1372176/151‑1 (MQ=255)

                                                                                        cAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTGATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCATAGCCGGc                                                                 >  1:10814/1‑150 (MQ=255)

                                                                                        cAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGACGGATGCGGCGTGAACGCCTTATCCGGc                                                                 >  1:1513100/1‑150 (MQ=255)

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

                                                                                            gAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGTCGTAATTTTCATATTGTAGACAACATAGTCTTTGTTAATGCAGGATGCGGCGTGAACGCGTTATcatacatga                                                             >  4:1425116/1‑142 (MQ=21)

                                                                                              aCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaa                                                           >  8:1644951/1‑150 (MQ=255)

                                                                                                gtttAGGGCTGATAATCGTAGCGCATGATTCAATGTTTAGTTTTTTATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATGTTGTAGACAAGGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATGCGGCATGAAAAcc                                                        <  4:1100169/148‑1 (MQ=17)

                                                                                                  gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAccc                                                       <  8:263331/150‑1 (MQ=255)

                                                                                                  gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCCCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAccc                                                       <  2:1325100/150‑1 (MQ=255)

                                                                                                       ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaaa                                                 <  4:1264759/151‑1 (MQ=255)

                                                                                                          gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAAt                                                <  6:45512/149‑1 (MQ=255)

                                                                                                              aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGAGGCGTGATCGGCGGATCCGGCATGAAAACCCTTCAAATCCAAt                                           >  3:1394506/1‑150 (MQ=18)

                                                                                                               gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCGTATCCGGCATGAGAGCCCTTCAAATCCAAt                                           >  2:268838/1‑149 (MQ=255)

                                                                                                                cGTAGCGCATCAGGCAATGTTGTGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATCCGGCGTGAACGACTAATCCGGCATGAAAACCCTTCAAATCCAAt                                           >  4:685197/1‑148 (MQ=255)

                                                                                                                 gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTCATCCGGCATGAAAACCCTTCAAATCCAATAGAt                                       <  1:579458/151‑1 (MQ=255)

                                                                                                                  tAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGAtt                                      >  3:46214/1‑151 (MQ=255)

                                                                                                                     cgcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTATAAGGGGTCATTTACATATTGTAGACAACGTAGGCTTTGTTCATCACGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAACATTAAAATCCAATAGATTgg                                    >  3:608906/1‑149 (MQ=18)

                                                                                                                      gcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAg                                  <  1:63238/151‑1 (MQ=255)

                                                                                                                       catcatGCAATGTTGAGTTTGTCATCATTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGt                                 <  6:702632/151‑1 (MQ=255)

                                                                                                                       aaTCATGCAATGTGGCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGt                                 <  1:132405/150‑1 (MQ=255)

                                                                                                                        aTCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAAGACTTTTCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTg                                >  5:1517992/1‑151 (MQ=255)

                                                                                                                        aTCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGAAGACAACGTAGCCTTTGTTCATGCCGGATGCGGCGTCAACGCCTTATCAGGCATGAAAACCCTTCAAATCCAATACAtttaaat                                 >  6:1223047/1‑145 (MQ=17)

                                                                                                                         tCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTTAAAGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTg                                >  3:901993/1‑150 (MQ=255)

                                                                                                                          gatttaatgtttatttttTCATCAGTTTAAAATGGCGCTTTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGaa                              <  4:904436/137‑1 (MQ=11)

                                                                                                                           aGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAATAGGCGGCATTTTGATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAAc                             <  4:26680/151‑1 (MQ=255)

                                                                                                                            ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                            <  3:702811/151‑1 (MQ=255)

                                                                                                                             ggaaTGTTGCTTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                            <  3:259926/148‑1 (MQ=255)

                                                                                                                               aaTGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgtgt                         <  5:275240/151‑1 (MQ=255)

                                                                                                                               aaTGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTCTTCATGCCCGATGCGGCTTGAAAGACGAATCCGGCATGAACACCCTTCAAATCCAATAGATTGCAGTGAACgtct                         >  7:649032/1‑149 (MQ=14)

                                                                                                                                  gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGAt                                       >  6:582048/1‑134 (MQ=255)

                                                                                                                                  gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAgg                      <  4:160394/151‑1 (MQ=255)

                                                                                                                                      cGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTg                  >  7:899316/1‑151 (MQ=255)

                                                                                                                                      cGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTg                  >  5:1219647/1‑151 (MQ=255)

                                                                                                                                          tttcatcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTATGCCTGATa               <  5:446223/148‑1 (MQ=255)

                                                                                                                                           gTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCAAATAGATTGCAGTGAACGTGTAGGCCTGATa               >  1:453808/1‑149 (MQ=255)

                                                                                                                                           cgcatcaGGTGAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa              <  2:151575/148‑1 (MQ=2)

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

                                                                                                                                                  gTTTCAAATGGCGCTGTAACAGCCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg      <  3:1008588/151‑1 (MQ=25)

                                                                                                                                                  gTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg      <  5:1022934/151‑1 (MQ=18)

                                                                                                                                                      cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATc  >  6:1358764/1‑151 (MQ=32)

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