BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,828,980	A→C	6.5%	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.5%	512.9 / 3.3	215	Q12H (CAA→CAC)	xanP	xanthine permease
Reads supporting (aligned to +/- strand): ref base A (107/88); new base C (3/16); total (110/105)
Fisher's exact test for biased strand distribution p-value = 1.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 221 total aligned reads displayed.

ATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828833‑3829129

ataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                                                                                                     >  1:3909738/1‑150 (MQ=255)

 taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                    <  1:2766601/150‑1 (MQ=255)

 taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                    <  2:732839/150‑1 (MQ=255)

  aaaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCgg                                                                                                                                                   >  1:3285991/1‑150 (MQ=255)

        gTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgcgc                                                                                                                                             >  2:833442/1‑150 (MQ=255)

           cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGa                                                                                                                                          >  1:1429506/1‑150 (MQ=255)

           cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGa                                                                                                                                          <  2:3274364/150‑1 (MQ=255)

           cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGa                                                                                                                                          <  2:3399159/150‑1 (MQ=255)

           cATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGa                                                                                                                                          >  1:787783/1‑150 (MQ=255)

            aTTTATCTCATGTAAATATTTAGATGCGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTAGAGCAGAc                                                                                                                                         >  2:3365194/1‑150 (MQ=255)

                aTCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaa                                                                                                                                     <  2:196300/150‑1 (MQ=255)

                 tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaa                                                                                                                                    >  1:1787952/1‑150 (MQ=255)

                  ctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTATTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaaa                                                                                                                                   >  2:2439960/1‑150 (MQ=255)

                   tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaaaa                                                                                                                                  >  1:1899116/1‑149 (MQ=255)

                   tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAc                                                                                                                                  >  1:2574923/1‑150 (MQ=255)

                   tcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGGTGCGCAGACTCAAAAc                                                                                                                                  >  1:1405757/1‑150 (MQ=255)

                      tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCCCCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCCGACTCAAAACAGc                                                                                                                               >  1:1449254/1‑150 (MQ=255)

                      tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                                                                                               >  2:3024703/1‑150 (MQ=255)

                      tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                                                                                               >  2:3776752/1‑150 (MQ=255)

                      tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                                                                                               >  2:2764191/1‑150 (MQ=255)

                      tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGGTGCGCAGACTCAAAACAGc                                                                                                                               >  2:2998790/1‑150 (MQ=255)

                        tAAATATTTTGTTGTGGTTGTGTCTAACTTTGGCTTTTCTCTTTCGGTCAGATAAACTGCCCGATTTTATACCCGGGTTTTTTTTTGG‑TTGTCTTTTTCCACCCTCTAGTCATAAATGGCGCAACGGGTGGCGCAGACTCAAAACAGCGa                                                                                                                             <  1:2088866/150‑1 (MQ=255)

                              tttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAt                                                                                                                       >  2:2298424/1‑150 (MQ=255)

                               ttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAtt                                                                                                                      <  2:433596/150‑1 (MQ=255)

                                 aGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAATATGCGCCACCGGTTGCGCAGACTCCAAACAGCGAACTGCtttt                                                                                                                    >  2:2394991/1‑149 (MQ=255)

                                       ggttgtgGCAAAGTTTTGCTTTTCTCTTGCGGTCAGATAAAATGCCCGCTTTGCTACCAGGGGTTTTTTTTCG‑TTGTCTTTTACCCCCCTCGATTAATAAAATGCGAAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCt                                                                                                              <  1:3233359/147‑1 (MQ=255)

                                        gATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTTCGCAGAATCAAAACAGCGAACTGATTTACCGTCtt                                                                                                             >  1:2617308/1‑150 (MQ=255)

                                          gggggCAACCGTTTGGTTTGCTCGTTCGGTCCGCTACACGGCCCGCTTTGCTCCCCCGGGTTTGTGTTCG‑ATGTCTTCTGCCCCCCGCGATCCAGAAAATGCGCACCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGa                                                                                                           <  1:2435664/147‑1 (MQ=255)

                                           gtgGCTAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaa                                                                                                          <  2:316960/150‑1 (MQ=255)

                                           ggggCAAACGTTTGCTGTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaa                                                                                                          <  2:1457214/148‑1 (MQ=255)

                                             ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCACCAGTCAGAAAAAGCGCCACCCGTTGCGCACACACAAAAAAGCGAACAGATTTAACCTCTTTaaaa                                                                                                        >  1:2792439/1‑148 (MQ=255)

                                              gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTTATTTACCGTCTTGAAGAt                                                                                                       >  2:2687246/1‑150 (MQ=255)

                                                aaaCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACACAAAACATCGAAATGATTTACAGTCTTGAATATCg                                                                                                     >  1:3063988/1‑150 (MQ=255)

                                                  atttttGCTTTTCTCGTTCGTCCAGCTAAAAGGCCCGTTTTGCTACCA‑GGTGTTTTTTTCGTATGTCTGTCACCACCCTCGAGTCATAAATTGCGCAACGGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTc                                                                                                   <  1:3341430/147‑1 (MQ=255)

                                                  aggtttgtttTTATCTTTCGGTCAGATAAACGGCCCGCTTTGCTACCTGGTGTTTGTTTTCG‑TTGTCTGTTCCCACCCTCGAGAAAGAACATGCGCACCGGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTc                                                                                                   <  1:1245597/148‑1 (MQ=255)

                                                     tttgctttgcTCTTTCGGTCAGATAAACTGCCCGATTTGCTACCTGGGTATTTTTTTCT‑TTTTCTGTTACCCCAATCGAGTCAGAATATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                                                                                                <  1:2526675/150‑1 (MQ=255)

                                                     tttGNTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACAGATTTAACGTCTTGAAGATCGTccgc                                                                                                >  1:2691732/1‑150 (MQ=255)

                                                     tttGCTTTTCTCGTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                                                                                                <  2:2635455/150‑1 (MQ=255)

                                                          ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCCGTTGCGCCGACACAAAACAGCCAACGGAGTTACCGTCTTGAAGATCGTCcgccgcgg                                                                                           >  1:3335405/1‑148 (MQ=255)

                                                          ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt                                                                                           >  1:955179/1‑150 (MQ=255)

                                                          ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt                                                                                           >  1:2325664/1‑150 (MQ=255)

                                                          ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAACGTCTTGAAGATCGTCCGCCGCtt                                                                                           >  1:3402485/1‑150 (MQ=255)

                                                            ttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTcc                                                                                         <  2:1927224/150‑1 (MQ=255)

                                                                 tttCGTTCAGATAAATTGCCCGCTTTGCTACCAGGGGTTGGTTTTCG‑TTTTCTGTTTCCACCCTCTAGTCACAAATTGCGCAACGGTTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                    <  1:1970052/150‑1 (MQ=255)

                                                                 tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                    >  2:3911915/1‑150 (MQ=255)

                                                                 tttCGGTCAGCTAAAATGCCCGCTTGGCTACCAGTGTTTTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGACCTAAAATGCGAAACGGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                    <  1:3320147/150‑1 (MQ=255)

                                                                  ttCGGTCAGCTAAAAGGCCCGTTTTGCTACCAGGGTTTTTTTTTCT‑TTTTGTGTTTCCACCCTCGAGTCATAAAAGGCGCACCGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc                                                                                   <  1:3033943/150‑1 (MQ=255)

                                                                      gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACTGATTTAACGTCTTGAAGATAGTCCGCCGCGTCCTCAAAcccag                                                                               >  1:3662511/1‑148 (MQ=255)

                                                                       tCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGt                                                                              >  1:222342/1‑150 (MQ=255)

                                                                        cAGCTAAACTGCCCGCTTTGCTACCTCGGGATTGTTTTCT‑TTTTCTGTTTCCCCCCGCGAGTCAGAAATTGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGtt                                                                             <  1:3709459/150‑1 (MQ=255)

                                                                        cAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGtt                                                                             <  2:2462663/150‑1 (MQ=255)

                                                                         aGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCCGACTCAAAACAGCGAACGGAATTACCGTCGTGAAGATAGTCCGCCCCGTCCGCAAACCCAGttt                                                                            >  1:377033/1‑150 (MQ=255)

                                                                           cTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgc                                                                          >  2:2563498/1‑150 (MQ=255)

                                                                              aaaTGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCTGTTTgccgc                                                                       >  1:2894643/1‑150 (MQ=255)

                                                                                ttGCCCGCTTTGCTACCAGGGTTTTGTTTTCG‑TTGTCTTTTTCCACCCTCTAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCt                                                                     <  1:1299665/149‑1 (MQ=255)

                                                                                  gCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGt                                                                   <  2:1169697/150‑1 (MQ=255)

                                                                                  gCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGt                                                                   <  2:2493201/150‑1 (MQ=255)

                                                                                        tttGCTACCACGGGTTTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCAt                                                             <  2:3384223/150‑1 (MQ=255)

                                                                                        tgtGCTACCACGGGATTGTTTTCG‑CTGTCTGTTTCCACCCTCGAGTCAGAAACTGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGAGCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCAt                                                             <  2:3501709/148‑1 (MQ=255)

                                                                                         ttGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAATCAAAACAACGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAAACAGTATGCCGCCTGTCAGCATc                                                            >  2:1486874/1‑150 (MQ=255)

                                                                                               ccccGGGATTGTTTTCG‑ATGTCTGTTTCCACCCGCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAACCCAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTgg                                                      <  2:2687804/147‑1 (MQ=255)

                                                                                                  cGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCCGTTGCGCAGACACAAAACAGCCAACAGATTTACCCGCTTGAAGATCCTCCCCCGCATCCACAAACCCGGGTTGCCGCCAGTCCGCAACGGCTGGCgc                                                   >  1:3137009/1‑149 (MQ=255)

                                                                                                  cGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCACAAACCCTGTTTGCCGCCTGTCAGCATCAGCTGGCGa                                                   >  2:3010572/1‑150 (MQ=255)

                                                                                                     gATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGt                                                >  2:1119019/1‑150 (MQ=255)

                                                                                                      aTTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGtt                                               <  2:1812738/150‑1 (MQ=255)

                                                                                                       ttGTTTTCG‑TTGTCTGTTTCCCCCCTCGAGTCAGAAATTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTc                                              <  1:94465/150‑1 (MQ=255)

                                                                                                       ttGTTTTCG‑TTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTc                                              <  1:596347/150‑1 (MQ=255)

                                                                                                        tGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCg                                             >  2:152247/1‑150 (MQ=255)

                                                                                                         gTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGt                                            <  2:2325664/150‑1 (MQ=255)

                                                                                                          ttttCG‑TTTTCTGTTTCCACCCTCGAGTCATAAAATGCGCAACCGTTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                           <  1:1131392/150‑1 (MQ=255)

                                                                                                          ttttCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                           <  1:3589348/150‑1 (MQ=255)

                                                                                                          ttttCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                           >  2:356390/1‑150 (MQ=255)

                                                                                                          ttttCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                           <  1:288368/150‑1 (MQ=255)

                                                                                                          gtgtCG‑AGGTCTGTTGCCCCCCGCGCGCCCGAAACGGCGCACCCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGtt                                           <  1:288572/147‑1 (MQ=255)

                                                                                                           tttCG‑TTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTg                                          <  1:289968/150‑1 (MQ=255)

                                                                                                           tttCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACACGGTTGCGCAGACACAAAAAAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCACAAAACCAGTTTGACGCCAGTCAGCAACTGCTGGCGATGTTAGTTg                                          >  1:2120406/1‑150 (MQ=255)

                                                                                                            ttCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTTCGCAGACTCAAAAAAGCGAACAGATTTACCCTCTTTAAGATAGTCCGCCGCGTCCTCAAACCCAGTTTTCCGCCTGTCAGCAACTGCTGGCGATGTTGGTTGc                                         >  1:2044014/1‑150 (MQ=255)

                                                                                                             tCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCg                                        <  2:2027830/150‑1 (MQ=255)

                                                                                                              cG‑ATGTGTGTTTCCACCCTCGAGTCAGAAAATGCGCAACGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                       <  1:3677597/150‑1 (MQ=255)

                                                                                                                 aTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTg                                     >  2:3105360/1‑150 (MQ=255)

                                                                                                                  tttCTGTTTCCACCCTCTAGTCAGAAATTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGa                                    <  1:2832603/148‑1 (MQ=255)

                                                                                                                  tGTCTGTTTCCACCATCGAGTAATAAAATGCGAAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGa                                    <  2:2154811/150‑1 (MQ=255)

                                                                                                                   gTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCGGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGAt                                   >  2:654173/1‑150 (MQ=255)

                                                                                                                    tCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATc                                  >  2:638061/1‑150 (MQ=255)

                                                                                                                    tCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATc                                  >  2:1616188/1‑150 (MQ=255)

                                                                                                                    tCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACAAAACAGCGAACAGATTTACCGGCATTAAGATCGTCCGCCGCGTCCGCAAAACCTGGTTTGCGCCTGTCAGCCTATGCTGGCGAAGTTAGGTGCGGTGATc                                  >  1:2237271/1‑150 (MQ=255)

                                                                                                                     cTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCa                                 >  2:2292541/1‑150 (MQ=255)

                                                                                                                      tGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGACGATCGTCCGCCGCTTCCTCAAACCCTGGTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCAc                                >  2:606952/1‑150 (MQ=255)

                                                                                                                      tGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCAc                                >  2:2794178/1‑150 (MQ=255)

                                                                                                                        tttCCCCCCTCGAGTCAGAAAAGGCGCAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGc                              <  1:3028567/150‑1 (MQ=255)

                                                                                                                        tttCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTAACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCGGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGc                              >  2:1531421/1‑150 (MQ=255)

                                                                                                                                     gTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAAt                 <  2:1230653/150‑1 (MQ=255)

                                                                                                                                        tgAATATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAACGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTg              <  2:1935940/149‑1 (MQ=255)

                                                                                                                                         gAATATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGc             <  1:1890028/150‑1 (MQ=255)

                                                                                                                                         gAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGc             <  1:3911915/150‑1 (MQ=255)

                                                                                                                                          aaaaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGcc            <  2:773657/150‑1 (MQ=255)

                                                                                                                                            aaTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAg          >  2:2079835/1‑150 (MQ=255)

                                                                                                                                                  cAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg    <  2:877951/150‑1 (MQ=255)

                                                                                                                                                   aaCGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   <  1:2944804/150‑1 (MQ=255)

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

ATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCG‑ATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828833‑3829129

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