BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,294	0	T	C	39.5%	201.5 / 55.1	129	intergenic (+100/+129)	yhdZ/rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand): ref base T (16/62); new base C (51/0); total (67/62)
Fisher's exact test for biased strand distribution p-value = 4.48e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

TAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCT  >  NC_000913/3423149‑3423439

tAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCGCCCCCCA‑CCCTCCcgg                                                                                                                                                  >  3:684965/1‑151 (MQ=255)

tAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCCCCCCCCC‑CCCTcc                                                                                                                                                     >  3:1057808/1‑148 (MQ=255)

 aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCggg                                                                                                                                                  <  1:248956/150‑1 (MQ=255)

  aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGcc                                                                                                                                                <  1:1030181/151‑1 (MQ=255)

   aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCg                                                                                                                                               <  3:775553/151‑1 (MQ=255)

   aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCC‑CCCCCCGGGcc                                                                                                                                                >  3:1160620/1‑150 (MQ=255)

   aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCA‑CCCCCCcgg                                                                                                                                                  >  2:573207/1‑148 (MQ=255)

      agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCA‑CCCCCCCGGCCggt                                                                                                                                             >  3:495883/1‑148 (MQ=255)

       gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGc                                                                                                                                           <  3:541154/151‑1 (MQ=255)

         gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTTCCCCCCA‑CCCCCCCGGcc                                                                                                                                                >  3:477344/1‑144 (MQ=255)

             aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑a                                                                                                                                  <  3:648962/151‑1 (MQ=255)

               gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCCGGCCGTTGCttttc                                                                                                                                      >  3:523351/1‑145 (MQ=255)

                  ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGtt                                                                                                                             <  2:226668/151‑1 (MQ=255)

                    cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCa                                                                                                                           <  1:1181142/151‑1 (MQ=255)

                      tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCGCCCCCCC‑CCCCCCCGGCCGGTGCTTT‑G‑C‑AACAtt                                                                                                                             >  3:1035097/1‑147 (MQ=25)

                      tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCAGGCCCCCACCCCCCCGGGCCGtgtctcccaaagtt                                                                                                                                 >  3:1317562/1‑134 (MQ=25)

                       ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCC‑CCCCCCCGGCCGGTGCCTT‑G‑C‑aaagattaa                                                                                                                          >  1:768896/1‑142 (MQ=255)

                       ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCCGCCCGCCC‑CCCCCCGGGCCGTTGGTTC‑T‑C‑AAAGTTCaaa                                                                                                                         >  2:104499/1‑150 (MQ=255)

                         ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCTTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGGGGGC‑GTCCTGCCCGCCT‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAAt                                                                                                                        <  1:126720/149‑1 (MQ=255)

                               cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GCCCCGTCCCCCC‑CCCCCCCCGCCGGTGGTTT‑C‑C‑AACCTTCCAAttccccccc                                                                                                                >  1:1102223/1‑143 (MQ=255)

                                aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCC‑CCCCCCGGCCCGTT‑CTTC‑T‑C‑CACGCTCAAATCCGCccccc                                                                                                               >  1:49409/1‑150 (MQ=255)

                                aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCcg                                                                                                               <  3:512369/151‑1 (MQ=255)

                                    aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCCTGCCCGCCC‑CCCCCCCGGCCGTTGCTTC‑G‑C‑AACGCTCAAATCCGCTCCccgcgg                                                                                                           >  1:177523/1‑151 (MQ=255)

                                    aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCC‑CCCCCCCGGGCGTTGCTTC‑G‑C‑AAAGTTTAAATcccccccccg                                                                                                              >  2:217709/1‑146 (MQ=21)

                                    aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGTCTGCCCCCCC‑CCCCCCCCGCCGGTGCTTC‑G‑C‑AAAATTCAAATTCGCTccccg                                                                                                              >  4:5107/1‑146 (MQ=255)

                                       ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGAt                                                                                                         <  1:1318507/150‑1 (MQ=255)

                                          ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGTTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGTCGGATTTg                                                                                                      <  1:562775/150‑1 (MQ=255)

                                          ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑CTCCTGCCCCCCG‑CCCTCCCGGCCGTTGCTCCGG‑C‑AAC‑TTCCAATACGCACCCTGCGGAatt                                                                                                       >  3:418880/1‑149 (MQ=255)

                                           tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCA‑CCCCCCCGGCCGCTGCTTT‑G‑C‑AACGTTCAAATCCGCTCCCCGCGGGTTTGt                                                                                                     >  3:263456/1‑150 (MQ=255)

                                            aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCATTCCTTGCGGATGGCCTTGCGCTTGATGTGATGTCTGGCAGTTTAGGGCGGGC‑GTCCTGCCGGCCG‑CCGTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGtg                                                                                                    <  4:167929/150‑2 (MQ=17)

                                             aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCA‑CCCCCCCGGCCGGTGCTTC‑G‑C‑AACGTTCAAATCCCCTCCccgcgggatt                                                                                                       >  3:653060/1‑142 (MQ=18)

                                              ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCCTTCCCGCCC‑CCCCCCCGGCCCTTGCTTC‑T‑C‑AAACATTAAATCCCCCCCccgcggggtt                                                                                                       >  4:393004/1‑141 (MQ=255)

                                               tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCC‑CCCCCCCGGCCGTTGCTTT‑G‑C‑CACGTTCAAAACCCCTCCccgcggggtttttcta                                                                                                 >  4:224828/1‑140 (MQ=18)

                                                 aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                              <  2:182375/151‑1 (MQ=255)

                                                  aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCA‑CCCCCCCGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCCCTCCCCGCGGGTTTGTCCttc                                                                                                >  2:565787/1‑146 (MQ=18)

                                                   tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCg                                                                                             <  3:1151971/150‑1 (MQ=255)

                                                    gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCgg                                                                                            <  3:373557/150‑1 (MQ=255)

                                                    gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGGCTGCCCCCCA‑CCCCCCCGGCCGGTGCTTC‑T‑C‑CACATTCAAAACCCCCCCCCGCGGGTTTTTGCTACTc                                                                                              >  3:289691/1‑148 (MQ=255)

                                                     ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCgg                                                                                            <  1:672821/149‑1 (MQ=255)

                                                      tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGga                                                                                          <  1:690896/150‑1 (MQ=255)

                                                      tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGga                                                                                          <  2:196694/150‑1 (MQ=255)

                                                       aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgag                                                                                         <  1:1109987/150‑1 (MQ=255)

                                                        gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCTGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGga                                                                                          <  1:298538/148‑1 (MQ=255)

                                                        gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGagag                                                                                       <  1:1402214/151‑1 (MQ=255)

                                                        gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCGCCCCCCC‑CCCCCCCGGCCGTTGCTTC‑G‑C‑AAAGTTCAAATCCCCCCCCCGCGGATTTTTCCTTCTCcgggg                                                                                         >  3:59779/1‑147 (MQ=17)

                                                          tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTCCCGCCCCCCC‑CCCCCCCGGGCGGTGCTTT‑G‑C‑AACGTTCAAATCCGCCCCCCGCGGGGTTGTTCCTCTCCGGagaggg                                                                                     >  1:128788/1‑149 (MQ=255)

                                                          tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgt                                                                                      <  4:549673/150‑1 (MQ=255)

                                                           cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCTTTTCGCTTGATTTGATGTCTGCCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCGCCCGGCGGATTTGTCCTACTCGGGAGAgtgt                                                                                    <  3:926433/151‑1 (MQ=255)

                                                           cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCGCCC‑CCCCCCCGGCCGGTGCTTC‑T‑C‑CACCCTCCAATCCCCCCCCCGCGGCTTTTTCCTCCTCCGGggaggg                                                                                     >  1:927526/1‑148 (MQ=255)

                                                            gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GATCTTCCCCCCACCCCCCCGGCCCTTTGCTT‑‑G‑C‑CAAGTTCCAATACCCCCCCGGCGGATTTTTCCttcctggggggac                                                                                      >  3:966827/1‑136 (MQ=255)

                                                              aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCC‑CCCCCCGGGCCGTTGCTTC‑G‑C‑AACGTTCCAATCCGCCCCCCGCGGGTTTGTCCTACCCCGGAGAGGGtt                                                                                   >  2:236255/1‑149 (MQ=17)

                                                                aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCCGCCCCCCA‑CCCCCCCGGCCGTTGCTTC‑G‑C‑CACGTTCAAATCCGCTCCCCGCGGGTTTGTTCTACTCGGGAGAGTGTT‑ccac                                                                              >  3:1378790/1‑148 (MQ=255)

                                                                   aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGac                                                                           <  1:647479/151‑1 (MQ=255)

                                                                    aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGaca                                                                          <  1:1306968/151‑1 (MQ=255)

                                                                      aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGGCCGCCC‑CCCTCCCGGCCGTTTCTTC‑GCC‑AAGTTTCAAAT‑CGCCCCCCGCGGATTTGTCCTACTCGGGAGAGGGGT‑CAcc                                                                              >  3:657307/1‑145 (MQ=17)

                                                                         gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGTCCCCCC‑CCCTCCCGGCCGTTGCTTC‑C‑CGAAAGTCCAAAACC‑CTCCCCCCGGACTTGTCCTACTCCGGAGGGGGTTCCACCg‑caaaca                                                                      >  1:403738/1‑150 (MQ=255)

                                                                           cATCCTTTCGGATGGCCTTTCGCTTGAGTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCGA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGacaaacaa                                                                     <  4:63055/149‑1 (MQ=255)

                                                                             tCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAaacaac                                                                    <  3:356286/148‑1 (MQ=255)

                                                                              ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGAt                                                                <  1:286357/151‑1 (MQ=255)

                                                                              ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCCGCCCCCCC‑CCCTCCCGGGCGTTGCTTT‑T‑C‑AACGCTCCAATTCGCCCCCCGCCGGTTTTTCCTCCTCCGGAGAGTGTT‑CTCCCacaaaacaaaca                                                                 >  4:58808/1‑143 (MQ=9)

                                                                               cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCGCCCGCCA‑CCCCCCCGGCCGTTGCTTC‑G‑C‑AACGTTTAAATCCGCCCCCCGCGGGGTTTTCCTACTCCGGAGAGGGTT‑TACCCacaaacaa                                                                     >  3:773420/1‑145 (MQ=14)

                                                                                tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCCCCCCCCCC‑CCCCCCCCGGCGTTGCTTC‑T‑C‑CACGTTTCAATCCCCGCCCCCCGGGGTTGTTCCACACGGGAGAGGGGT‑CACCGacaaaa                                                                       >  4:332424/1‑141 (MQ=9)

                                                                                   cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATaaaa                                                            <  4:409356/150‑1 (MQ=255)

                                                                                   cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATaaaa                                                            <  3:1248102/150‑1 (MQ=255)

                                                                                   cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATaaaa                                                            <  1:929792/150‑1 (MQ=255)

                                                                                    ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACa                                                          <  3:194241/151‑1 (MQ=255)

                                                                                    ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACa                                                          <  1:893507/151‑1 (MQ=255)

                                                                                       tGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GCCCTGCCCGCCA‑CCCTCCGGGCCCTTTCTTC‑T‑C‑AAACATCAAATACCCCCCCCGCGGGGTTGTTTTTCACTGGGGAGGGTT‑TACCGCCAAAAAACAGATaaaaaaaa                                                        >  1:811574/1‑150 (MQ=9)

                                                                                        ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAg                                                      <  1:517438/151‑1 (MQ=255)

                                                                                             ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGcc                                                   <  3:1168163/149‑1 (MQ=255)

                                                                                               cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCT‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                <  3:123147/150‑1 (MQ=255)

                                                                                               cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                <  2:34327/150‑1 (MQ=255)

                                                                                               cGCTTGATTTGATGCCTGGCATTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AATGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGcc                                                   <  1:226795/147‑1 (MQ=255)

                                                                                                gCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTCCCTCCCGCCA‑CCCCCCCGGGCCTGGCTTC‑G‑C‑AACGTTCAAATCCGCCCCCCGCGGGCTTGTCCTACTCCGGGGGGTGTT‑TACCGCCCAACaaaaaaaaaaaaaaaaagacccgcc                                              >  1:1062135/1‑142 (MQ=2)

                                                                                                 cTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCC‑CCCCCCCGGCCGTTGCTTT‑G‑C‑AACCTTCAAAACCCCTCCCCGCGGGGTTTTCCTTCTCGGGGGAGGGTT‑CACCGACAAACAACAGATAAAAAAAAAAGCCCAgg                                               >  1:584603/1‑148 (MQ=14)

                                                                                                 cTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCC‑CCCCCCCGGCCGGTGCTTT‑C‑C‑CACCTTCAAAACCCCCCCCCGCGGGTTTTTCCTTCTCCCGGGGGGGGT‑TACCCCCAAACAACAAATAaaaaaaaaaaccccacc                                              >  1:554834/1‑141 (MQ=2)

                                                                                                 cTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCC‑CCCCCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCCGCGGGGTTGTCCTACTCGGGGGAGTGTT‑TCCCGACAAACAAAAGATAAAAAAAAAAGCCCacgc                                              >  3:38149/1‑147 (MQ=14)

                                                                                                  ttaatttgatGCCTGGCAGTTTATGGCGGGC‑GTCCCTCCCGCCA‑CCCCCCGGGCCGCTGCTTT‑T‑C‑AAAGTTTCAATTCGCCCCCCGCGGGTTTATCCCACCCCGGGGAGGGAT‑TACCCGCAAAAAAAACAaaaaaaaaaaaaccccgactt                                            >  3:1062945/4‑140 (MQ=0)

                                                                                                     atttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTcc                                          <  4:7099/150‑1 (MQ=38)

                                                                                                     atttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                         <  1:1199806/151‑1 (MQ=38)

                                                                                                     atttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                         <  1:1313923/151‑1 (MQ=38)

                                                                                                     atttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                         <  4:319000/151‑1 (MQ=38)

                                                                                                                gtcAGTTTATGCCGGGC‑GCCCTGCCCGACA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCTGATTTGTTCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAACCAAAAGGCACAGCCTTCCGACTGAGccct                               <  1:1348886/148‑3 (MQ=11)

                                                                                                                ggCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                 <  2:90109/148‑1 (MQ=35)

                                                                                                                 gCAGTTTATGGCGGGG‑GTACTGCCCCCCA‑CCCCCCCGGCCGGTGCTTC‑G‑C‑AACGGTCCAATCCCCCCCCCGCGGGTTTTTTCTACTCGGGGGAGGGTT‑TACCGACAAACAACAGATAAAAAAAAAAGACCAGTCCTCCCGCTGAGCCcttt                             >  4:549358/1‑151 (MQ=2)

                                                                                                                 gCAGTTTATGGCGGGC‑GTCCCGCCCCCCA‑CCCCCCCCGCCGTTGCTTC‑G‑C‑AAAGTTCAAATCCCCCCCCGGCGGGTTTTTCCTACTCCGGGGAGGGTT‑CACCGCCAAAAAACAGATAAAAAAAAAAACCCCGTCTTCCCCCTGAGCCtttt                             >  1:73828/1‑151 (MQ=2)

                                                                                                                  cAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                            <  2:119783/151‑1 (MQ=37)

                                                                                                                   aGTTTATGGCGGGC‑GTCCTGCCCCCCC‑CCCCCCCGGCCGTTGCTTC‑G‑C‑AACGTTCAAAACCCCCCCCCGCGGGTTTGTCCTACTCCGGAGAGTGTT‑TACCCACAAACAAAAGATAAAACAAAAAGCCCCCTCTTTCGACTGCGCCcttttt                           >  3:1340754/1‑149 (MQ=2)

                                                                                                                        aTGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                       <  1:1334124/150‑1 (MQ=38)

                                                                                                                          ggcgggc‑gTCCTGCCCGCCA‑CCCTCCGGGTCGTTGCTTC‑G‑C‑AACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTCTTTTa                       <  4:176583/148‑1 (MQ=255)

                                                                                                                             gggc‑gTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                 <  2:138943/151‑1 (MQ=38)

                                                                                                                                  gTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt             >  3:1411147/1‑151 (MQ=40)

                                                                                                                                  gTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt             >  1:209035/1‑151 (MQ=40)

                                                                                                                                  gTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt             >  3:102060/1‑151 (MQ=40)

                                                                                                                                    ccTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg            <  2:336697/150‑1 (MQ=40)

                                                                                                                                    ccTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg            <  2:358221/150‑1 (MQ=40)

                                                                                                                                      tGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc          >  1:535056/1‑150 (MQ=40)

                                                                                                                                       gCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc          <  3:2975/149‑1 (MQ=40)

                                                                                                                                         ccGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGtt      >  2:504688/1‑151 (MQ=40)

                                                                                                                                             cA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt  >  1:1337780/1‑151 (MQ=40)

TAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTC‑G‑C‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTT‑CACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCT  >  NC_000913/3423149‑3423439

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