BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,068,557	T→A	7.0%	K428M (AAG→ATG)	rutG ←	pyrimidine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,068,557	0	T	A	7.0%	1138.1 / 11.4	355	K428M (AAG→ATG)	rutG	pyrimidine permease
Reads supporting (aligned to +/- strand): ref base T (147/171); new base A (31/3); total (181/174)
Fisher's exact test for biased strand distribution p-value = 2.46e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 365 total aligned reads displayed.

ACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAA  >  NC_000913/1068410‑1068706

  tGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCt                                                                                                                                                   >  2:6604874/1‑150 (MQ=255)

  tGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCt                                                                                                                                                   >  2:4029817/1‑150 (MQ=255)

  tGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCt                                                                                                                                                   >  1:5699179/1‑150 (MQ=255)

  tGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCGGCCCCGCCGGTCACCCCGGCGGTTTTCTGTGTATGGCGGCTGATGTGCACCGTCTGGCGGGGGAACGTCAACCAACTTTCt                                                                                                                                                   >  2:2524667/1‑150 (MQ=25)

       tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTa                                                                                                                                              >  2:1536088/1‑150 (MQ=255)

        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTaa                                                                                                                                             >  2:3042841/1‑150 (MQ=255)

        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACATTCTGCTTaa                                                                                                                                             >  2:3290360/1‑150 (MQ=255)

            gCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGGTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACATTATGCTTaacaac                                                                                                                                         >  1:1045632/1‑150 (MQ=255)

             ccGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACg                                                                                                                                        >  2:487061/1‑150 (MQ=255)

              ctgATAAGTCGTTTACGCCGCATCCGGCATGAAAAATGAGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTATGTTTATGTCTCCTGATTAACAACTTATGGCGGTGGAACGTGTACCAACTTTCTGCTTaaaaaagg                                                                                                                                       >  2:3888564/3‑147 (MQ=255)

              cGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGc                                                                                                                                       <  1:4934301/150‑1 (MQ=255)

                   aaGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGa                                                                                                                                  <  2:5231535/150‑1 (MQ=255)

                   aaGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGa                                                                                                                                  <  2:2907621/150‑1 (MQ=255)

                        gTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTGTGTTTATTGCTGCAGAAGAACAAATTCATGCGGTGGAAAGTCAACCAAAATTCTGCTTAACAACGCAATGATTaaa                                                                                                                             >  1:1609062/1‑150 (MQ=37)

                         tttACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTGTTCTGTTTATTGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACATTATGCTTAAAAACGCATTGAGTaaaa                                                                                                                            >  1:4227304/1‑150 (MQ=255)

                           tACGCAGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACAGCCGGTAACCCCGGCGGTTTTCTGTTTAAGGCTACTGATGAACAACATCTGGCGGTGGAACGTCAACCAACAATCTGCATAAAAACGCAATGAGTAAAAta                                                                                                                          >  2:6460488/1‑149 (MQ=255)

                              gccgcATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACGTCTGGCGGTGGAACGGCAACCAACTTTCTGCTTAAAAACGCATTGAGTAAAATcgcg                                                                                                                       >  1:3635714/1‑150 (MQ=255)

                                cgcATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATTGCTCCTGATGAACAACTTCAGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAAAGCATTTAGTAAAATCGCGcc                                                                                                                     >  2:2581065/1‑150 (MQ=255)

                                      cGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTGTGTTTATTGGTCCTGATTTACAACATATTGCGGTGGGAAGTCAACCAACATTGTTCTTTAAAAAGCAATTAGTAAAATAGCGCCCAAGGt                                                                                                               >  1:5151281/1‑150 (MQ=255)

                                      cGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTGTGTTTATTGCTCCAGATGTACAACTTATGGCGGTGGGACGTCAAACAACTTTCTGCTTAACAAAGCATTGATGAAAATAGCGCCAAAGGt                                                                                                               >  1:2644317/1‑150 (MQ=255)

                                       ggCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTGTTCGGTTTATGGCGCCGGGTGAACAACAGCAGGCGGTGGCACCGCCACCAACCTTCAGCGTAACAACGCATTGAGTAAAATCGCGCCACAgggt                                                                                                              >  1:3613470/1‑148 (MQ=255)

                                        gCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTTTTTTTGGCTCCAGATGAACAAATTCTGGCGGTGGAACGGAAACCCACATAATGCTTAAAAACGCAGTGAGTTAAAAAGCGCCAAAGGTTg                                                                                                             >  1:995041/1‑150 (MQ=255)

                                        gCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTg                                                                                                             <  1:1410228/150‑1 (MQ=255)

                                          aTGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCt                                                                                                           <  1:711098/150‑1 (MQ=255)

                                                   gCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACGTCTGGCGGTGGAACGTCAAACAAATTTCTGCTTAACAACGCAATTAGTAAAATCGCGCCAAAGGTTGCTGTACCAAtc                                                                                                  >  2:3392631/1‑150 (MQ=255)

                                                     gTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCAGATGAACAAAATCTGGCGGTGGAACGTCAACCAACATTCTGCTTAAAAACGCCTTGAGTAAAATAGCGCCAAATGTTGCTGTACCAAaccc                                                                                                >  2:1429939/1‑150 (MQ=255)

                                                     gTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCAGTTTATGGATCCTTATGAACAACTTCTTGCGGTGGAACGTCAACCAACATTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAAtccc                                                                                                >  2:5386693/1‑150 (MQ=255)

                                                      tACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAAtccct                                                                                               >  2:1565212/1‑150 (MQ=255)

                                                      tACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTAATGGTTATGGCTCCTGATGGACAACATCTGGCGGTGGAACGGCAACCAACTTTCTGCTTAAAAACGCATTGAGGAAAAACGCCCCAAAGGTTGCCGGACCCAtccct                                                                                               >  2:5618231/1‑150 (MQ=255)

                                                        cGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAAtccctcc                                                                                             <  1:5651241/150‑1 (MQ=255)

                                                        cGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAAtccctcc                                                                                             <  1:2484405/150‑1 (MQ=255)

                                                        cGTTGTCAACAATCGGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAAtccctcc                                                                                             <  1:2548425/150‑1 (MQ=255)

                                                          ttGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCa                                                                                           <  2:2283783/150‑1 (MQ=255)

                                                           tGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCaa                                                                                          <  1:2140150/150‑1 (MQ=255)

                                                           tGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCaa                                                                                          <  1:1728284/150‑1 (MQ=255)

                                                           tGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCaa                                                                                          <  2:6188339/150‑1 (MQ=255)

                                                           tGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCaa                                                                                          >  2:2124020/1‑150 (MQ=255)

                                                              caacaaTCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGt                                                                                       <  1:1426480/150‑1 (MQ=255)

                                                                acaaTCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTaa                                                                                     <  2:1181906/150‑1 (MQ=255)

                                                                 caaTCTGCACCGCCGGTAACCCCGGCGGTTTTATGTTTATGGCTCCTGTTGAACAACTTCTGGCGGTGGAAAGTCAAACAACTTTATGCATAAAAACCCATTTAGTAAAATAGCGCCAAAGGTTGCTGTACCAATCCCACCCAACCTaaa                                                                                    >  2:1884371/1‑150 (MQ=255)

                                                                          ccgccgGTAACCCCGGCGGTTTTCTGTTTATGGCGCCGGGTGGACAACTGCTGGCGGTGGAACGGCAACCAACATTCTGCGTAAAAACCCCTGGAGTAAAATCGCGCCAAAGGTTGCGGTTCCAATCCCTCCCAACCGAAAAACGccccc                                                                           >  2:1703053/1‑148 (MQ=255)

                                                                          ccgccgGTAACCCCGGCGGTGTTCTGTTTATGGCGCCTGGTGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCAGTGAGTAAAATAGCGCCAAAGGTTGCTGTACCAATACCTCCCAACGGAAAAACGCCCAg                                                                           >  2:4593439/1‑150 (MQ=255)

                                                                              cgGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTg                                                                       <  2:6154845/150‑1 (MQ=255)

                                                                               ggTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTga                                                                      <  1:743913/150‑1 (MQ=255)

                                                                                gTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTgag                                                                     <  1:6505233/150‑1 (MQ=255)

                                                                                gTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTgag                                                                     <  1:2964730/150‑1 (MQ=255)

                                                                                 tAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTgaga                                                                    <  2:4305158/150‑1 (MQ=255)

                                                                                  aaCCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGagag                                                                   <  2:137069/150‑1 (MQ=255)

                                                                                   aCCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAACACCGCCCAGCGTGAGAGc                                                                  <  1:1565212/150‑1 (MQ=255)

                                                                                     cccGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACAAACTTTCTGCTTAACAACGCTTTGAGTAAATTCTCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGcga                                                                <  1:4808187/150‑3 (MQ=255)

                                                                                      ccGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAACACCGCCCAGCGTGAGAGCaaa                                                               <  2:5526945/150‑1 (MQ=255)

                                                                                       cggcggTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTACAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCaaaa                                                              <  2:6529711/150‑1 (MQ=255)

                                                                                       cggcggTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCaaaa                                                              >  2:5474844/1‑150 (MQ=255)

                                                                                       cggcggTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCaaaa                                                              <  1:4179422/150‑1 (MQ=255)

                                                                                           ggTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCAc                                                          <  2:262390/150‑1 (MQ=255)

                                                                                           ggTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCAc                                                          >  1:2682349/1‑150 (MQ=255)

                                                                                           ggTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCAc                                                          <  2:1988345/150‑1 (MQ=255)

                                                                                           ggTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCAc                                                          <  1:5203948/150‑1 (MQ=255)

                                                                                            gTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTATAACCGCCCAGCGTGAGAGCAAAATCAcc                                                         <  1:6179629/150‑1 (MQ=255)

                                                                                            gTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCAcc                                                         <  1:1248569/150‑1 (MQ=255)

                                                                                              tttCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCcg                                                       >  2:5881469/1‑150 (MQ=255)

                                                                                                 cTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCgcgc                                                    >  1:3942375/1‑150 (MQ=255)

                                                                                                   gTTTATGTCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGATTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGccc                                                  <  1:4565104/150‑1 (MQ=255)

                                                                                                   gTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGccc                                                  >  2:3579684/1‑150 (MQ=255)

                                                                                                    tttATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCa                                                 <  1:5994336/150‑1 (MQ=255)

                                                                                                     ttATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAg                                                <  2:6124073/150‑1 (MQ=255)

                                                                                                         ggCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACt                                            >  1:3421246/1‑150 (MQ=255)

                                                                                                         ggCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACt                                            >  1:5442641/1‑150 (MQ=255)

                                                                                                         ggCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACt                                            >  2:1077611/1‑150 (MQ=255)

                                                                                                           cTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTaa                                          <  1:5094092/150‑1 (MQ=255)

                                                                                                              cTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGt                                       >  2:3649046/1‑150 (MQ=255)

                                                                                                               tgatgaACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGtt                                      >  2:3110879/1‑150 (MQ=255)

                                                                                                                 atgaACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTAc                                    <  2:3677820/150‑1 (MQ=255)

                                                                                                                   gaACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTg                                  <  2:6274757/150‑1 (MQ=255)

                                                                                                                   gaACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTg                                  >  2:1647059/1‑150 (MQ=255)

                                                                                                                   gaACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTg                                  <  2:4690021/150‑1 (MQ=255)

                                                                                                                     acaacTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCg                                >  1:407471/1‑150 (MQ=255)

                                                                                                                      caacTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGa                               <  1:401529/150‑1 (MQ=255)

                                                                                                                      caacTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGa                               <  1:2267728/150‑1 (MQ=255)

                                                                                                                       aacTTCTGGCTGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCACTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGAc                              <  2:6408028/150‑1 (MQ=255)

                                                                                                                       aacTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGAc                              <  2:5093898/150‑1 (MQ=255)

                                                                                                                       aacTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGAc                              >  2:2064087/1‑150 (MQ=255)

                                                                                                                           tCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATa                          >  2:3502219/1‑150 (MQ=255)

                                                                                                                             tGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAAt                        <  1:6101169/150‑1 (MQ=255)

                                                                                                                              ggCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTg                                                                                                             >  1:3143634/1‑64 (MQ=255)

                                                                                                                                cGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCGGCGCCCAGCACTAAGGTTACTGCGACCATAATCaa                     <  1:5883245/150‑1 (MQ=255)

                                                                                                                                cGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCaa                     <  2:27438/150‑1 (MQ=255)

                                                                                                                                    ggAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTa                 >  2:4199999/1‑150 (MQ=255)

                                                                                                                                    ggAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTa                 >  2:3798345/1‑150 (MQ=255)

                                                                                                                                      aaCGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTAcc               <  2:5711525/150‑1 (MQ=255)

                                                                                                                                      aaCGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTAcc               >  1:1880858/1‑150 (MQ=255)

                                                                                                                                        cGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGt             >  1:875931/1‑150 (MQ=255)

                                                                                                                                         acgaaccaacTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGtt            <  2:4900567/147‑1 (MQ=255)

                                                                                                                                            aaccaacTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTg         <  2:3001427/150‑1 (MQ=255)

                                                                                                                                             accaacTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTgg        >  2:5956053/1‑150 (MQ=255)

                                                                                                                                               caacTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCt      >  2:1807194/1‑150 (MQ=255)

                                                                                                                                                aacTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCtt     <  1:6567310/150‑1 (MQ=255)

                                                                                                                                                   tttCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTaaa  >  2:4769406/1‑150 (MQ=255)

                                                                                                                                                   tttCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTaaa  >  2:4494105/1‑150 (MQ=255)

                                                                                                                                                   tttCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTaaa  <  2:686677/150‑1 (MQ=255)

ACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACCAATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCAGCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAA  >  NC_000913/1068410‑1068706

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