BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,612,260	A→T	5.9%	S45C (AGT→TGT)	ytjA →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,612,260	0	A	T	5.9%	1559.3 / 4.8	514	S45C (AGT→TGT)	ytjA	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base A (271/187); new base T (2/49); total (274/240)
Fisher's exact test for biased strand distribution p-value = 1.03e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 522 total aligned reads displayed.

TGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGA  >  NC_000913/4612111‑4612409

tGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGa                                                                                                                                                       <  2:4657231/150‑1 (MQ=255)

 ggTAAAGGAGACGCTTATTTTTAGTTTGGGCATAAAATTTATGGTTATAGCGTTAAAAGCCGCCGCAAGATGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAAATTTCTTTTTCGTCGGGATTATTATGTTACTTGTGAg                                                                                                                                                      >  1:10568362/1‑150 (MQ=255)

  gTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGt                                                                                                                                                     >  2:1125411/1‑150 (MQ=255)

   tAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCAGTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTGGGGATTATTCTGTTCCTGGTGAGtt                                                                                                                                                    >  2:3488841/1‑150 (MQ=255)

        gagaCGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTc                                                                                                                                               >  2:2677979/1‑150 (MQ=255)

          gaCGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCAt                                                                                                                                             >  1:10892351/1‑150 (MQ=255)

           aCGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATg                                                                                                                                            <  2:5601841/150‑1 (MQ=255)

            cGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTCATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATgg                                                                                                                                           <  1:5824493/150‑1 (MQ=255)

            cGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCAttg                                                                                                                                           <  1:2313740/150‑3 (MQ=255)

            cGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATgg                                                                                                                                           <  1:5573673/150‑1 (MQ=255)

             gCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATggg                                                                                                                                          >  1:9779412/1‑150 (MQ=255)

             gCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATggg                                                                                                                                          >  2:8242378/1‑150 (MQ=255)

              cTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGc                                                                                                                                         >  1:7493259/1‑150 (MQ=255)

               ttATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGcc                                                                                                                                        <  2:5824899/150‑1 (MQ=255)

                         ttGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGAcc                                                                                                                              <  2:708782/150‑1 (MQ=255)

                           ggggAATAATATTTCTGTTTATCGCTTTAATCGCCGCCGCACTTGGTTTTGTTGGTCTGGCCGTTACCGCTGCAGGCGAAGATAAAATTTTCTTTGTCGTCGGGATTATTCTTTTCCTGTTGTGTTTTTTCATGGGCCGAAAACGACCCt                                                                                                                            <  2:3183275/150‑1 (MQ=255)

                            gggCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTa                                                                                                                           >  2:752388/1‑150 (MQ=255)

                                     atTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGa                                                                                                                  >  1:3610610/1‑150 (MQ=255)

                                        tcgggTTATCGCGTTAATCGCCGCCGCCCTGGGGTTGGGGGGTCGGGCGGGTACCGCTGCCGGCGCCGCTAAAATTGTCTTGGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCa                                                                                                               <  2:4328031/147‑1 (MQ=255)

                                         cTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCAt                                                                                                              >  2:7810909/1‑150 (MQ=255)

                                          tgtttATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTGGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCAtt                                                                                                             <  2:2453311/147‑1 (MQ=255)

                                          gggTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTGTGGCCGGTACCGCTGCCGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTCCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCAtt                                                                                                             <  2:4887596/149‑1 (MQ=255)

                                             ttATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAg                                                                                                          >  1:2905787/1‑150 (MQ=255)

                                              tATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGa                                                                                                         <  2:3682543/150‑1 (MQ=255)

                                               aTCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGaa                                                                                                        >  2:9711526/1‑150 (MQ=255)

                                                tCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGaaa                                                                                                       >  2:2314461/1‑150 (MQ=255)

                                                 cgcgTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGaaaa                                                                                                      >  1:1332543/1‑150 (MQ=255)

                                                   cgTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAAtt                                                                                                    >  1:607106/1‑150 (MQ=255)

                                                     ttAATCGCCGCCGCACGTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCAGTACGAAAAtttt                                                                                                  >  2:2614434/1‑150 (MQ=255)

                                                       aaTCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTAt                                                                                                >  2:21325/1‑150 (MQ=255)

                                                        aTCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAAGTTTATc                                                                                               >  1:5485120/1‑150 (MQ=255)

                                                               ccgcACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGcc                                                                                        >  2:2699849/1‑150 (MQ=255)

                                                                   ccTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCCGGCGCAGCTAAAATTGTCTTGGTCGTGGGGATTATTCTGTTCCTGGTGTGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTc                                                                                    <  2:9098453/149‑1 (MQ=255)

                                                                   agttGGTTTTGGTGGTCGGGCGGGTCCCGCTGCCGGCGCTGATAATTTTTTCTTTGTCGTGGTTTTTATTCTGTTCGTGGTGTTTTTGTTCAGGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTc                                                                                    <  1:6868878/148‑1 (MQ=255)

                                                                     ttGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCa                                                                                  >  2:1100298/1‑150 (MQ=255)

                                                                       gggTTGGGTGGTGGGGCGGGAACCGGTGCCGGCGCCGCTAAATTGGTGTTTGTGGGGGGTTTTATTCGGTTCCTGGTGTTTGTGTTCTGGGGCCTAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGc                                                                                <  1:7943368/150‑1 (MQ=255)

                                                                        gtttgggtgttCTGGCCGGTACGGCTGCTGGCGCTGCTAAATTTTTCTTTGTCGGCGGGATTATTGTGTTCCTGGTGAGTTTGTTCAGGGCCCTAAAACGACCCTAGATTTCAGACAATTAAGAAAATTTTATCCAAAGCCAGTCCAGCg                                                                               <  1:10051423/145‑1 (MQ=255)

                                                                        ggTTTGGTGGTCGGGCCGGTACCGCTGCAGGCGAAGCTAAAATTGTCTTTGTCGTCGGGTTTTTTCTGTTCCTGGTGATTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCg                                                                               <  2:8432068/150‑1 (MQ=255)

                                                                         gTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCgg                                                                              >  1:4337155/1‑150 (MQ=255)

                                                                          tttGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGa                                                                             >  1:8082326/1‑150 (MQ=255)

                                                                          tttGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGa                                                                             >  1:2469324/1‑150 (MQ=255)

                                                                            gggtggtCGGGCCGGTACCGCTGCCGGCGAATCTAAAATTGTCTTTGTCGTGGGTATTATTCTTTTCCTGGTGATTTTGTTCATGGGCCCAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                           <  1:10425431/149‑1 (MQ=255)

                                                                             ggtggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCGGTTCCTGGTGAGTTTGTTCATGGGCCGAACACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTg                                                                          <  2:10508881/150‑1 (MQ=255)

                                                                              gtggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTgg                                                                         >  1:7349692/1‑150 (MQ=255)

                                                                               tggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCGTGGTGAGTTTGTTCATGGGCCGAAATCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAATGCCAGTCCAGCGGACTggg                                                                        >  2:5661501/1‑149 (MQ=255)

                                                                               tggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGc                                                                        >  1:4059453/1‑150 (MQ=255)

                                                                               tggtCTGGCCGGTACCGCAGCAGGCGCAGCTAAAATAGTATTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGAATTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGCCTGGc                                                                        >  1:1559146/1‑150 (MQ=255)

                                                                                 gtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCtt                                                                      >  1:6267372/1‑150 (MQ=255)

                                                                                 gtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTGGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCtt                                                                      <  2:1554892/150‑1 (MQ=255)

                                                                                  tcgggCCGGTACCGCTGCCGGCCCAGCTAAAATTGGCTTTGTCGCGGGGATTATTCTGTGCCTGGTGTGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCttt                                                                     <  2:10643298/147‑1 (MQ=255)

                                                                                  tCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCttt                                                                     >  2:6071960/1‑150 (MQ=255)

                                                                                     ggCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGc                                                                  >  2:3902289/1‑150 (MQ=255)

                                                                                     ggCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGc                                                                  >  1:9006277/1‑150 (MQ=255)

                                                                                      gCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCAGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCg                                                                 >  1:1136126/1‑150 (MQ=255)

                                                                                         ggTACCGCTGCAGGCTCACATAAAATTGTCTTTGTCGTGGGTATTATTCTGTGCCGGGTTTTTTGGTCCATGGCCCCAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGtt                                                              <  1:1950608/150‑1 (MQ=255)

                                                                                         ggTACAGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATATTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCgggt                                                              >  1:8538061/1‑148 (MQ=255)

                                                                                            aCCGCTGCCGGCGACGCTAAAATTGTCTTTGTCGTGGGTATTATTCTTTGCCGGGTGATTTTGTTCATGGGCCTAACACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTa                                                           <  1:10595430/150‑1 (MQ=255)

                                                                                            aCCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAATGCCAGTCCAGCGGACTGGCTGTTGCGGTTTTa                                                           >  1:8440928/1‑150 (MQ=255)

                                                                                                cTGCCGGCGCAGCTAAAATTGTCTTTTTCTTGGGGATTTTTCTGTTCCTGGTGAGTTTGTTCAGTGGCCGAAAACGACCATAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGa                                                       <  2:10312844/150‑1 (MQ=255)

                                                                                                 tGCTGGCGCTGATAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaa                                                      <  1:11061290/150‑1 (MQ=255)

                                                                                                 ggCAGGCGCAGATAAAATTGTCTGGTGCGGCGGTTTTATTCTGTTCCTGTTGTGTTTGTACAGGGGCCTAACACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCATTTGCGGTTTTAGCGaa                                                      <  1:4743239/149‑1 (MQ=255)

                                                                                                      gcgcAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAataata                                                 >  1:4935425/1‑150 (MQ=255)

                                                                                                        gcAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATaaa                                               >  1:387479/1‑150 (MQ=255)

                                                                                                         cAGCTAAACTTGTGTGTGTGGGGGGGTTTTTTGTGTTCGGGGTGTGTTTGTTCTGGGGCCGAACACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAAt                                              <  1:3161313/150‑1 (MQ=255)

                                                                                                           gCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTg                                            >  1:3312090/1‑150 (MQ=255)

                                                                                                             tAAATTTTTCTTTGTCGTCGTGATTATTCTTTTCCTGGTTTTTTTGTTCTTGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCg                                          <  2:1372279/150‑1 (MQ=255)

                                                                                                             tAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCGGGAGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCg                                          <  2:6381446/150‑1 (MQ=255)

                                                                                                                atttGTCTTTGTGGTCGGTTTTATTCTTTTCGTGGTGTGTTTTTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGTCTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTa                                       <  2:3349649/148‑1 (MQ=255)

                                                                                                                aaTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTa                                       >  1:8861467/1‑150 (MQ=255)

                                                                                                                    gTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACttt                                   <  1:6527194/150‑1 (MQ=255)

                                                                                                                       tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTc                                >  1:6723293/1‑150 (MQ=255)

                                                                                                                        ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCa                               >  1:9916559/1‑150 (MQ=255)

                                                                                                                        ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCa                               <  2:5504962/150‑1 (MQ=255)

                                                                                                                        ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCa                               <  2:7851412/150‑1 (MQ=255)

                                                                                                                        ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCa                               >  2:2012788/1‑150 (MQ=255)

                                                                                                                           tcgtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTc                            >  1:699074/1‑150 (MQ=255)

                                                                                                                            cgtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCAGGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCt                           <  2:7278677/150‑1 (MQ=255)

                                                                                                                                ggtattattCTGTTCGTGTTGATTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTca                       <  1:5311464/147‑1 (MQ=255)

                                                                                                                                gggATTATTCTTTTCCTTGTGAGTTTGTTCATGGGCCTAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTca                       <  2:4130636/150‑1 (MQ=255)

                                                                                                                                  tattattCTGTTCCTGGTGTGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCATACCATTAAGAATATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcaca                     <  2:8714427/149‑1 (MQ=255)

                                                                                                                                  gATTATTGTGTTCCTGGTGAGTTTGTTCATGGGCCGAAATCGACGCTAGATTTCAGAGCATTATGAAAATTTTATCCAAAGCCAGTTCAGGGGAGTGGGTTTTGGGGTTTTAGGGTATAATAAATTGCGTTAGTTTGTCATTCTTTctca                     >  2:1885783/1‑150 (MQ=255)

                                                                                                                                  gATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcaca                     >  2:8522344/1‑150 (MQ=255)

                                                                                                                                  gATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcaca                     >  1:7723455/1‑150 (MQ=255)

                                                                                                                                        tCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaa               >  2:6966896/1‑150 (MQ=255)

                                                                                                                                        tCTGTTCCTGGTGAGTTTGTTCATGGGCCAAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaa               >  2:8635844/1‑150 (MQ=255)

                                                                                                                                          tGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAAc             >  2:4522891/1‑150 (MQ=255)

                                                                                                                                          ggTTCCGGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAAc             <  1:1828272/149‑1 (MQ=255)

                                                                                                                                           gTTCCTGGTTAGTTTGTTAATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa            <  2:3266385/150‑1 (MQ=255)

                                                                                                                                            ttCCTGGTGTGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCGTTCTTTCACAACAAAAACAg           <  2:8872357/150‑1 (MQ=255)

                                                                                                                                            ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTATGTAAATTTTATCCATTGCCAGTCCAGCGGAGTGGCTTTTGCGGTTTTAGCGAATATTAAATTGCGTTACTTTGTCATTCTTTCAGAACAAAAACAg           >  2:5196575/1‑150 (MQ=255)

                                                                                                                                            ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAg           >  1:8185051/1‑150 (MQ=255)

                                                                                                                                            ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAg           >  1:2664055/1‑150 (MQ=255)

                                                                                                                                            ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAg           <  1:9305566/150‑1 (MQ=255)

                                                                                                                                            ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAg           >  1:8582741/1‑150 (MQ=255)

                                                                                                                                             tCCTGGTGCGTTTGTTCTTGGGCCTAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAgg          <  2:2933017/150‑1 (MQ=255)

                                                                                                                                                tGGTGAGTTTGTTCATGGGCCGAACACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg       <  2:7645449/150‑1 (MQ=255)

                                                                                                                                                tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg       >  2:6817758/1‑150 (MQ=255)

                                                                                                                                                tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg       >  1:10649833/1‑150 (MQ=255)

                                                                                                                                                tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg       >  2:9450058/1‑150 (MQ=255)

                                                                                                                                                 ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAgg      >  2:10126504/1‑150 (MQ=255)

TGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGA  >  NC_000913/4612111‑4612409

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