BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,828,980	0	A	C	8.1%	897.5 / 11.2	359	Q12H (CAA→CAC)	xanP	xanthine permease
Reads supporting (aligned to +/- strand): ref base A (173/147); new base C (6/32); total (180/179)
Fisher's exact test for biased strand distribution p-value = 8.69e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 370 total aligned reads displayed.

ATATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGG  >  NC_000913/3828831‑3829128

atataAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCaa                                                                                                                                                      <  2:4473764/150‑1 (MQ=255)

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

   taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAACTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCg                                                                                                                                                   <  2:2420222/150‑1 (MQ=255)

    aaaaaaGTGCAATTATCACAATTAAATATTTAGATGAGGATGTGGCAAACGATAGCTATTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCACGAGTCAGAAAATGCGCAACCgg                                                                                                                                                  <  2:6991334/150‑1 (MQ=255)

     aaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGt                                                                                                                                                 >  1:6273909/1‑150 (MQ=255)

         aGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgcg                                                                                                                                             <  2:5011820/150‑1 (MQ=255)

         aGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgcg                                                                                                                                             <  2:3952220/150‑1 (MQ=255)

         aGTGCATTTATCTCATTTAAAGATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAACATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgcg                                                                                                                                             <  2:2879477/150‑1 (MQ=255)

          gTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgcgc                                                                                                                                            <  2:5072219/150‑1 (MQ=255)

           tGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCa                                                                                                                                           >  2:2317679/1‑150 (MQ=255)

                ttATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTc                                                                                                                                      >  2:1133669/1‑150 (MQ=255)

                  aTCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCaa                                                                                                                                    <  2:1582038/150‑1 (MQ=255)

                  aTCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTTCGCAGACTCaa                                                                                                                                    >  1:3315583/1‑150 (MQ=255)

                    ctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaaa                                                                                                                                  >  1:4632705/1‑150 (MQ=255)

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

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

                       aTTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGGTGCGCAGACTCAAAACAg                                                                                                                               >  2:5473770/1‑150 (MQ=255)

                       aTTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCAACCCACGAGTCAGAAAATGCGCAACAGGTTGCGCAGAAACAAAAAAg                                                                                                                               >  2:6322398/1‑150 (MQ=255)

                       aTTTAAATATTAAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAg                                                                                                                               >  1:4417763/1‑150 (MQ=255)

                        tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                                                                                              <  2:2549449/150‑1 (MQ=255)

                           aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGaa                                                                                                                           >  2:1910128/1‑150 (MQ=255)

                            aaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACAAAACAGCGAAc                                                                                                                          >  1:702828/1‑150 (MQ=255)

                              aatTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTATTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGACTGTTTCCACCCTCGACTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTg                                                                                                                        <  2:3565407/149‑1 (MQ=255)

                               atTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTTCGCAGACTCAAAAAATCGAACAGa                                                                                                                       >  1:5482728/1‑150 (MQ=255)

                                 ttAGATGAGGATGTGGCAAACGTTTGCTTTTCACTTTCGGTCAGCTAAACTGCCCGCTTTGCTACCACGGGATGGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAtt                                                                                                                     <  2:489117/150‑1 (MQ=255)

                                  tAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAttt                                                                                                                    <  2:5318009/150‑1 (MQ=255)

                                   aGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGCGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGATTTa                                                                                                                   >  1:538953/1‑150 (MQ=255)

                                     aTGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAcc                                                                                                                 <  2:1328663/150‑1 (MQ=255)

                                        tggttgtgCCAAACTTGTGATTTTCTCTTGCGGACAGAAAAAAGGCCCGCTTTGCTACCACGGTTTTGTTTTCGTTGTCTGTTTCCACCCTCGAGTCATAAATTGCGCACCGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTc                                                                                                              <  1:4113553/146‑1 (MQ=255)

                                        aGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCTACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTc                                                                                                              <  2:111564/150‑1 (MQ=255)

                                         ggATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCt                                                                                                             <  2:1446311/150‑1 (MQ=255)

                                          gATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACAAAACAGCGAACTGATTTACCGTCtt                                                                                                            >  2:3387737/1‑150 (MQ=255)

                                             gtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACACAAAAAAACGAACAGATTTAACGTCTTGaa                                                                                                         >  1:1588594/1‑150 (MQ=255)

                                                 caaagttttGCTTTTCTCTTTCGGTCATATAAATTGCCCGTTTTGCTACCAGGGTTTTTTTTTCGTTGTCTGTTCCCACCCTCGAGTCATAAAATGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                     <  1:3919860/144‑1 (MQ=255)

                                                     cGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTCCTACCCCGGTTTTTTTTTCTTTGTCTGTTTCCACCCTCTAGTAATAAAATGCGCAACCGGTTGCCCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTcc                                                                                                 <  1:5243307/150‑1 (MQ=255)

                                                        ttGCTTTTCTCTTTCTGTCAGCTAAATTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGAAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgcc                                                                                              <  2:837663/150‑1 (MQ=255)

                                                            ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCcgccgcgt                                                                                          >  1:2002427/1‑148 (MQ=255)

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

                                                            ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAAATGATTTAACGTCTTGAAGATAGTCCGCCGCtt                                                                                          >  1:4993068/1‑150 (MQ=255)

                                                             tttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTc                                                                                         <  2:2179449/150‑1 (MQ=255)

                                                             ttgcccgttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCGCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTc                                                                                         <  2:5703986/143‑1 (MQ=255)

                                                              gtctcTGTCGGTCAGCTAACAGGCCCGCTGTGCTACCACGGGATTGTTTGCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAACACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTcc                                                                                        <  2:4638813/149‑1 (MQ=255)

                                                                 tctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTTCGCAGACACAAAACAGCGAACTGAGTTAACCTCTTGAAGATAGTCCGCCGCTTCCTCa                                                                                     >  1:2680623/1‑150 (MQ=255)

                                                                   tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                   <  2:4084618/150‑1 (MQ=255)

                                                                   tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCaaa                                                                                   >  1:2295102/1‑150 (MQ=255)

                                                                   tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGGCCGCCGCTTCCGCaaa                                                                                   >  1:3410901/1‑150 (MQ=255)

                                                                   ttgcGGCCAGCTAAAAGGCCCGCTTTGCTACCAGGGTTTTTTTTTCGTTGTCTGTTTCCACCCTCGAGACATAAATTGCCCACCGTGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                   <  1:6222556/147‑1 (MQ=255)

                                                                    ttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc                                                                                  >  2:2653955/1‑150 (MQ=255)

                                                                     tCGGTCAGCTAAATTGCCCGCTTTGCTACCTGGGTTTTTTTTTCGATGTGTGTTTCCACCCTCGAGTCATAAATTGCGCACCGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAcc                                                                                 <  1:5741561/150‑1 (MQ=255)

                                                                           aGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCACAAACCCTGttt                                                                           >  1:375577/1‑150 (MQ=255)

                                                                             cTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgc                                                                         <  2:7207289/150‑1 (MQ=255)

                                                                             cTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAACGTCTTGAAGATCGTCCGCCGCTTCCGCAAACCCAGttttc                                                                         >  1:1874809/1‑148 (MQ=255)

                                                                              tAAAAGGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgcc                                                                        <  2:1002952/150‑1 (MQ=255)

                                                                                  aTGCCCGCTTTGCTACCACGGTATTGTTTTCTTTGTCTTTTTCCACCCTCGAGTCAGAAATGGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCt                                                                    <  1:3596602/150‑1 (MQ=255)

                                                                                      ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGATTTAACGTCTTTAAGATCGTCCGCCGCTTCCACAAACCCGGTTTGCCGCCTGTCa                                                                >  1:2983919/1‑150 (MQ=255)

                                                                                       cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCAGTTTGCCGCCAGTCAg                                                               >  1:81378/1‑150 (MQ=255)

                                                                                         ctgtGCTACCAGGGGATTGTTTTCGATTTCTGTTACCACCCTATAGAAAGAATATGCTAAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCTTGTTGAAGATCGTCCGCCGCTTCCTCATACCCTGTTTGCAGCCGGTCAGCa                                                             <  2:313154/147‑1 (MQ=255)

                                                                                          tttGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCAt                                                            <  1:3453770/150‑1 (MQ=255)

                                                                                           ttGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATc                                                           <  2:6788338/150‑1 (MQ=255)

                                                                                             gCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACAAAAAAACGAAAAGATTTAACGTCATTAAGATCGTCCGCCCCATCCTCAAAACCAGTTTTCCGCCGGTCAGCaacag                                                         >  1:6821582/1‑146 (MQ=255)

                                                                                              cTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATAGTCCGCCGCGTCCACAAACCCTGTTTGCCGCCTGTCAGCATctgc                                                        >  1:5937748/1‑150 (MQ=255)

                                                                                               tACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgct                                                       >  2:1885971/1‑150 (MQ=255)

                                                                                                 ccccGGTATTGTTTTCGCTGTCTGTTTCCACCCTCGAGACACAAATTGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTgg                                                     <  1:7208407/147‑1 (MQ=255)

                                                                                                 ccACGGGTTTGTTTTCGTTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTgg                                                     <  1:1304907/150‑1 (MQ=255)

                                                                                                 ccACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAACGAACTGATTTAACGTCATGAAGATCGTCCGCCGCTTACTCAAAACCAGTTTGCCGCCAGTCAGCATAAGCTgg                                                     >  1:4783374/1‑150 (MQ=255)

                                                                                                 ccACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAAGGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTgg                                                     <  1:170131/150‑1 (MQ=255)

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

                                                                                                     gggATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGAt                                                 >  2:389510/1‑150 (MQ=255)

                                                                                                        aTTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGCGTGCGCAGACACAAAACACCGAACAGCGTTACCGTCTTGAAGATCGGCCCCCGCGGCCCCAAACCCGGTTTGCCGCCGGTCACCAGCTGCTGGCGATGtg                                              >  1:2786504/1‑149 (MQ=255)

                                                                                                            ttttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCCGGTGCGCAGACTCAAAACACCGAACAGATTTAACGACATGAAGAACGTCCGCCCCATCCGCAAACCCAGTTTGCCGCCCGTCAGCAACTGCAGGCGCTGTTCGtt                                          >  1:5130511/1‑150 (MQ=255)

                                                                                                            ttttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTTCGCAGACACAAAACAGCGAACTGATTTAACCTCTTTAAGATAGTACCCCCCATCCACAAAACCAGTTTTGCGCCGGTAAGCATAAGCTGGCGATGTTAGtt                                          >  1:4588636/1‑150 (MQ=255)

                                                                                                             tttCGTGGTCTGTTACCACCCTCGAGGCACAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTg                                         <  1:2725703/150‑1 (MQ=255)

                                                                                                             tttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTg                                         <  2:2800370/150‑1 (MQ=255)

                                                                                                              ttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCACAAACCCGGTTTTCCGCCGGTCAGCATCTGCTGGCGATGTTAGTTGc                                        >  1:1096911/1‑150 (MQ=255)

                                                                                                               tCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGGGCAGACTCAAAAAACCGAACAGATTTAACGTCGTGAAGATAGTCCGCCCCGTCCTCAAACCCAGTTTTCCCCCGGTCAGCAACTGCTGGCGAAGTTCGTTggg                                       >  1:1466719/1‑148 (MQ=255)

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

                                                                                                                    gTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGAt                                  <  2:3027707/150‑1 (MQ=255)

                                                                                                                    gTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGAt                                  <  2:4417706/150‑1 (MQ=255)

                                                                                                                       tttttCCACCCTCGATTATTATATGGCGCAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCAc                               <  1:1421105/148‑1 (MQ=255)

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

                                                                                                                       tGTTTCCACCCTCGAGCCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCAc                               <  2:1385899/150‑1 (MQ=255)

                                                                                                                             cACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCGGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgc                         >  2:289095/1‑150 (MQ=255)

                                                                                                                              aCCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcg                        <  1:3635616/150‑1 (MQ=255)

                                                                                                                                ccTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCt                      <  2:5626258/150‑1 (MQ=255)

                                                                                                                                ccTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGAGGTTCGTTGCGGTGATCACGCCAGCGCt                      >  2:6199469/1‑150 (MQ=255)

                                                                                                                                 cTCGAGTCAGAAAATGCGCACCGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTa                     <  1:4336017/150‑1 (MQ=255)

                                                                                                                                 cTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTa                     >  2:7035663/1‑150 (MQ=255)

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

                                                                                                                                   cGAGTCCGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTAtt                   <  1:2130761/150‑1 (MQ=255)

                                                                                                                                    gAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTa                  <  1:4509882/150‑1 (MQ=255)

                                                                                                                                      gacAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAAt                <  1:1012459/148‑1 (MQ=255)

                                                                                                                                      gTCAGAAATTGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAAt                <  1:4255691/150‑1 (MQ=255)

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

                                                                                                                                       tCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGGTCGTTGCGGTGATCACGCCAGCGCTATTAATc               >  2:5167585/1‑150 (MQ=255)

                                                                                                                                         aGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTg             >  2:5786162/1‑150 (MQ=255)

                                                                                                                                                  gcACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCt    <  1:1568326/150‑1 (MQ=255)

                                                                                                                                                  gcAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTGCCGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCt    >  1:91433/1‑150 (MQ=255)

                                                                                                                                                  gaaaCGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCt    <  1:610030/148‑1 (MQ=255)

                                                                                                                                                   cAACCGTTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg   <  1:6450815/150‑1 (MQ=255)

                                                                                                                                                    aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg  >  2:6389474/1‑150 (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: `‑`