BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,286,944	T→G	7.4%	intergenic (+418/+122)	narI → / ← rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,286,944	0	T	G	7.4%	815.5 / 42.0	312	intergenic (+418/+122)	narI/rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand): ref base T (104/167); new base G (24/2); total (142/170)
Fisher's exact test for biased strand distribution p-value = 5.72e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

GTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCG  >  NC_000913/1286795‑1287092

gTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCt                                                                                                                                                      >  1:4442572/1‑150 (MQ=255)

 tCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGGGAATTCCGGGGCTTTGGAGAGTGGTGGTGGTGGTGGAAGGATTAATCAGCTCTGCGCGCGTCGCCCTTAGGGTCGTTGGCTg                                                                                                                                                     >  2:5721122/1‑150 (MQ=14)

 tCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                                                                                                                                                     >  1:4822716/1‑150 (MQ=255)

   gAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTTCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCg                                                                                                                                                   >  2:884902/1‑150 (MQ=255)

   gAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCg                                                                                                                                                   <  1:936616/150‑1 (MQ=255)

    caCCTTGGTGTAAGCTTCTCATCCCTCCCCGCTTTGGCAGATTAATAGATTTCTGATTAGTATGAGAATTCCTGGGCTTATGTAAGTGATGGTGGGGGGTGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                                  <  1:6741176/149‑1 (MQ=14)

    caCCTGGGTCGAAGCTTCTCAGCCGTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAGGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                                  <  2:2320863/149‑1 (MQ=14)

    aaCCTTGGTCGAAGCTTCGCATCCTTCCCCGCGTTGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                                  <  1:3340185/150‑1 (MQ=255)

       cTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                               <  1:2681217/150‑1 (MQ=255)

       cTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                               <  1:5746000/150‑1 (MQ=255)

         tGGTCGAAGGTTCTGATCGGTCCCCGCTTGGGCAGATTTTTTTATTGCGGAGTAGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGCAGGCAACg                                                                                                                                             <  2:5897339/150‑1 (MQ=14)

          ggTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc                                                                                                                                            <  1:3782879/150‑1 (MQ=255)

                 gCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCg                                                                                                                                     <  1:6537335/150‑1 (MQ=255)

                  cgtctcATCCTGCCCCGCGGGGGCAGAATTATTGATTGCGGCTTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGc                                                                                                                                    <  1:5061138/148‑1 (MQ=14)

                  cTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGc                                                                                                                                    >  2:8263217/1‑150 (MQ=255)

                    tctcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                                                                                                                  >  1:1397956/1‑150 (MQ=255)

                     atcATCCGTCCCCGCTTGGGCAGATTATTTGATTGCGGATTCGTTAGAGAATTCCGGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGa                                                                                                                                 <  2:4259268/149‑1 (MQ=255)

                      tcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCGCCCACGCTgtc                                                                                                                                >  1:49238/1‑148 (MQ=14)

                        aTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGc                                                                                                                              >  1:4195534/1‑150 (MQ=255)

                        aTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCTCTCACGCTGACGc                                                                                                                              >  1:364992/1‑150 (MQ=255)

                            ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCGCACGCTGACGCTCGa                                                                                                                          >  1:72637/1‑150 (MQ=255)

                            ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAATTCCGGGGCTGGTGAAAGAGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCAGCGGCACCGCGCGCTCGCGGAcgcgcgc                                                                                                                          >  2:7851035/1‑149 (MQ=14)

                             tCCCCGCTGGGGCAGAATATTTGAGTGCGGAGTCGTTTGAGAATTCCGGGGCGTTTGAAAGTGATGGTGGTGGGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                                                                                         <  1:5805518/150‑1 (MQ=14)

                               cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                       <  2:6687926/150‑1 (MQ=255)

                                ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                                      >  1:693032/1‑150 (MQ=255)

                                 cGCTTGGGCAGAATCTATGATTGCGGATTCGTTTGAGAATTCCGGGGATTTAGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                     <  1:1274086/150‑1 (MQ=255)

                                 cGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                     <  2:3955266/150‑1 (MQ=255)

                                   cTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCATGCGGCAACGCACTCTCGCTGTCGCTCGAGACGAAc                                                                                                                   >  1:8241647/1‑150 (MQ=255)

                                    ttGGGCAGAAGATTTGATTGCGGATTCGTTTGAGTATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAcc                                                                                                                  <  1:7937460/150‑1 (MQ=255)

                                         ctgAATTTTTGCGTGCGGTTTGGTTTGAGTATTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGt                                                                                                             <  1:1846464/148‑1 (MQ=14)

                                          aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACAAACATCTGCGCCCTTCTCCCAACGGGACGATTCAGACCGCAACGCACTCAAGCGGACCCCCGAGCCCCACCGTAGTc                                                                                                            >  1:3231031/1‑150 (MQ=9)

                                            aaTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACACGCTGACGCTCAAGCCGTAAAATAGTCGa                                                                                                          >  2:2778581/1‑150 (MQ=14)

                                              tatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCATAGTCGAAg                                                                                                        >  2:4357996/1‑150 (MQ=255)

                                               atTTGATTGCGGATTAGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGc                                                                                                       <  2:6413124/150‑1 (MQ=255)

                                                tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCACACGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                                                                                      >  2:3022781/1‑150 (MQ=255)

                                                     tttcGGATTCGTTTGAGTATTCCGTGTCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                                                                                 <  2:5113406/147‑1 (MQ=14)

                                                     ttGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                                                                                 <  1:2562904/150‑1 (MQ=255)

                                                     ttGCGGATTCGTTTGAGAAATCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTTGCCAGCGGCAACGCACACACCCAGACGCACGAGCCGAACCCTAGTCGAAACTtcacc                                                                                                 >  1:7588377/1‑147 (MQ=14)

                                                       gCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATc                                                                                               <  1:4002688/150‑1 (MQ=255)

                                                          gATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACACATCGCACCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACTCTCTCTCGCTGGCCCACGAGACGAACCAAAGTCGAAGCTTCTCATCCtt                                                                                            >  2:4794926/1‑150 (MQ=14)

                                                             tCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                         <  1:5936318/150‑1 (MQ=255)

                                                                tttGAGTATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGc                                                                                      <  2:6232081/150‑1 (MQ=255)

                                                                tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGACGCTCGAGCCGAACCTTATTAGAAGCATCTCAACCATcccccc                                                                                      >  1:1591945/1‑148 (MQ=14)

                                                                tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCACGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTATAAACCATCCCCGc                                                                                      >  2:1832425/1‑150 (MQ=14)

                                                                  tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                    <  2:2997423/150‑1 (MQ=255)

                                                                  tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                    <  2:6326890/150‑1 (MQ=255)

                                                                  tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCAGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                    >  2:6904480/1‑150 (MQ=255)

                                                                  tGAGAATTCCGGGGATTTTGTTATTGTTGGTGGTGGGTGAAGGACTACTCTGCGCTGCGAGCTTCGCCCTTCGTGTCGTTTCCTTAGTCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                    <  2:3410063/150‑1 (MQ=14)

                                                                   gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTAGCTGACGCTCCAGCCGAACCGTAGTCGAAGCTTATCATACATCCCCGCATg                                                                                   >  2:1457194/1‑150 (MQ=14)

                                                                   gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACAAACAGCTGCGCCCTTCGCCCATCTGTTCGTTGCATGAGGCACAGCGCTCTCGCGGACGCGCGATCCGGACCTAGGTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                   >  2:6137880/1‑150 (MQ=14)

                                                                    agaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                  <  1:5546361/150‑1 (MQ=255)

                                                                      aaTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGAGCTTCGACCTTCGGGTCGATTCCTTCCTCAAAGCTCTATAGCTGGCGCTCCACCCGAACCTTAATAGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                >  2:2420602/1‑150 (MQ=11)

                                                                            ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                          <  1:708588/150‑1 (MQ=39)

                                                                             gggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAatat                                                                         >  2:5369314/1‑150 (MQ=255)

                                                                              ggCTTTTGTAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                        <  2:2669813/150‑1 (MQ=255)

                                                                                 ttttGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                     <  2:8324990/150‑1 (MQ=38)

                                                                                 ttttGATAGTTATGGTGGTGGGGGATGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                     <  1:1519627/150‑1 (MQ=255)

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                     <  1:5361236/150‑1 (MQ=38)

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCATCGGCAACCCTCTAACCCAGACGCACGAGCCCAACAATAGTCGAAGCATATCAAAAATCACCCCAAGGGCAGAATATTTGa                                                                     >  1:4353062/1‑150 (MQ=2)

                                                                                   ttGTAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                   <  1:2222154/150‑1 (MQ=255)

                                                                                       aaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCgg                                                               <  1:1958710/150‑1 (MQ=39)

                                                                                       aaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTAGTAGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTGATTGCgg                                                               >  2:2618962/1‑150 (MQ=11)

                                                                                           gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                           <  1:7884196/150‑1 (MQ=39)

                                                                                           gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTc                                                           >  2:6739482/1‑150 (MQ=21)

                                                                                           gATGGTGGTGGGGGAAGGACTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                           <  2:7188566/150‑1 (MQ=255)

                                                                                            aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCg                                                          <  2:2536540/150‑1 (MQ=39)

                                                                                            aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAAAATTTGATTGctattctt                                                          >  2:1966362/1‑143 (MQ=255)

                                                                                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                    <  2:316267/150‑1 (MQ=255)

                                                                                                   tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                   <  1:719341/150‑1 (MQ=255)

                                                                                                        gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                              <  2:6199489/150‑1 (MQ=2)

                                                                                                        gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                              <  2:6131515/150‑1 (MQ=255)

                                                                                                        gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                              <  1:2736958/150‑1 (MQ=255)

                                                                                                         aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                             <  2:5700238/150‑1 (MQ=255)

                                                                                                         aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                             <  1:7851035/150‑1 (MQ=255)

                                                                                                         aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTACCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                             <  1:8400648/150‑1 (MQ=255)

                                                                                                          aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                            <  2:5018954/150‑1 (MQ=255)

                                                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                          <  1:4138679/150‑1 (MQ=255)

                                                                                                                 cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTc                                                                                                            <  1:5608648/79‑1 (MQ=255)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGACGAACCTTAGTCGGAGCTTCTCATCCTTCCCCGCATGGGCAGTATATTTGATTGCGGATTCGGTGGGGAGTTCAGGGACGTTTg                                    >  2:3291688/1‑150 (MQ=11)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTATCGCGGACGCGCGACCCGAACCATAGTCGAAGCGTATCAACCGTCCCCGCATGGGGAGAAGATGTGAGTGCGGAGTCGCGTGAGAGGGAAGGGACTGTTg                                    >  2:7067121/1‑150 (MQ=9)

                                                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCGCTCTCGCTGACGCGCGAGCCGAACCGTAGTCGAAGCTTCTCAACCGTACCCCCAGGGGCAGAAAATTTGATTGCGGATTCGCGTGGGAGTTCAGGGACGTTTgta                                  >  2:2152617/1‑148 (MQ=14)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                       >  2:5298619/1‑60 (MQ=40)

                                                                                                                          gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGGTGGGGAGTTCAGGGACTTTTGAAAGTGAT                            >  2:3497476‑M1/1‑150 (MQ=255)

                                                                                                                           cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCCGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACTTTTGAAAGTGATG                           >  2:7657382‑M1/1‑150 (MQ=255)

                                                                                                                           cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATG                           <  1:5702499‑M1/150‑1 (MQ=255)

                                                                                                                                    cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGG                  <  1:6176690‑M1/150‑1 (MQ=255)

                                                                                                                                      cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGA                <  2:5565590‑M1/150‑1 (MQ=255)

                                                                                                                                        tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG              <  1:2987034‑M1/150‑1 (MQ=255)

                                                                                                                                        tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG              <  2:1846362‑M1/150‑1 (MQ=255)

                                                                                                                                         gagagcgttgcCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGAAAGTGATGTTGGTGGGGGAAGG             >  2:8125435‑M1/6‑150 (MQ=255)

                                                                                                                                          agagcgttgCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGA            >  1:8068774‑M1/5‑150 (MQ=255)

                                                                                                                                            gTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATT          >  2:2009458‑M1/1‑150 (MQ=255)

                                                                                                                                            gTCGTTGCCTGCGGCAACGCTCTCTAGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTTGGGAGTTAATGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATT          >  2:4454788‑M1/1‑150 (MQ=255)

                                                                                                                                            agcgtTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGGGAGTTCAGGGACTGTTGAAAGTGATGGTGGGGGGGGAAGGATa          >  2:8190741‑M1/3‑149 (MQ=255)

                                                                                                                                             tCGTTGCCTGCGGCACCGCACTCCCGCTGTCTCTCGAGCCGAACCTTAGTCGACGCTTCTCATCGCTCCCCGGATGGGCAGACTTTTTGATTGGGGCTTCGATTGAGCGTTGTTTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa         <  2:751113/150‑1 (MQ=14)

                                                                                                                                             gcgTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGGTCGTTCAGGGACGTGTGAAAGTGATGGTGGTGGGGGAAGGATTc         >  2:8024691‑M1/2‑149 (MQ=255)

                                                                                                                                             gcGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGCCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa         >  2:3832491/2‑150 (MQ=255)

                                                                                                                                              cGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTATAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCTCTTGATAGTTCATGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAc        <  1:2402636/150‑1 (MQ=17)

                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTATTGAAAGTGATGGTGGTGGGGGAAGGATTACTc      >  2:4414017/1‑150 (MQ=255)

                                                                                                                                                  gCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAg    <  1:3054085/150‑1 (MQ=255)

GTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCG  >  NC_000913/1286795‑1287092

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