BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,122	T→G	9.4%	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	9.4%	300.7 / 2.7	147	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 (67/54); new base G (16/2); total (90/57)
Fisher's exact test for biased strand distribution p-value = 8.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

CCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286974‑1287271

ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                                                                                                                                                      >  2:3581045/1‑150 (MQ=255)

 cgcaCCTTCGTCGAAGCGGCTCAGCCGGCCCCGCATGGGCAGCCTCTTTGCTTGCGGCTTCGCTTGAGCGTTCCGGGACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGc                                                                                                                                                     <  2:2683966/147‑1 (MQ=11)

  gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTTATGTTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCg                                                                                                                                                    >  1:4298664/1‑150 (MQ=255)

  gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCGTGTGAGTTCAGGGACTTTTGAAAGTGGTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCAACGGGTCGTTGCCGGCg                                                                                                                                                    >  1:1463343/1‑150 (MQ=14)

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

      cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                                >  1:2044329/1‑150 (MQ=255)

      cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGAGCGCTTCGCCCGTCGGGTCGTTTCCGGCCGCaa                                                                                                                                                >  1:2044312/1‑150 (MQ=14)

       tttgTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAc                                                                                                                                               <  1:2250533/147‑1 (MQ=255)

       ttggTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAc                                                                                                                                               <  1:315293/147‑1 (MQ=21)

       ttAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCATCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAAc                                                                                                                                               >  1:350514/1‑150 (MQ=255)

       ttAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAc                                                                                                                                               <  1:2519539/150‑1 (MQ=255)

       ttAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGAGCTTCGCCCTTCGGGTCGTTTCATGAGGCAAc                                                                                                                                               >  2:3229739/1‑150 (MQ=14)

             gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                                                                                                                                         >  1:1922416/1‑150 (MQ=255)

                   tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                                                                                                                   >  1:1840788/1‑150 (MQ=255)

                   tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                                                                                                                   >  1:3051756/1‑150 (MQ=255)

                   tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCGCTCTCGcgg                                                                                                                                   >  1:259435/1‑148 (MQ=255)

                     tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCTCCCTACGGGTCGTGGCAGGGGGCAAAGCGCTCTCGCGGgc                                                                                                                                 >  1:1965350/1‑150 (MQ=14)

                     tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTGCCTGAGGCAACGCGCTCTCGCTgtc                                                                                                                                 >  2:4186355/1‑148 (MQ=255)

                      cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCGTCGGGTCGTTGCCGGCCGCAAAGCTCTCTCGCTGgcg                                                                                                                                >  2:719995/1‑150 (MQ=255)

                      cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAAAAAGCGCTGCGCGCTTCGCCATTCGGTTCGTTTCCTGCGCCAAAGCTCTCTCGCTGgcc                                                                                                                                >  2:719981/1‑149 (MQ=14)

                        tCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACACAGAGCTGCGCGCTTCGCCCATCGGGTCGTTGCCTGACGCAACGCTCTCTCGCTGgcgcg                                                                                                                              >  2:3234520/1‑149 (MQ=14)

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

                           ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                           >  2:3127354/1‑150 (MQ=255)

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

                           ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCAGGCGCTCGa                                                                                                                           >  2:1886432/1‑150 (MQ=255)

                           ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAACCGCTGCGCGCATCGCCCTTCGGGTCGGTGCCTGAGCCAACGCGCTCTCGCTGGCCCACGa                                                                                                                           >  2:2857272/1‑150 (MQ=14)

                              cccGCCTGGGCAGACTCTTTGCGTGCGGCGGCGCTGGAGCGGTCCGGGACTTTTGACAGTGAGGGGGGTGGGGGTAGGAGTCCTCCGCGCTGGGCGCTTCGCCCTTCGGGTCGGGGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTc                                                                                                                        <  1:2373418/150‑1 (MQ=2)

                               ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCAGCGGCAACGCACTCTCGCTGGCGCTCGAGTCg                                                                                                                       >  1:3312712/1‑150 (MQ=255)

                               ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTCCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCg                                                                                                                       >  2:3270958/1‑150 (MQ=255)

                                    tGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGGGGCAAAGCTCTCTCGCTGGCGCTCGATTCggaacg                                                                                                                  >  2:403626/1‑145 (MQ=14)

                                       gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAATCGAACCTTgg                                                                                                               >  1:2105743/1‑150 (MQ=255)

                                        cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCTCTCTCGCTGGCGCACGAGTCGaaacattgt                                                                                                              >  1:1278281/1‑143 (MQ=14)

                                            atatTAGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTCATCGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGaa                                                                                                          <  1:3962497/150‑1 (MQ=255)

                                             tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAg                                                                                                         <  1:1068354/150‑1 (MQ=255)

                                               tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGGAACCCACTCAACCAGGCCCGCGACTCGACCCTTGGTCGACGCt                                                                                                       >  1:4173445/1‑150 (MQ=14)

                                                ttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCtt                                                                                                      <  1:1703262/150‑1 (MQ=255)

                                                   aTTGCGGCTTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                   <  1:3576090/150‑1 (MQ=255)

                                                   aTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                   >  2:3529426/1‑150 (MQ=255)

                                                    ttGCGGATTCGCTTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                  <  1:2937592/150‑1 (MQ=255)

                                                    ttGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCGCGCACGCTGGCGCACGAGTCGAACCAGGGTCGAAGCATCTCa                                                                                                  >  1:3834480/1‑150 (MQ=14)

                                                      gCGGATTCGCTTGAGAGTTAAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCCGCAAAGCGCTCTCGCTGGCGCTCCAGTAGAACCATGGTCGAAGCTACTCATc                                                                                                >  2:1483173/1‑150 (MQ=14)

                                                       cGGATTCGCTTGAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                               <  2:4276490/150‑1 (MQ=255)

                                                       cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCAAAAACGCTGGCGCTCGAGTCGAAACATGGTCGAATCCTCTCAAcc                                                                                               >  2:1486777/1‑150 (MQ=14)

                                                        ggATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCt                                                                                              >  2:114801/1‑150 (MQ=255)

                                                        ggATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCTATCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTAAACCt                                                                                              >  2:114776/1‑150 (MQ=14)

                                                        ggATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACCGCTGCGCCCTTCGCCCATCCGGTCCATGCCGCCCGCAACGCCCGCTCGCGCGCCCCCCAGGCGAACCGTCGGCGAAGCCTCTCATccg                                                                                              >  2:792398/1‑149 (MQ=9)

                                                            tCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc                                                                                          <  1:2445537/150‑1 (MQ=255)

                                                                ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGcc                                                                                      >  2:1413363/1‑149 (MQ=255)

                                                                      tttCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                <  2:1674525/149‑1 (MQ=255)

                                                                          aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCACTCGCTGGCGCTCGACTCGAACCATGGTAGAAGCTTCTCATACATACCCCCATGGGCAGaaa                                                                            >  2:820679/1‑149 (MQ=14)

                                                                           gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                                   <  1:1848676/78‑1 (MQ=34)

                                                                            ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGGCGCACGAGTCGAAACTTGGTCGAAGCATCACATCCATACCCGCAAGGGCAGAatat                                                                          >  2:3917819/1‑150 (MQ=14)

                                                                            ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAACGCTCTCTCGCTGGCGCGCGCGTCGAAACTTGGTCGAATCAGATAATCCTTCACCCCAAGGGCACAatat                                                                          >  2:779172/1‑150 (MQ=14)

                                                                                ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                      <  2:1889843/150‑1 (MQ=38)

                                                                                ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                      <  2:2606455/150‑1 (MQ=38)

                                                                                ttttGAAAGTGATGGAGGTCGGGGAAGCATTACTCAGCGCTGCGCGCTTCGCCCTTCGGCTCGTAGCCTGCGGCAACGCTCAATCGCTGGCGCTCGAGACGAACCATGGTCGAAGCATCACATCCTTCCCCGCATGAGCACAAAATTaaa                                                                      >  2:3107846/1‑150 (MQ=11)

                                                                                 tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTAAAGCGTCCCCGCCTGGGCCGGAGATTTAAt                                                                     >  2:2077714/1‑150 (MQ=14)

                                                                                    gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCAGCGGCAAACCTCTATAGCTGGCGCTCCAGTAGAACCAAGGTAGAAGCTTATAAAACATACCCCCAAGGGCCGAATATTTaaattc                                                                  >  1:3035759/1‑146 (MQ=2)

                                                                                    gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                  <  2:1338220/150‑1 (MQ=39)

                                                                                    gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                  <  1:3206026/150‑1 (MQ=39)

                                                                                     aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCAAGCGGCAAAGCAATCTCGCTGGCGCTCGATTCGAAACTTGGTCGAATCATCACATACTTACCCCCAAGGGCAGAATATTTaaattcg                                                                 >  2:3366024/1‑145 (MQ=11)

                                                                                     aaaGTGATGGTGGTGGGGGAAGGATTACACAGCGCGACGACCTTCGCCCATCGGGTCCTTGCATCCCCCAAAGCTATCTCTCTAACGAACGAGACGAACCTAGGTCCAATCTTATCATACATCCCCGCATGCGCCGAATAATTAAATGCg                                                                 >  2:2467010/1‑150 (MQ=2)

                                                                                       agcgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                                                                                               >  2:2383885/4‑102 (MQ=1)

                                                                                       agcgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                                                                                               <  1:2383885/99‑1 (MQ=1)

                                                                                           aTGGTGGTGGGGGAAGGATTACTCATCGCTGAGCGCTTCTCCCTTCTGGTCGTTTCCTGGGGCAAGGCTGTCTGGGTGGCGGTGGGGTGGATCATTGTTAGAAGCTTGTGATGCTTCCCCGGATGGGGGGTATTTTTATTTGGGGAttgg                                                           >  2:3397051/1‑148 (MQ=2)

                                                                                             tgtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCTGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                         <  2:2726518/149‑1 (MQ=255)

                                                                                             ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGATATAGCTTGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATACTTCCCCGCATGGGCAGAAAAATTAATTGCGGATTAGtt                                                         >  2:3695564/1‑150 (MQ=11)

                                                                                               tggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTgg                                                       >  2:2490702/1‑150 (MQ=38)

                                                                                               tggtggGGGAAGGATTACTCAGCGAAGCGAGATACGCCCTTCGGGTCCTTGCCTCAGGCAACGCTCTCACGAAGGAGCTCGACTCGAACCTTGGACGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTgg                                                       >  2:2490565/1‑150 (MQ=11)

                                                                                                ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTggg                                                      <  1:316745/150‑1 (MQ=39)

                                                                                                 gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                  <  1:3877447/137‑1 (MQ=21)

                                                                                                 gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                  >  2:4328407/1‑137 (MQ=37)

                                                                                                 gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                  <  1:4328407/137‑1 (MQ=37)

                                                                                                 gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                  >  2:3877447/1‑137 (MQ=17)

                                                                                                    ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCAGGGGCAGAATATTTAAGTGGGGATTGGGTGGGGAGt                                                  >  2:4035016/1‑150 (MQ=11)

                                                                                                    ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCAGGGGCAGAAGATTTAATTGCGGATTCGATGGGACGt                                                  >  2:3987044/1‑150 (MQ=14)

                                                                                                           gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTCGCGCCGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGCTTCGGTGGGGAGTTCAGGGa                                           >  2:884134/1‑150 (MQ=11)

                                                                                                           gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTATCGCAGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATACATCCCCCCATGGGCAGAATATTTAATTTCGGATTAGTTGGGAAGTTCAGGGa                                           >  2:371042/1‑150 (MQ=17)

                                                                                                                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                           <  1:1234588/96‑1 (MQ=33)

                                                                                                                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                           >  2:1234588/1‑96 (MQ=33)

                                                                                                                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTg                                     >  2:1640216/1‑150 (MQ=255)

                                                                                                                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAATCATCTCATCCTTACCCCCAAGGGCAGAATATTTAATTGCGGCTTCGTTGGGGAGGTAAGGGTCttttt                                     >  2:1222585/1‑149 (MQ=14)

                                                                                                                    gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                  <  1:3210562/150‑1 (MQ=255)

                                                                                                                      gcTGCGCGCTTCGCCCTTGGGGTAGTTGCCGGGGGCATCGGTCTGTCGCTGTCTCTCGAGTAGTACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGt                                <  1:2844488/150‑1 (MQ=14)

                                                                                                                       cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTCTTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTg                               >  2:1461840/1‑150 (MQ=255)

                                                                                                                       cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTg                                                        >  2:299013/1‑125 (MQ=39)

                                                                                                                       cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTg                                                        <  1:299013/125‑1 (MQ=39)

                                                                                                                           gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtgg                           <  1:1231512/150‑1 (MQ=255)

                                                                                                                            cgccacGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTCTCGCACGAGTCTAAACATGGTCGAAGCTTCTAATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggt                          <  1:752800/145‑1 (MQ=14)

                                                                                                                            cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggt                          <  1:2494508/150‑1 (MQ=255)

                                                                                                                             gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtg                         <  1:1651206/150‑1 (MQ=255)

                                                                                                                               ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATACATCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGTAAGTGAtggtggt                       >  2:561698/1‑150 (MQ=255)

                                                                                                                               ttCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTAGAAGCTTCTCATACATACCCGCATGGGCAGAAAATTTAATTTGGGATTCGTTGGGAAGTTCATGGACTTTTTAAAGTGAtggtggt                       >  2:2927809/1‑150 (MQ=17)

                                                                                                                                   cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                   <  1:4110435/150‑1 (MQ=255)

                                                                                                                                        cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg              <  1:3270958/150‑1 (MQ=255)

                                                                                                                                              gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCGTCTCATCCGTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa        <  2:3834480/150‑1 (MQ=255)

                                                                                                                                              gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTAGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGTCAGAATATTTTTTTGCGGATTCGTTTGTAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa        <  1:3937826/150‑1 (MQ=14)

                                                                                                                                               ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAAGATTTAATTGCGGAATCGTTGGGAACTTCAGGGACGTTTGACAGTGATGGTGGTGGGGGAAGGATTCgga       >  1:2763379/1‑148 (MQ=17)

                                                                                                                                                 gCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAcc     <  1:666265/150‑1 (MQ=255)

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

CCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286974‑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: `‑`