BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,828,980	A→C	8.0%	Q12H (CAA→CAC)	xanP →	xanthine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,828,980	0	A	C	8.0%	1607.2 / 22.6	626	Q12H (CAA→CAC)	xanP	xanthine permease
Reads supporting (aligned to +/- strand): ref base A (314/240); new base C (8/61); total (323/303)
Fisher's exact test for biased strand distribution p-value = 2.72e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 650 total aligned reads displayed.

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828831‑3829129

 tataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAc                                                                                                                                                       <  1:3445534/150‑1 (MQ=255)

 tataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAc                                                                                                                                                       <  2:5561039/150‑1 (MQ=255)

 tataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAc                                                                                                                                                       <  1:967309/150‑1 (MQ=255)

 tataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAc                                                                                                                                                       <  1:2168940/150‑1 (MQ=255)

  ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCAcc                                                                                                                                                      >  2:4280444/1‑150 (MQ=255)

  ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                      >  2:5205899/1‑150 (MQ=255)

   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCCCAcccc                                                                                                                                                     >  1:1883095/1‑149 (MQ=255)

    aaaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCgg                                                                                                                                                    >  2:4509900/1‑150 (MQ=255)

     aaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCG‑t                                                                                                                                                  >  1:2388900/1‑150 (MQ=255)

      aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑tt                                                                                                                                                 >  2:3273954/1‑150 (MQ=255)

           tGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCa                                                                                                                                            >  1:5494090/1‑150 (MQ=255)

           tGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCa                                                                                                                                            <  2:4934359/150‑1 (MQ=255)

             cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGa                                                                                                                                          <  2:2019425/150‑1 (MQ=255)

             cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGa                                                                                                                                          <  2:5751826/150‑1 (MQ=255)

              aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTATCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGAc                                                                                                                                         >  1:5586501/1‑150 (MQ=255)

                ttATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGATTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGAATGTTTTCGATGTATGTTACCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTc                                                                                                                                       >  2:2469947/1‑150 (MQ=255)

                   tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCaaa                                                                                                                                    >  1:4811050/1‑150 (MQ=255)

                    ctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGTTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCaaaa                                                                                                                                   <  2:1232039/150‑1 (MQ=255)

                    ctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCaaaa                                                                                                                                   <  2:1416824/150‑1 (MQ=255)

                     tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTATTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAAc                                                                                                                                  >  1:1487514/1‑150 (MQ=255)

                      cATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACa                                                                                                                                 <  2:5372426/150‑1 (MQ=255)

                        tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCaaaaaacc                                                                                                                               >  1:4912007/1‑146 (MQ=255)

                        tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAAAAGc                                                                                                                               >  1:2053842/1‑150 (MQ=255)

                           aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGaa                                                                                                                            >  2:5287676/1‑150 (MQ=255)

                           aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTGGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCTAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGaa                                                                                                                            <  2:794023/150‑1 (MQ=255)

                             atatTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACACTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAAAAGCGAACt                                                                                                                          >  2:546555/1‑150 (MQ=255)

                                     aTGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTAcc                                                                                                                  <  2:5025883/150‑1 (MQ=255)

                                       gAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGt                                                                                                                >  2:972446/1‑150 (MQ=255)

                                       gAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCG‑GTGCGCAGACACAAAACAACGCACTCATTTAACGt                                                                                                                >  1:4660128/1‑150 (MQ=255)

                                        tggttgtgTCAATCGTTTGATTTTCTCTTTCGGTCAGATAAAAGGCCCGCTTTGCTACCTGGGTTTTTTTTTCGATGTCTGTTTCCACCCTCGTGCCATAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTc                                                                                                               <  1:4331522/146‑1 (MQ=255)

                                        aGGATGTGGCAAACGTTTGCTGTTCTCTTTCGGTCAGCTAAAATGCCCGCGTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCACGAGACAGAAAATGCGCCCCCGG‑TTGCGCAGACTCAAAACAGCGCCCTGATTTACCCTc                                                                                                               >  2:2629364/1‑150 (MQ=255)

                                          gATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCCCCCG‑TTGCGCAGACACAAAACAGCGAACAGAGTTAcccgctt                                                                                                             >  1:3820492/1‑145 (MQ=255)

                                          gATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCaaaaa                                                                                                                                  >  1:3736956/1‑128 (MQ=255)

                                          gATGTGGCAAACGTATGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCtt                                                                                                             >  2:3846729/1‑150 (MQ=255)

                                           aTGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTg                                                                                                            >  2:6385239/1‑150 (MQ=255)

                                                gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCTCCCTCGAGTCAGAAAATGCGGAACCGG‑TTGCGCAGACTCAAAACAGCGAAGTGATTTACCGTCTTTAAGAt                                                                                                       >  2:3191621/1‑150 (MQ=255)

                                                gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGAt                                                                                                       >  2:5110385/1‑150 (MQ=255)

                                                gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCCACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGAt                                                                                                       >  2:3057233/1‑150 (MQ=255)

                                                 cAAACGTGTGCTGTTCTCTTGCGGTCAGCTAACAGGCCCGCGTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                      <  2:6578403/150‑1 (MQ=255)

                                                 aaatcttttGTTTTTCTCTTTCGGTCATCTAAACTGCCCGCTTTGCTACCCGGGTTTTTTTTTTGTTGTTTTTTTCCACCCTCTAGTAAGAAAATGCGCACCGGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                      <  1:1350851/144‑1 (MQ=255)

                                                      gTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccg                                                                                                 <  2:1359285/150‑1 (MQ=255)

                                                       tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACAGATTTACCGTCTTTAAGATAGTacgc                                                                                                >  1:6136426/1‑150 (MQ=255)

                                                        ttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgcc                                                                                               <  2:2594893/150‑1 (MQ=255)

                                                        ttGCTTCTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgcc                                                                                               <  2:3882831/150‑1 (MQ=255)

                                                             tttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTc                                                                                          <  2:3920037/150‑1 (MQ=255)

                                                                 tctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCa                                                                                      >  1:2091691/1‑150 (MQ=255)

                                                                    ttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc                                                                                   <  2:4053753/150‑1 (MQ=255)

                                                                     tCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAcc                                                                                  >  1:6441437/1‑150 (MQ=255)

                                                                       ggTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCt                                                                                >  2:869125/1‑150 (MQ=255)

                                                                        gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTg                                                                               >  1:571139/1‑150 (MQ=255)

                                                                         tCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑GTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGt                                                                              >  1:3210787/1‑150 (MQ=255)

                                                                             aaaaaaTCCCCCCTTTCAAACCACGGCATTTTTTACCATCTCTGTAACCACCCTCCACAAAAAAAATGCCCAACCGC‑TCCCCCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgc                                                                          <  1:3058375/148‑1 (MQ=255)

                                                                              tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgcc                                                                         >  2:229179/1‑150 (MQ=255)

                                                                                aaaTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgccgc                                                                       <  2:5865792/150‑1 (MQ=255)

                                                                                  aTGCCCGCTTTGCTACCTCGGTATTGTTTTCGTTGTCTGTTGCCACCCTCGAGTCATAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCt                                                                     <  1:6285819/150‑1 (MQ=255)

                                                                                  aTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCt                                                                     <  2:2870937/150‑1 (MQ=255)

                                                                                     cccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGT‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTc                                                                  <  2:568002/150‑1 (MQ=255)

                                                                                      ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCa                                                                 <  2:484819/150‑1 (MQ=255)

                                                                                         cTTTGCTACCTGGGGTTTGTTTTCGTTTTCTGTTTCCCCCCTGGAGTCATAAAATGCG‑AACCGGTTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCa                                                              <  1:4600193/150‑1 (MQ=255)

                                                                                          tttGCTACCACGTTTTTGTTTTCTATGTCTTTTTCCACCCTCGAGACATAAAATGCGCAACCGT‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCAt                                                             <  1:3284581/150‑1 (MQ=255)

                                                                                          tttGCTACCACGGGTTTGTTTTCGAGGTCTGTTTCCCCCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCAt                                                             <  2:1883095/150‑1 (MQ=255)

                                                                                            ggcgaCCCCGGGATTGTTGTCGAGGTCTGTTTCCCCCCTCGAGTCCGAAAAGGCGCCCCCGG‑TGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATct                                                           <  1:3940000/146‑1 (MQ=255)

                                                                                               tACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACAGATTTACCGTCTTGAAGATAGTCCGCCGCGTCCTCAAACCCAGTTTGGCGCCAGTCAGCATctgct                                                        >  1:4243405/1‑150 (MQ=255)

                                                                                               tACCACGGGATTGGTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGGTTGCCGCCAGTCAGCATctgct                                                        >  2:4412666/1‑150 (MQ=255)

                                                                                                aCCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGAGTTACCGTCTTGTAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgctg                                                       >  2:1687594/1‑150 (MQ=255)

                                                                                                   aCGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCG‑TTGCGCAGACACAAAAAAACGAACAGATTTAACCTCTTGAAGATAGTCCGCCCCATCCTCAAAACCAGTTTGCCGCCTGTCAGCATCTGCTGGCg                                                    >  1:2581059/1‑150 (MQ=255)

                                                                                                           gTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCCG‑GTGCGCAGACACAAAACACCCACCGGAGTTACCCGCTTGAAGATCGTCCGCCGCAGCCACACACCCCGTTTGCCGCCAGTCAGCATCTGCGGGCGCTGGTCgg                                            >  1:789611/1‑149 (MQ=255)

                                                                                                            ttttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                           >  2:3404287/1‑150 (MQ=255)

                                                                                                            ttttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                           >  2:6106181/1‑150 (MQ=255)

                                                                                                             tttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTg                                          <  2:4200190/150‑1 (MQ=255)

                                                                                                             tttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCG‑TTTCGCAGAATAAAAAAAGCGAAAGGATTTAACCTCTTTAATAAAGTACGCCGCGTCCACAAACCCAGTTTTCCGCCAGTCAGCAACAGCTGGCGATGTTCGTTg                                          >  1:5174959/1‑150 (MQ=255)

                                                                                                              ttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGc                                         >  2:727777/1‑150 (MQ=255)

                                                                                                               tCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCg                                        <  2:3553689/150‑1 (MQ=255)

                                                                                                                ctaTGTCTTTTACCACCAACTATTCATAAAATGCTAAACCGT‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                       <  1:2688071/148‑1 (MQ=255)

                                                                                                                cGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                       <  2:2974421/150‑1 (MQ=255)

                                                                                                                  aTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCCTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTg                                     >  2:2422655/1‑150 (MQ=255)

                                                                                                                   tGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGa                                    >  2:1118538/1‑150 (MQ=255)

                                                                                                                    gTCTGTTCCCACCCTCGAGACAGAAAAGGCGAAACGGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGAt                                   <  1:5328567/150‑1 (MQ=255)

                                                                                                                      cTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCa                                 <  2:4624025/150‑1 (MQ=255)

                                                                                                                       tGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCAc                                <  2:2053842/150‑1 (MQ=255)

                                                                                                                       tGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCAc                                >  2:5996035/1‑150 (MQ=255)

                                                                                                                       tGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCAc                                >  2:1578839/1‑150 (MQ=255)

                                                                                                                         tttCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGc                              >  2:6333441/1‑150 (MQ=255)

                                                                                                                           tCCACCCTCTATTCATAAAATGCTCAACGGG‑TTGCGCAGACTCAAAACAGCGTACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCa                            <  1:5951808/150‑1 (MQ=255)

                                                                                                                              aCCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAATATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCAGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcg                         >  2:6222812/1‑150 (MQ=255)

                                                                                                                              aCCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATAGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcg                         >  2:1807942/1‑150 (MQ=255)

                                                                                                                                 cTCGAGTCAGAAAATGCGCACCCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTa                      <  1:172834/150‑1 (MQ=255)

                                                                                                                                 cTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTa                      >  2:3495386/1‑150 (MQ=255)

                                                                                                                                 cTCGAGACAGAAAATGCGCACCCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTa                      <  1:6612096/150‑1 (MQ=255)

                                                                                                                                    gAGTCAGACAAGGCGCACCGGG‑TGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTa                   <  1:1214048/150‑1 (MQ=255)

                                                                                                                                    gAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTa                   <  1:1645949/150‑1 (MQ=255)

                                                                                                                                      gTCAGAAACTGCGCAACCGG‑TTGCGCAGACTCAAACCAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAAt                 <  2:4737568/150‑1 (MQ=255)

                                                                                                                                      gTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAAt                 <  2:3792533/150‑1 (MQ=255)

                                                                                                                                           aaaaTGCGCCCGGGG‑TGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGcc            <  1:101908/150‑1 (MQ=255)

                                                                                                                                           aaaaTGCGCAACCGG‑TTGCGCAGACTCATAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGcc            <  1:3594509/150‑1 (MQ=255)

                                                                                                                                                  gcAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCt     >  2:2689022/1‑150 (MQ=255)

                                                                                                                                                  gcAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCt     <  2:1034699/150‑1 (MQ=255)

                                                                                                                                                  gcAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCt     >  1:4556862/1‑150 (MQ=255)

                                                                                                                                                   cAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg    >  1:3507082/1‑150 (MQ=255)

                                                                                                                                                     aCCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTggg  <  2:3693474/150‑1 (MQ=255)

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGG‑TTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828831‑3829129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`