BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,618,947	A→T	5.6%	I405N (ATC→AAC) ‡	uraA ←	uracil permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,618,947	0	A	T	5.6%	635.8 / 19.5	178	I405N (ATC→AAC) ‡	uraA	uracil permease
Reads supporting (aligned to +/- strand): ref base A (75/93); new base T (2/8); total (77/101)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 185 total aligned reads displayed.

ACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTTTCGAT  >  NC_000913/2618798‑2619096

acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                       >  1:583867/1‑150 (MQ=255)

acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                       >  2:2912400/1‑150 (MQ=255)

acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                       >  1:2037588/1‑150 (MQ=255)

 caAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGAt                                                                                                                                                      <  2:3740463/150‑1 (MQ=255)

   aaaTCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCa                                                                                                                                                    <  2:1762176/150‑1 (MQ=255)

     aTCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTAGATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCaaa                                                                                                                                                  <  2:1971233/150‑1 (MQ=255)

      tCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAAc                                                                                                                                                 <  1:1347661/150‑1 (MQ=255)

          gagaCGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                             <  2:4129815/150‑1 (MQ=255)

          gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATCTCCGCGTCTTCTGCGTCCAGCACCTCTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                             <  2:4129849/150‑1 (MQ=255)

          gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGATCTTCAAACtgat                                                                                                                                             <  1:3785413/150‑4 (MQ=255)

          gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                             <  1:2781960/150‑1 (MQ=255)

          gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTAGATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGATCTTCAAACtgat                                                                                                                                             <  1:1578509/150‑4 (MQ=255)

          gagaAGATTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                             <  1:2883691/150‑1 (MQ=255)

              aGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCg                                                                                                                                         >  1:2138720/1‑150 (MQ=255)

               gTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGa                                                                                                                                        <  2:1366049/150‑1 (MQ=255)

                tttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGAt                                                                                                                                       >  1:2502232/1‑150 (MQ=255)

                tttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGAt                                                                                                                                       >  1:1006623/1‑150 (MQ=255)

                 ttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATa                                                                                                                                      <  2:2592377/150‑1 (MQ=255)

                 ttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATa                                                                                                                                      <  1:3496340/150‑1 (MQ=255)

                     cacaGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcga                                                                                                                                  <  2:1066054/150‑1 (MQ=255)

                     cacaGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcga                                                                                                                                  <  1:2078595/150‑1 (MQ=255)

                     cacaGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcga                                                                                                                                  <  2:2646748/150‑1 (MQ=255)

                       caGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacg                                                                                                                                >  2:1222232/1‑150 (MQ=255)

                        aGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacga                                                                                                                               <  2:1015500/150‑1 (MQ=255)

                               cATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGc                                                                                                                        >  1:2197550/1‑150 (MQ=255)

                                  gatgTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTTTCTGTTATATAAGCGTCTTATGCGTACAGCACCACATCTTAAGGGCGCAGCACCAAGATAATTTTGAAGATAAAAATTAAGCCGATACCGAAAATGTTCGCCac                                                                                                                     >  1:354020/1‑149 (MQ=255)

                                  gatgTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAg                                                                                                                     >  1:2981625/1‑150 (MQ=255)

                                     gTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAgcgc                                                                                                                  >  1:2153312/1‑150 (MQ=255)

                                       ttGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCa                                                                                                                <  2:3321115/150‑1 (MQ=255)

                                       ttGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCa                                                                                                                <  1:3691841/150‑1 (MQ=255)

                                         gAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATa                                                                                                              <  2:3595799/150‑1 (MQ=255)

                                               ggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCtttt                                                                                                        <  1:1025133/150‑1 (MQ=255)

                                               ggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCtttt                                                                                                        >  1:3872975/1‑150 (MQ=255)

                                               ggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCtttt                                                                                                        <  2:703318/150‑1 (MQ=255)

                                                gcagcaACACTGCCCGGTCGGGACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTa                                                                                                       >  1:3584189/1‑150 (MQ=255)

                                                gcagcaACACTGCCAGGTCGGTACCATATTTGTGTGTTATATCCGCGTATTCAGCGTCCAGCACCACTTCTTCAGGGCGCCGCACACTGATCAGTTTGAAGATAAAACTTAAGCCGATACCCACGATGGTAGCCAGCGCCAAACCTtata                                                                                                       >  1:3584225/1‑150 (MQ=255)

                                                           gCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCcgccgcg                                                                                            <  1:3522366/150‑1 (MQ=255)

                                                            cccGGTCGGTACATTATTTTTCTGTTATATCCGCGTCTTCTTCGTCCAGCACCACATCTTAAGGGCGCAGCACACAGATCAGTTTGAATATAAAACATAAGCCGATACCGACGAATGTCGCCACCGCCAAACCAATTAACTAccccccgc                                                                                           >  1:2452282/1‑150 (MQ=255)

                                                            cccGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCcgccgcgc                                                                                           <  2:2697723/150‑1 (MQ=255)

                                                                    gtACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGttt                                                                                   >  1:3975829/1‑150 (MQ=255)

                                                                      aCATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTAc                                                                                 <  2:1246957/150‑1 (MQ=255)

                                                                          tattTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTc                                                                             >  1:1912946/1‑150 (MQ=255)

                                                                          tattTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTc                                                                             <  2:3707372/150‑1 (MQ=255)

                                                                           attTGTCTGTTATCTCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTcg                                                                            <  2:180460/150‑1 (MQ=255)

                                                                            tttGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTcgc                                                                           <  1:453792/150‑1 (MQ=255)

                                                                            cttGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTcgc                                                                           <  1:95381/149‑1 (MQ=255)

                                                                              tGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCgcgc                                                                         <  1:158594/150‑1 (MQ=255)

                                                                              tGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCACCGCGCCGATGTTTACCTTCgcgc                                                                         <  1:814412/150‑1 (MQ=255)

                                                                              tGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTGCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCgcgc                                                                         <  2:2853813/150‑1 (MQ=255)

                                                                                tCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCa                                                                       <  1:3344394/150‑1 (MQ=255)

                                                                                 cTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCCGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCAc                                                                      <  1:1133470/150‑1 (MQ=255)

                                                                                  tGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACt                                                                     >  2:2505394/1‑150 (MQ=255)

                                                                                    ttATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGa                                                                   >  2:1463810/1‑150 (MQ=255)

                                                                                    ttATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGa                                                                   <  2:3406219/150‑1 (MQ=255)

                                                                                     tataTCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGAc                                                                  <  2:2731129/150‑1 (MQ=255)

                                                                                        atCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGcc                                                               <  2:1850287/150‑1 (MQ=255)

                                                                                               tcttctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATc                                                        <  2:3850278/150‑1 (MQ=255)

                                                                                               tcttctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATAAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATc                                                        <  2:2439147/150‑1 (MQ=255)

                                                                                                cttctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCa                                                       >  1:930698/1‑150 (MQ=255)

                                                                                                 ttctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGAGGTTTACCTTCGCGCCACTGACGCCGATGATCaa                                                      >  2:1908626/1‑150 (MQ=255)

                                                                                                   ctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCaaaa                                                    >  1:2741053/1‑150 (MQ=255)

                                                                                                    tGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAAt                                                   <  2:3062822/150‑1 (MQ=255)

                                                                                                      cGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCa                                                 <  2:964276/150‑1 (MQ=255)

                                                                                                         ccAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACgg                                              >  1:2760527/1‑150 (MQ=255)

                                                                                                           aGCACCACTTCTTCAGGGCGCAGCACTCTGATCTGTTTGTAGATCAAAGTTATGCCGTTAGCGACGATGGTCGCCAGCGCCATAGCTTTTATCTCCGCCGCGCCGTTGTTTTCCTTCGCGCCTCTGACGGCGATGTTCTTATTCTCGGaa                                            <  2:86254/150‑1 (MQ=255)

                                                                                                           aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGaa                                            >  2:1028396/1‑150 (MQ=255)

                                                                                                           aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGaa                                            >  1:2574449/1‑150 (MQ=255)

                                                                                                           aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGaa                                            >  1:2688114/1‑150 (MQ=255)

                                                                                                              accacTTCTTCCGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTc                                         <  1:2443773/150‑1 (MQ=255)

                                                                                                              accacTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTc                                         <  2:2196655/150‑1 (MQ=255)

                                                                                                               ccacTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCa                                        <  2:2385587/150‑1 (MQ=255)

                                                                                                               ccacTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAAACACGGAAGTCa                                        >  2:3101924/1‑150 (MQ=255)

                                                                                                                cacTTCTTCAGGGCGCAACACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAg                                       >  2:391781/1‑150 (MQ=255)

                                                                                                                 acTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAgg                                      >  1:1207140/1‑150 (MQ=255)

                                                                                                                      ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACTATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                                 <  2:16894/150‑1 (MQ=255)

                                                                                                                      ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCAGGGAAGTCAGGATCAg                                 >  2:730719/1‑150 (MQ=255)

                                                                                                                      ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                                 >  1:3528450/1‑150 (MQ=255)

                                                                                                                      ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                                 >  1:2512902/1‑150 (MQ=255)

                                                                                                                       tcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAgg                                >  2:3857133/1‑150 (MQ=255)

                                                                                                                           ggCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCt                            <  2:1208388/150‑1 (MQ=255)

                                                                                                                             cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                          <  2:152869/150‑1 (MQ=255)

                                                                                                                             cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                          <  1:1694359/150‑1 (MQ=255)

                                                                                                                             cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                          <  1:1694634/150‑1 (MQ=255)

                                                                                                                             cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                          <  1:2369493/150‑1 (MQ=255)

                                                                                                                              gcagcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGAGGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgtg                         <  2:3284587/150‑1 (MQ=255)

                                                                                                                                 gcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCtt                      <  1:3967786/150‑1 (MQ=255)

                                                                                                                                 gcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCtt                      >  1:1251180/1‑150 (MQ=255)

                                                                                                                                  cacaCTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCttt                     <  1:1501792/150‑1 (MQ=255)

                                                                                                                                     acTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTttgtt                  >  2:1297960/1‑150 (MQ=255)

                                                                                                                                     acTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTttgtt                  <  2:583867/150‑1 (MQ=255)

                                                                                                                                     acTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTttgtt                  >  2:217563/1‑150 (MQ=255)

                                                                                                                                        cttaaGTTTGATCTTCAAACTTAAGCCGATACCGACGATGGTCGCAAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTa               >  2:4043879/5‑150 (MQ=255)

                                                                                                                                         aTCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTaa              >  1:3182341/1‑150 (MQ=255)

                                                                                                                                          tCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAAt             <  2:3404868/150‑1 (MQ=255)

                                                                                                                                           cAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATc            <  1:90942/150‑1 (MQ=255)

                                                                                                                                                tgtaGATCAAAGTTAAGCCGATACCGACGTTGGTCGCCTGCGCCATACGTTTTAACTGGGCCGCGCCGATGTTTACCTTCGCGCGACTGACGCCGATGATCAAATTCTGGGTTGTCAGGATCTGGTTCTGTGCTTTGTTGTATTCCAGtt       >  2:3097507/4‑150 (MQ=255)

                                                                                                                                                tGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACtt       >  2:2127383/1‑150 (MQ=255)

                                                                                                                                                tGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTATGTTGTAATCCACtt       >  2:1427117/1‑150 (MQ=255)

                                                                                                                                                  aaGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTttc     <  2:3975829/150‑1 (MQ=255)

                                                                                                                                                    gATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTttcga   >  2:2274775/1‑150 (MQ=255)

ACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTTTCGAT  >  NC_000913/2618798‑2619096

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