BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,122	T→G	5.8%	noncoding (115/171 nt) R19R (AGG→CGG)	rttR ← tpr ←	rtT sRNA, processed from tyrT transcript protamine‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,287,122	0	T	G	5.8%	831.4 / 13.4	350	noncoding (115/171 nt) R19R (AGG→CGG)	rttR tpr	rtT sRNA, processed from tyrT transcript protamine‑like protein
Reads supporting (aligned to +/- strand): ref base T (161/153); new base G (26/3); total (193/157)
Fisher's exact test for biased strand distribution p-value = 4.95e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 375 total aligned reads displayed.

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286973‑1287271

gCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCt                                                                                                                                                       >  2:8764283/1‑150 (MQ=255)

   gAACCCAAGTCGACCCTTCTCGTCCATCCCCGGATCTGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCg                                                                                                                                                    <  2:5080915/150‑1 (MQ=11)

    aaCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                                   <  1:10966874/150‑1 (MQ=255)

        ctcgTCGAAGCTCCTCCTCCTGCCCCGCATGGGCAGAATATTAGATTGCGGATTCGCTTGAGAGTTCAGGGACTATTGCAAGGGATGTTGGTGGGGGAAGTATTACTCCGGGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAc                                                                                                                                               <  1:7309662/147‑1 (MQ=14)

           gTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGc                                                                                                                                                     <  2:11270182/141‑1 (MQ=255)

            tCGAAGCTTCTCATCCTTCCCCGCTTGGGCTTAATTTTTGATTGCGGTTTGGCTTGTGTGTTCAGGGACTTTTTAATGTTTTGGTGGTGGGGGATGTTTTACTCAGCGGTGCGCGCTTCGCCCTTGGGGTCGTTGCCGCAGGCAACgacc                                                                                                                                           <  2:8150880/150‑4 (MQ=2)

            tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGGGGGGTTACTCTGCGCGGCGCGCTTCGCGCTGCGGGGCGTAGCCggcggcgccgcgc                                                                                                                                           >  2:4017274/1‑143 (MQ=14)

              gAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                                                                                                                                         >  1:6587032/1‑150 (MQ=37)

              gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCGGGCAATTTTGAAATTGATGCTGGTGGTTGAAGGATTACTAAGAGCAGCGCGCTTCGCCCTACTGGTCGTTTCCTGCGGCAACGcgcgc                                                                                                                                         >  2:7314214/1‑146 (MQ=14)

                 gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCTTGCGGCAACGCTCTCTCg                                                                                                                                      >  1:5965326/1‑150 (MQ=255)

                   ttCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGcg                                                                                                                                    >  1:3960731/1‑149 (MQ=255)

                      tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgc                                                                                                                                 >  2:6300726/1‑150 (MQ=255)

                      tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgc                                                                                                                                 >  2:11076728/1‑150 (MQ=255)

                      tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGgc                                                                                                                                 >  2:10855455/1‑150 (MQ=255)

                        aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcgc                                                                                                                               >  2:12105421/1‑150 (MQ=255)

                          ccgtCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                             <  2:8973591/147‑1 (MQ=255)

                           cTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCg                                                                                                                            >  1:5857023/1‑150 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                           >  1:824396/1‑150 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                           <  2:7620186/150‑1 (MQ=255)

                              ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt                                                                                                                         <  1:8463378/150‑1 (MQ=255)

                                 cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGa                                                                                                                      <  1:8516503/150‑1 (MQ=255)

                                  gCATGGGCCGAATATTTGATAGCGGATTCGCTTGAGACTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCAGCGCTGCGCGCTTCTCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                     <  1:6398835/150‑1 (MQ=14)

                                      gggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCtt                                                                                                                 <  1:3055949/150‑1 (MQ=255)

                                      gggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCtt                                                                                                                 >  1:7012021/1‑150 (MQ=255)

                                       ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg                                                                                                                >  1:8852676/1‑150 (MQ=255)

                                       ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCACGCTGGCGCACGCGTCGAACCtgg                                                                                                                >  1:4350857/1‑148 (MQ=255)

                                             atatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGaa                                                                                                          <  2:6587032/150‑1 (MQ=255)

                                              tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAg                                                                                                         >  2:10088768/1‑150 (MQ=255)

                                                tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAAAGCACAATCGCTTGCGCTCGAGTAGAACCTTGGTCGAAGCt                                                                                                       >  2:2292051/1‑150 (MQ=14)

                                                      tGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCTCACGCTTGCGCTCGAGTCGAACCATTGTCGAAGCTTCTCaa                                                                                                 >  1:9138344/1‑149 (MQ=14)

                                                       gCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                <  1:1613343/150‑1 (MQ=255)

                                                        cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCGCGAGTCGAACCTTGGTAGAAGCTTATAATcc                                                                                               >  2:12318/1‑150 (MQ=255)

                                                           aTTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCGGCGGCAACGCTCGCTAGCTTGCGCACGAGTCGAAACATGGTCGAAGCTTATCATACtaa                                                                                            >  1:4282056/1‑148 (MQ=14)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                      >  2:7444290/1‑150 (MQ=255)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTACGGGTCGTTTCCAGCGGCAACGCTCACTCGCAGGCGCTCGAGTCGAACCATAGTAGAAGGTTCACACCCTTcccccaa                                                                                      >  2:3812839/1‑147 (MQ=14)

                                                                  tGAGAGTTGAGGTACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                     <  1:1251898/150‑1 (MQ=255)

                                                                      atttCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                 <  1:12178032/148‑1 (MQ=255)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCATCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCAGCGGCAAAGCGCAATCCCATGCGCACGAGTCCAACCGTGGTCGAAGCATCTCATCCATACCCCCAAGGGCa                                                                                >  1:10223259/1‑150 (MQ=14)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCAGGGGCa                                                                                >  2:11189009/1‑150 (MQ=255)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCAGGCGCGCGAGTCGAAACATGGTCGAAGCTTAAAAAACTTCCCCGCAGGGGCCGAat                                                                            >  1:7598239/1‑150 (MQ=14)

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGGGGCAACGCTCTCTCGCTGTCGCACGAGTCGAACCTTGGACGAAGCTTCTCAACCTTCCCCCCATCGGCCGaaaaaaaaa                                                                      >  1:7926089/1‑143 (MQ=11)

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCATACCCGCATGGGCAGAATATTTaa                                                                      >  2:5973643/1‑150 (MQ=255)

                                                                                     gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                  <  1:3738462/150‑1 (MQ=21)

                                                                                      aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTAGCTGGCGCTAGAGTCGAACCTTGGTCGTATCTTATCATCCTTCCCCGCATGGGCAGAAATTTTaaattcg                                                                 >  2:11631136/1‑145 (MQ=11)

                                                                                      aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCACACACCCAGGCGCTCGACTCGAACCAAGGTCGAAGCATCAAAAACAACACCCCAAGGGCAGAATAATTaaattag                                                                 >  1:8756127/1‑145 (MQ=2)

                                                                                          cgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                     >  2:7922298/2‑94 (MQ=17)

                                                                                               gtgtcgggggCAGGAGTTCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTg                                                        <  1:9858061/145‑1 (MQ=11)

                                                                                               gtggtggGGGAAGGAGTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTg                                                        <  1:9630979/150‑1 (MQ=38)

                                                                                                tggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTgg                                                       <  1:2515835/150‑1 (MQ=12)

                                                                                                tggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTgg                                                       <  2:7012021/150‑1 (MQ=38)

                                                                                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                      >  2:2952025/1‑133 (MQ=37)

                                                                                                   tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGGGTCGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTAATTGCGGTTTCGGTGGGaa                                                    >  2:533548/1‑150 (MQ=11)

                                                                                                   tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                     <  1:10572342/101‑1 (MQ=32)

                                                                                                    gggggCAGGAGTACTCAGCGCGGCGCGCGTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCGGGCGCTCGCGTCGCACCGGGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTGGGGAAg                                                   <  2:10332074/150‑1 (MQ=14)

                                                                                                         aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                              <  2:6796020/150‑1 (MQ=255)

                                                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCGCAATCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATCCGTCCCCGCCTGGGCAGAAAATTTAATTGCGGTTTCGTTGGGCAGTTCAGGGa                                           >  2:8373283/1‑150 (MQ=14)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGTAGTTCAg                                              >  2:962847/1‑141 (MQ=255)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                              <  1:6191475/141‑1 (MQ=255)

                                                                                                                     gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                  <  2:4634196/150‑1 (MQ=255)

                                                                                                                      cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAg                                 <  2:2590052/150‑1 (MQ=255)

                                                                                                                      cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATACGTCCCCGCATGGGCCGAATATTTAATTGCGGATTCGGTGGGGAGTTCAGGGACTTTTGAAAg                                 >  2:3432318/1‑150 (MQ=11)

                                                                                                                        cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTg                               <  2:10545323/150‑1 (MQ=255)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                              <  1:5152094/150‑1 (MQ=255)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCAACCTTCCCCGCATGGGCAGAAAATTTAATTTCGGGTTCGTTTGGAAGTTCAGGGACTTTTGTAAGTgt                              >  2:5424281/1‑149 (MQ=14)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGGATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                              <  2:8469471/150‑1 (MQ=255)

                                                                                                                              gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGGCGCTCGGGCCGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGCCGGGGCCGAATATTTAATTGCGGATTCGGGGGGGAGTTCAGGGACTTTTGAAAGTGAggggg                         >  2:9461692/1‑148 (MQ=11)

                                                                                                                                ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                       <  2:4256624/150‑1 (MQ=255)

                                                                                                                                ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGt                                                                                                              <  1:11874952/63‑1 (MQ=35)

                                                                                                                                ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                       <  2:6108036/150‑1 (MQ=255)

                                                                                                                                  cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAATTTCAGGGACTTTTGAAAGTGAtggtggtgg                     <  1:11289374/150‑1 (MQ=255)

                                                                                                                                    cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGTACCTTGGTCGAATCTTATCATACTTACCCGCATTGGCAGAATATTTAATTGAGGTTTAGTTGGGAATTTGAGGGACTTTTGAAATTGATTGTGGTgggg                   >  2:9068150/1‑150 (MQ=17)

                                                                                                                                     ccTTCGGGTCGTTGCCGGCGGCAAAGCTCTCTCCCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTATCATCCTTCCCCCCAAGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGCAAGTGATGGTGGTggggg                  >  2:1429247/1‑150 (MQ=17)

                                                                                                                                     ccTTCGGGTCGGTGCCGTCGGGAACGCTCTCGCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                  <  1:10396406/150‑1 (MQ=17)

                                                                                                                                       ttCGGGTCGTTGCCTTCGGCAACGCTCTCTCGCTGGCGCTCGAGTGGAACCTTGGTCTTATCTTCTGTTCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTTGTAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                <  1:11148023/150‑1 (MQ=14)

                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                <  1:6741788/150‑1 (MQ=255)

                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                <  1:3494417/150‑1 (MQ=255)

                                                                                                                                          gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa             <  1:9461692/150‑1 (MQ=255)

                                                                                                                                          gggTCGTGGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTGTGTTTGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa             <  1:12103739/150‑1 (MQ=11)

                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc          <  1:2238363/150‑1 (MQ=255)

                                                                                                                                             gcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCCCATGGGCAGAAGATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTGTTGAAAGTGATGGTGGTGGGGGGAGGATTc          >  1:3527144/2‑150 (MQ=17)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa        <  1:4543154/150‑1 (MQ=255)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa        <  1:4685491/150‑1 (MQ=255)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTATCATCCTTCCCCGCATGGGCCGAATATTTAATTTCGGATTCGTTGGGAAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGTATGATTAGa        >  1:8140898/1‑150 (MQ=17)

                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                              <  1:10041488/127‑1 (MQ=255)

                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa       >  2:6440728/1‑150 (MQ=255)

                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCCCCTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa       >  1:7560811/1‑150 (MQ=255)

                                                                                                                                                 tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCAGCGCATGGTCATAATATTTTATTGCGGATTCGTTGGGAAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAc      <  1:11666934/150‑1 (MQ=17)

                                                                                                                                                 ggCCTGCGGCACCGCGCTCCCGCTGGCGCTCGAGTCGAACCGTGGTCGACGCTTCTGATCCTTCCCCGCTTGGGCAGAATATTTTTTTGCGGATTCGTTGGGAAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAc      <  1:493452/149‑1 (MQ=14)

                                                                                                                                                   ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                     <  2:11778961/37‑1 (MQ=40)

                                                                                                                                                   ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCt    >  2:5913081/1‑150 (MQ=255)

                                                                                                                                                   ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCt    >  2:211550/1‑150 (MQ=255)

                                                                                                                                                   ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATACTTCCCCGCCTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTCGAACCt    >  1:7685321/1‑150 (MQ=255)

                                                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCtt   >  2:8943650/1‑150 (MQ=255)

                                                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTg                               <  1:3816790/122‑1 (MQ=255)

                                                                                                                                                     tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  <  1:533666/150‑1 (MQ=255)

                                                                                                                                                     tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  <  1:8574452/150‑1 (MQ=255)

                                                                                                                                                     tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  >  2:366130/1‑150 (MQ=255)

                                                                                                                                                     tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  >  2:1728617/1‑150 (MQ=255)

                                                                                                                                                     tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCATCCCCGCATGGGCAGAATATTTAATTTCGGATTCGGTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGGTGGATTAGAACCTTc  >  1:8074941/1‑150 (MQ=17)

                                                                                                                                                     tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  <  2:1158731/150‑1 (MQ=255)

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286973‑1287271

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