BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,828,980	A→C	6.9%	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	6.9%	2117.5 / 23.5	814	Q12H (CAA→CAC)	xanP	xanthine permease
Reads supporting (aligned to +/- strand): ref base A (400/336); new base C (12/64); total (412/402)
Fisher's exact test for biased strand distribution p-value = 6.51e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 834 total aligned reads displayed.

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828831‑3829129

 tataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAc                                                                                                                                                       >  1:9636393/1‑150 (MQ=255)

  ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                      <  2:5629725/150‑1 (MQ=255)

  ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGAAAcc                                                                                                                                                      >  2:1728230/1‑150 (MQ=255)

   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                     >  2:9807997/1‑150 (MQ=255)

   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                     >  2:9485902/1‑150 (MQ=255)

   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTATATTTCGGTCAGCTAAAATGACCGCTTTGCTACCACGGGATTGTTTTCGATTTCTGTTTCCACCCACGAGTCAGAAAATGCGCAACCg                                                                                                                                                     >  1:14900122/1‑150 (MQ=255)

    aaaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCgg                                                                                                                                                    >  2:3945706/1‑150 (MQ=255)

        aaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgc                                                                                                                                                <  1:487763/150‑1 (MQ=255)

           tGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCa                                                                                                                                             >  1:6535595/1‑150 (MQ=255)

           tGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCa                                                                                                                                             >  1:3222462/1‑150 (MQ=255)

            gCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCat                                                                                                                                            <  2:2173521/150‑2 (MQ=255)

             cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCatt                                                                                                                                           <  2:5217365/150‑3 (MQ=255)

              aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGAc                                                                                                                                          >  1:13562118/1‑150 (MQ=255)

              aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAc                                                                                                                                          <  2:2067429/150‑1 (MQ=255)

              aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAc                                                                                                                                          <  2:7355238/150‑1 (MQ=255)

               tttATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAc                                                                                                                                          >  2:8746277/1‑150 (MQ=255)

                 tATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattttct                                                                                                                                       <  2:3807226/150‑7 (MQ=255)

                  tcctcATCTAAATATTTAAATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaa                                                                                                                                      >  2:6071495/3‑150 (MQ=255)

                   tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaa                                                                                                                                     >  2:14242672/1‑150 (MQ=255)

                     tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAc                                                                                                                                   >  1:15433910/1‑150 (MQ=255)

                     tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACaaaaa                                                                                                                                   >  1:3983471/1‑149 (MQ=255)

                      cATTTAAATATTTAGATGTGGTTGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGTATTGTTTTCGTTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGAACCCGGTTGCGCAGACTCAAAACa                                                                                                                                  <  2:6609960/150‑1 (MQ=255)

                        tttAAATATTTAGATGAGGATGTGGCAAACGTATGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCAACCACGGGATTGTTTTCTATCTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGCAATCGCAAACTCAAAACAGc                                                                                                                                >  1:8288473/1‑150 (MQ=255)

                           aaaTTTTTAGATGAGGATGTGGCAATCGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGaa                                                                                                                             <  2:8424274/150‑1 (MQ=255)

                             atatTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACt                                                                                                                           <  2:6484225/150‑1 (MQ=255)

                               atTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGa                                                                                                                         <  2:11992746/150‑1 (MQ=255)

                               atTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGa                                                                                                                         <  2:7791565/150‑1 (MQ=255)

                                tttAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAt                                                                                                                        >  2:3451260/1‑150 (MQ=255)

                                       gAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGATTTACCGt                                                                                                                 >  1:13085388/1‑150 (MQ=255)

                                               ggCAAACTTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGa                                                                                                         <  2:7056700/150‑1 (MQ=255)

                                               ggCAAACGTTTGC‑TTTTCTCTGTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTTTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAATTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGa                                                                                                         <  2:14074892/150‑1 (MQ=255)

                                                gCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGAt                                                                                                        >  2:9192898/1‑150 (MQ=255)

                                                gCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGAt                                                                                                        >  2:12757198/1‑150 (MQ=255)

                                                gCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGAATCAAAACAGCGAACAGATTTAACGGCTTGAAGAt                                                                                                        >  1:11259017/1‑150 (MQ=255)

                                                         tGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgccg                                                                                               >  2:4157654/1‑150 (MQ=255)

                                                           c‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACAAAAAAGCGAACTGATTTACAGTCATGAAGATCGTCCGCCGCt                                                                                             >  1:5931816/1‑150 (MQ=255)

                                                             ttttCTCTTTCGGTCCGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAACGAACTGCTTTACCGTCTTTAAGATCGTCCGCCGCtg                                                                                            >  1:10582344/1‑149 (MQ=255)

                                                             ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt                                                                                            >  2:7561131/1‑150 (MQ=255)

                                                             ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt                                                                                            >  2:1306239/1‑150 (MQ=255)

                                                             ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt                                                                                            >  1:12807852/1‑150 (MQ=255)

                                                              tttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTc                                                                                           >  2:13571716/1‑150 (MQ=255)

                                                                tctcTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCt                                                                                         <  2:10798967/150‑1 (MQ=255)

                                                                 ctctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTc                                                                                        >  2:12408710/1‑150 (MQ=255)

                                                                  tctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑a                                                                                      >  2:6223024/1‑150 (MQ=255)

                                                                   ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑aa                                                                                     >  1:2338600/1‑150 (MQ=255)

                                                                    tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAAAACCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTC‑aaa                                                                                    >  1:4538736/1‑150 (MQ=255)

                                                                    tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑aaa                                                                                    <  2:14430173/150‑1 (MQ=255)

                                                                            atcTAAAATGCCAGCTTTGATACCACGGTATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCATAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATAGTCCGCCGCTTCCTC‑AAACCCTGttt                                                                            <  2:5450893/148‑1 (MQ=255)

                                                                             gCTAAAATGCCCGCTTTGCTACCACGGGTTTGTTTTCGTTGTCTGTTTCCACCCTCGAGTCAGAAATTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTg                                                                           <  1:9839657/150‑1 (MQ=255)

                                                                                aaaaTGCCCGCTTTGCTCCCACGGTTTTGTTTTCGTTGTGTGTGTCCACCCTCGAGTAAGAAACGGCGCACCCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCTTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTgccg                                                                        <  1:8936155/150‑1 (MQ=255)

                                                                                 aattGCCCGCTTTGCTACCAGGGTTTTTTTTTCTTTTTTTGTCTCCACCCTCGTTTCATAAAATGCGCAACGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTgccgc                                                                       <  1:14316748/147‑1 (MQ=255)

                                                                                 aaaTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTgccgc                                                                       >  2:5776618/1‑150 (MQ=255)

                                                                                   aggCCCGTTTTGCAACCAGGGTATTGTTTTCGAGGTCTGTTCCCACCCTCGTGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCt                                                                     <  1:2747054/148‑1 (MQ=255)

                                                                                   aTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCt                                                                     <  2:3138022/150‑1 (MQ=255)

                                                                                      cccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTc                                                                  <  2:8239710/150‑1 (MQ=255)

                                                                                       ccGCTTTGCTACCAGGGGTTTGTTTTCTATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCa                                                                 <  1:13831364/150‑1 (MQ=255)

                                                                                       ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCa                                                                 >  2:10996281/1‑150 (MQ=255)

                                                                                         gCTTTGCTACCAGGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGc                                                               <  2:690430/150‑1 (MQ=255)

                                                                                            ttGCTACCACGGTTTTTTTTTCGTTGTCTGTTTCCACCCTCTATTAACAAAATGCGCAACGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATc                                                            <  1:15260665/150‑1 (MQ=255)

                                                                                               cTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATctgc                                                         >  2:4468863/1‑150 (MQ=255)

                                                                                                tACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGCTTGAGCAGACTCAAAACAGCGAACGGATTAAACGGAATGAAGAACGTCCGCCCCGTCCAC‑AAAACCAGTTTGACGCCAGTGAGCATcagat                                                        >  1:12271661/1‑146 (MQ=255)

                                                                                                tACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATctgct                                                        >  2:2623554/1‑150 (MQ=255)

                                                                                                 aCCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATctgctg                                                       >  2:1252992/1‑150 (MQ=255)

                                                                                                 aCCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATctgctg                                                       <  1:7682494/150‑1 (MQ=255)

                                                                                                  ccACGGGTTTGTTTTCGATGTCTGTTTCCACCCTCGCTTCAGAAAATGCGCATCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTgg                                                      <  2:14845060/150‑1 (MQ=255)

                                                                                                  ccACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACACAAAACAGCGAACAGATTTAACGGCTTGAATATAGTCCGCCGCTTCCTC‑AAACCCAGTTTGCCGCCGGTCAGCAACTGCTgg                                                      >  1:3756125/1‑150 (MQ=255)

                                                                                                   cACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAAATGATTTACCGTCTTGAAGATCGTCCGCCCCGTCCAC‑AAACCCAGTTTTCCGCCTGTCAGCAACTGCTGGc                                                     >  1:847136/1‑150 (MQ=255)

                                                                                                   cACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAACCGAAATGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTC‑AAACCCAGTTTTCCGCCTGTCAGCATATTCTGGc                                                     >  1:10056088/1‑150 (MQ=255)

                                                                                                   cACGGGATTGTTTTCGATGTCTGTTTACACCCTCGAGTCAGAAAATGCGCAACCCGGTTAGCAGACAAAAAAAATCGAAAAGATTTAACTTCTTTAAGATCGTCCGCCGCATCCTC‑AAAACCTGTTTGCCGCAAGTCAGCATCTGCAGGc                                                     >  1:11019966/1‑150 (MQ=255)

                                                                                                     cGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGa                                                   >  2:13312209/1‑150 (MQ=255)

                                                                                                      ggtaTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTATTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGAt                                                  <  2:8952652/147‑1 (MQ=255)

                                                                                                      gggATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCCCCCGGTGCGCCGACACAAAAAAGCGAACAGATTTACAGTCTTTTATATCGTCCGCCGCGTCATCAAAACCAGGTTTGCC‑CCGGGCAGCAACTGCTGGCggg                                                  >  1:5739217/1‑148 (MQ=255)

                                                                                                         aTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACAAAAAAAAGCGAACTGATTTAACGTCTTGAAGATCGTCCGCCCCGTCCTC‑AAACCCAGTTTTCCGCCAGTCAGCAACTGCTGGCGATGtt                                               >  1:14144763/1‑150 (MQ=255)

                                                                                                          tttttttCGTTGTCTGTTTCCCCCCCCGAGACAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTc                                              <  1:5524134/147‑1 (MQ=255)

                                                                                                          ttGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTc                                              <  1:8622545/150‑1 (MQ=255)

                                                                                                           tGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCg                                             <  2:15091768/150‑1 (MQ=255)

                                                                                                             ttttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                           <  1:15552243/150‑1 (MQ=255)

                                                                                                              tttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTg                                          >  2:11098513/1‑150 (MQ=255)

                                                                                                                tCGTTGTCTGTTACCACCCTCGAGCCAGAAAAGGCGCACCCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCg                                        <  1:7009817/150‑1 (MQ=255)

                                                                                                                tCGATGTCTGTTTCCCCCCTCGAGTCAGAAAATGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCg                                        <  1:1838142/150‑1 (MQ=255)

                                                                                                                 cgttGTCTGTTTCCACCCTCGAGTCATAAATTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                       <  1:9169333/147‑1 (MQ=255)

                                                                                                                  gATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCCGTTGCGCAGACTCAAAACAGCGAACAGAATTTACGTCTTGAAGATAGTCCGCCCCGTCCTC‑AAAACCAGTTTTGCGCCAGTCAGCCTCAGCTGGCGATGTTAGTTGCggg                                      >  1:6149796/1‑149 (MQ=255)

                                                                                                                      tCTTTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATc                                  <  1:5173535/150‑1 (MQ=255)

                                                                                                                      tCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATc                                  >  2:2377883/1‑150 (MQ=255)

                                                                                                                       cTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCa                                 >  2:10269985/1‑150 (MQ=255)

                                                                                                                           ttCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGcc                             >  2:2876303/1‑150 (MQ=255)

                                                                                                                            tCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCa                            >  2:13738278/1‑150 (MQ=255)

                                                                                                                             ccACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                           >  2:8088073/1‑150 (MQ=255)

                                                                                                                              cACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgc                          >  2:11586345/1‑150 (MQ=255)

                                                                                                                               aCCCTCGAGTCAGAAAAGGCGCACCCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcg                         <  1:9337300/150‑1 (MQ=255)

                                                                                                                                       gTCAGAAATTGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAAt                 <  2:10558096/150‑1 (MQ=255)

                                                                                                                                       gTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAAt                 <  1:1078972/150‑1 (MQ=255)

                                                                                                                                        tCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATc                <  1:790246/150‑1 (MQ=255)

                                                                                                                                        acAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATc                <  1:1158249/149‑1 (MQ=255)

                                                                                                                                         cAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACGGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCt               >  2:4862332/1‑150 (MQ=255)

                                                                                                                                          aGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTg              >  2:10089588/1‑150 (MQ=255)

                                                                                                                                           gaaatggcgcACCCGGTGGCGCAGACTCAACACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGc             <  1:10376872/144‑1 (MQ=255)

                                                                                                                                            aaaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGcc            <  1:7409939/150‑1 (MQ=255)

                                                                                                                                                tGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGgc        <  2:15467501/150‑1 (MQ=255)

                                                                                                                                                      aCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTggg  <  1:11983040/150‑1 (MQ=255)

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGC‑TTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTC‑AAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  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: `‑`