BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,828,980	A→C	9.2%	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	9.2%	503.5 / 5.1	205	Q12H (CAA→CAC)	xanP	xanthine permease
Reads supporting (aligned to +/- strand): ref base A (121/60); new base C (5/19); total (126/79)
Fisher's exact test for biased strand distribution p-value = 2.36e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 211 total aligned reads displayed.

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGG  >  NC_000913/3828831‑3829128

 tataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAc                                                                                                                                                     >  1:301709/1‑150 (MQ=255)

   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCg                                                                                                                                                   >  2:3842390/1‑150 (MQ=255)

   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                   <  2:1204056/150‑1 (MQ=255)

   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                   <  1:2619203/150‑1 (MQ=255)

     aaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGt                                                                                                                                                 >  1:1222076/1‑150 (MQ=255)

      aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGtt                                                                                                                                                <  1:3454361/150‑1 (MQ=255)

      aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGtt                                                                                                                                                >  2:2894142/1‑150 (MQ=255)

           ggCATTTATCTCAGTTCCAGATTGAGATGAGGATGGGGCAACCGTTTGCTTTTCTCTGTCGGTCAGCTCAACTGCCCGCTGTGCTACCACGGGATTGTTCTCGATGGCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCa                                                                                                                                           <  2:4085833/149‑1 (MQ=255)

               tttATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACt                                                                                                                                       >  2:2006862/1‑150 (MQ=255)

                   tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAcaaaaa                                                                                                                                   >  2:1773014/1‑145 (MQ=255)

                     tcATTTAAATCTTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTGTCCACCCTCGAGTCAGAAACTGCGCAACCGGGTGCGCAGACTCAAAAc                                                                                                                                 >  2:2719395/1‑150 (MQ=255)

                     tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATTCGCAACCGGGTGCGCAGACTCaaaaa                                                                                                                                 >  1:3472034/1‑149 (MQ=255)

                     tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAc                                                                                                                                 >  1:123018/1‑150 (MQ=255)

                     tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACAAAAc                                                                                                                                 >  1:1809726/1‑150 (MQ=255)

                       aTTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAg                                                                                                                               >  1:3414983/1‑150 (MQ=255)

                        tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                                                                                              >  2:2069041/1‑150 (MQ=255)

                           aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACACAAAACAGCGaa                                                                                                                           >  1:2704033/1‑150 (MQ=255)

                           aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTAGATGTCTGTTTCCACCCACGAGTCAGAAAATTCGCAACAGGTTGCGCAGAATCAAAACACCGaa                                                                                                                           >  1:4090850/1‑150 (MQ=255)

                                  tAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGAGCACACACAAAAAAGCGAACAGAttt                                                                                                                    >  1:1844802/1‑150 (MQ=255)

                                   aGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAACATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAATTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTa                                                                                                                   <  2:1056306/150‑1 (MQ=255)

                                       gAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGt                                                                                                               >  2:52568/1‑150 (MQ=255)

                                        aGGTTGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTAACACGGGATTGTTTTCGATGTGTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTGAAAAATGGGTACTGTTTTACCGtg                                                                                                              >  1:2765703/1‑149 (MQ=255)

                                           ttgtgTCAATCTTTTGCTTTTCTCTTTGGGTCAGCTAAATGGCGGGCTTTGCTACCAGGGGTTTTTTTTCGTTGTCTGTTTCCACCCTCGATTCATAAAATGCGCACCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTg                                                                                                           <  1:3043532/149‑1 (MQ=255)

                                           aTGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTg                                                                                                           >  2:562801/1‑150 (MQ=255)

                                           aTGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTg                                                                                                           >  2:3599548/1‑150 (MQ=255)

                                                 cAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                     <  2:449850/150‑1 (MQ=255)

                                                 cAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                     <  2:1365075/150‑1 (MQ=255)

                                                    tcGTTTGCTTTTCTCTTTCGGTCACATAACATGCCCGCTTTGCTACCTGGGTTTTGTTTTCTATGTCTGTTTCCCCCCTCGAGCCACAAAATGCGCAACGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTc                                                                                                  <  1:2239552/149‑1 (MQ=255)

                                                    aCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTTCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTc                                                                                                  <  1:2398780/150‑1 (MQ=255)

                                                      gttggtttttCTCTTTGGGTCAGCTAAATTGCCCGCTTTGCTACCAGGGTTTTTTTTTCGTTGTCTGTTTCCACCCTCGCGTCATAAATTGCGCACCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccg                                                                                                <  1:1275545/144‑1 (MQ=255)

                                                       tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                                                                                               <  2:2348334/150‑1 (MQ=255)

                                                       tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                                                                                               >  2:435250/1‑150 (MQ=255)

                                                       tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACACCGAACTGATTTAACGTCTTGAAGATCGTccgc                                                                                               >  1:752147/1‑150 (MQ=255)

                                                            ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATAGTCCGCCGCtt                                                                                          >  2:387841/1‑150 (MQ=255)

                                                               tctcTTTGGGTCAGCTAAATGGCCCGCTTTGCTACCCCGGGTTTGTTTTCGTTGTCTGTTTCCCCCCTCGAGTCAGAAATGGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCt                                                                                       <  1:3764178/150‑1 (MQ=255)

                                                                ctctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGGTGCGCAGACTCAAAAAAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTc                                                                                      >  1:2045919/1‑150 (MQ=255)

                                                                 tctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGTTTTACCGTCTTGAAGTTCGTCCGCCGCGTCCGCa                                                                                     >  1:3291727/1‑150 (MQ=255)

                                                                 tctTTCGGACAGCTAAAATGCCCGCTTTGATACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCa                                                                                     <  2:1440029/150‑1 (MQ=255)

                                                                  ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaa                                                                                    >  2:2856180/1‑150 (MQ=255)

                                                                   tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                   >  1:3551881/1‑150 (MQ=255)

                                                                     tCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACGGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCGCAAAcc                                                                                 >  1:2932921/1‑150 (MQ=255)

                                                                          cAGATAAAATGCCCGCTTTGCTACCCCGGGATTGTTTTCGTTGTCTGTTCCCACCCTCGAGACAGAAATTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGtt                                                                            <  1:1416350/150‑1 (MQ=255)

                                                                            gCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTg                                                                          >  2:3163738/1‑150 (MQ=255)

                                                                             cTAAATTCCCCGGTGTGCTACGACGGTATTGTTTTCTTTGTCTGTTTCCACCCTCGACGATCAACTTGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgc                                                                         <  1:2524710/150‑1 (MQ=255)

                                                                             atAAAAGGCCCGCTTTGCTACCAGGGGTTTGTTTTCGCTGTCTGTTACCACCCCCGCGTCAGAAAATGCGCAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgc                                                                         <  1:119337/149‑1 (MQ=255)

                                                                              tataaTGCCCGCTTTGCTACCAGGGGATTGTTTTCGAGGTCTGTTTCCACCCTCGAGTCATAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgcc                                                                        <  1:2337519/147‑1 (MQ=255)

                                                                              tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACTGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgcc                                                                        >  1:160257/1‑150 (MQ=255)

                                                                              tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACTGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgcc                                                                        >  1:2359927/1‑150 (MQ=255)

                                                                              tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgcc                                                                        >  2:2593741/1‑150 (MQ=255)

                                                                               aatttGCCCGCGTTGCTACCACGGTATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCCGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgccg                                                                       <  1:4025126/146‑1 (MQ=255)

                                                                               aaaaTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACACAAAAAAGCGAACTGATTTACCGTCTTGTAGATAGTCCGCCGCGTCCTCAAACCCGGGTTgccg                                                                       >  1:3823496/1‑150 (MQ=255)

                                                                               aaaaTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAACCTGTTTgccg                                                                       >  2:3957159/1‑150 (MQ=255)

                                                                                  acaccagCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGATTTAACGTCTTTAATATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCt                                                                    >  1:281094/7‑150 (MQ=255)

                                                                                  aTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAAGGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCt                                                                    <  2:3467620/150‑1 (MQ=255)

                                                                                      ccGCTTTGCTACCAGGGTATTGTTTTCGATGTCTGTGTCCACCCGCGAGACAGAAATTGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCa                                                                <  1:831926/150‑1 (MQ=255)

                                                                                      ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCAGTTTGCCGCCTGTCa                                                                >  2:2400222/1‑150 (MQ=255)

                                                                                      ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGATTTAACCTCTTTAAGATCGTCCGCCGCTTCCACAAACCCAGTTTGCCGCCTGTCa                                                                >  1:56786/1‑150 (MQ=255)

                                                                                       cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTTTTTGCCGCCTGTCAg                                                               >  2:756568/1‑150 (MQ=255)

                                                                                         cTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCACACACAAAACAGCGAACTGATATACCCTCTAGAAGAAAGTCCCCCCCTTCCTCAAACCCTGTATGCCGCCTGACAGCa                                                             >  1:2229509/1‑150 (MQ=255)

                                                                                            tGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATct                                                          <  1:2633570/150‑1 (MQ=255)

                                                                                              cTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgc                                                        >  2:3989519/1‑150 (MQ=255)

                                                                                               tACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACACAAAAAAGCGAACAGATTTAACGTCTTGAAGATAGTCCGCCGCGTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgcg                                                       >  1:51807/1‑149 (MQ=255)

                                                                                                 ccACGGTTTTGTTTTCGCTGTCTGTTTCCCCCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCGGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTgg                                                     <  1:3175933/150‑1 (MQ=255)

                                                                                                 ccACGGGTTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTgg                                                     <  2:668537/150‑1 (MQ=255)

                                                                                                 ccACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGGTGCGCAGACACAAAACAGCCAACAGATTTAACGTCTTGAAGATCGTACGCCCCGTCCACAAACCCAGTTTTCCGCCAGTAAGCAtaagcagg                                                     >  1:1650663/1‑143 (MQ=255)

                                                                                                  cACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGc                                                    >  2:705356/1‑150 (MQ=255)

                                                                                                    cGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTTCGCAGACTCAAAAAATCGAACAGATTTAACGTCTTTAAGATAGTCCGCCGCGTACTCAAACACAGTTTTCCGCCAGTAAGCATCTGAAGGCGa                                                  >  1:3671270/1‑150 (MQ=255)

                                                                                                      ggATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGATGGCGATg                                                >  2:52656/1‑150 (MQ=255)

                                                                                                      ggATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCGTGAAGATCGTCCGCCGCTTCCTCAAACCCGGTTTGCCGCCTGTCAGCATCAGCTGGCGATg                                                >  2:2548762/1‑150 (MQ=255)

                                                                                                        aTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGtt                                              >  2:2486572/1‑150 (MQ=255)

                                                                                                         ttGTTTTCGATTTCTGTTCCCACCCTCGAGTCATAAAATGCGCAACGTGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTc                                             <  1:644733/150‑1 (MQ=255)

                                                                                                         ttGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTc                                             <  1:1345129/150‑1 (MQ=255)

                                                                                                         tggTTTTTGTTGTCTTTTTCCACCCTCGCGTCAGATTATGCGCAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTc                                             <  1:3866629/148‑1 (MQ=255)

                                                                                                            ttttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                          >  2:1654255/1‑150 (MQ=255)

                                                                                                                cGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                      <  2:3783544/150‑1 (MQ=255)

                                                                                                                cGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                      <  1:2894142/150‑1 (MQ=255)

                                                                                                                cGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                      >  2:2099379/1‑150 (MQ=255)

                                                                                                                   tGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGa                                   <  2:3030267/150‑1 (MQ=255)

                                                                                                                   tGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGa                                   <  2:1645231/150‑1 (MQ=255)

                                                                                                                         tttCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGc                             >  2:3909377/1‑150 (MQ=255)

                                                                                                                         tttCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGc                             >  2:3616419/1‑150 (MQ=255)

                                                                                                                         tttCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAATCAAAACAGCGAACTGATTTACCGTATTGAAGATCGTCCGCCGCTTCCACAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGc                             >  2:3909240/1‑150 (MQ=255)

                                                                                                                           tCCACCCTCGAGTCAGAAAAGGCGCAACCGGTGGCGCAGACTCTAAACAGCGACCTGATTTACCGTCTTGAAGATCGTCCGCCCCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCa                           >  2:806608/1‑150 (MQ=255)

                                                                                                                                 cTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGACTTTAAGATCGTACGCCGCTTCCTCAAACCCGGTTTTCCGCCTGTCAGCATCAGCAGGCGATGTTCGTTGCCGTGATCACGCCAGCGCTa                     >  2:1257314/1‑150 (MQ=255)

                                                                                                                                  tCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTAt                    >  2:3307882/1‑150 (MQ=255)

                                                                                                                                  tCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTAt                    >  2:1243623/1‑150 (MQ=255)

                                                                                                                                  tCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTAt                    >  2:1327070/1‑150 (MQ=255)

                                                                                                                                  tCGAGACAGAAAATGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTAt                    <  1:2134545/150‑1 (MQ=255)

                                                                                                                                   cGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTAtt                   <  2:301709/150‑1 (MQ=255)

                                                                                                                                     acacacaaaaTCCGCAACCGGTAGCCCAGACACAAAACAGCGAACTGATATACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTaa                 <  2:1865268/144‑1 (MQ=255)

                                                                                                                                       tCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATc               >  2:1269455/1‑150 (MQ=255)

                                                                                                                                         aGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTg             <  1:4060558/150‑1 (MQ=255)

                                                                                                                                            gtctGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa          >  2:3743431/4‑150 (MQ=255)

                                                                                                                                            aaaTGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa          <  1:2561498/150‑1 (MQ=255)

                                                                                                                                            aaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa          >  2:1508482/1‑150 (MQ=255)

                                                                                                                                                  gcACCCGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCt    <  2:2609428/150‑1 (MQ=255)

                                                                                                                                                   cAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCCGCGCTATTAATCTGCCAGgcgcgg   >  2:1129094/1‑148 (MQ=255)

                                                                                                                                                   cAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg   <  2:1580419/150‑1 (MQ=255)

                                                                                                                                                    aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg  >  2:2808832/1‑150 (MQ=255)

                                                                                                                                                    aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg  <  2:3511355/150‑1 (MQ=255)

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGG  >  NC_000913/3828831‑3829128

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