BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,122	T→G	7.5%	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	7.5%	523.5 / 10.9	224	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 (92/100); new base G (20/2); total (121/103)
Fisher's exact test for biased strand distribution p-value = 7.80e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 228 total aligned reads displayed.

CCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286974‑1287271

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

   aaCCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                                   <  2:1650402/150‑1 (MQ=255)

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

       ttAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAc                                                                                                                                               >  1:4034727/1‑150 (MQ=255)

          ttCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAGTCAGCGCTGCGCGCTACGCCCTGCGGGTCGTTGCCTGCGGGAGCGct                                                                                                                                            >  2:2269003/2‑150 (MQ=14)

             gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGTAAGTGTTTGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCGGGTCGTTTGCTGAGGCAActcaatc                                                                                                                                         >  1:1243038/1‑144 (MQ=14)

             gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGACAGCTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACACAGCTCTGCGCCCATCACCCTTCGGGTCCATGCCTGAGGCAAAGccctc                                                                                                                                         >  2:5110409/1‑150 (MQ=14)

                   tctcATCCTTCCCCGCTTGGGCATAATTTTTGATTGCGGATTCGCTTGAGATTTCAGGGACTTTTGATAGTGATGGTGGTGGGGTAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                                                                                                                   <  1:6375297/150‑1 (MQ=11)

                    ctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATAAGCGCGTCGCGCATCGCCCAACGGGTCGATGCCGGCGGCAACGCACGCTCGGTgg                                                                                                                                  >  2:1139948/1‑150 (MQ=14)

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

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

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

                     tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGAGCTTCGCGCTTCTCCCATCGGTTAGTTTCGTGCCCCAAAGCACTCTCgagttc                                                                                                                                 >  2:6154092/1‑145 (MQ=14)

                       aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCTCTGAGAGCTTCGCCCATCGGTTCCTGGCGTGACGCAAAGCGCTCTCGcggcccc                                                                                                                               >  2:1172178/1‑144 (MQ=14)

                           ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCATCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                           <  1:5496547/150‑1 (MQ=255)

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

                           ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGGCGCTCGa                                                                                                                           >  1:4437704/1‑150 (MQ=255)

                            tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCAGGCGCTCGAg                                                                                                                          >  2:3076607/1‑150 (MQ=255)

                            tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTTCGCCCTTCGGGTCGTTGACGCAGGCAACCCTCTCTCGCTCGCGCGCGAg                                                                                                                          >  2:411582/1‑150 (MQ=14)

                               ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCGCACTCGCTGGCGCGCGAGTCg                                                                                                                       >  2:5481110/1‑150 (MQ=255)

                               ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAACCTCTGCGACCATCTCCCTTCTGGTCGATGCCTCAGGCAAAGGACTCTCGCTATCGCACGAATCg                                                                                                                       >  2:5507832/1‑150 (MQ=14)

                                 gCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                     <  1:3148464/150‑1 (MQ=255)

                                  cttGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAAc                                                                                                                    <  2:2044804/148‑1 (MQ=255)

                                   aTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAAcc                                                                                                                   >  2:5696125/1‑150 (MQ=255)

                                   aTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCACTCCCTGGCGCTCGAGTCGAAcc                                                                                                                   >  1:3330062/1‑150 (MQ=255)

                                      ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg                                                                                                                <  2:5443952/150‑1 (MQ=255)

                                       gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCACTCTCGCTGGCGCACGAGTCGTACCGTgg                                                                                                               >  2:5151559/1‑150 (MQ=14)

                                       gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTACGCCCTTCGGGTCGTTGCCGGCGGCAACCCACTCACGCACGCGCTCGAATCCAAccaggg                                                                                                               >  1:2406939/1‑146 (MQ=14)

                                         aGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACACTCGCTGGCGCTCGAGTCGAACCTTGGTc                                                                                                             >  1:6669294/1‑150 (MQ=255)

                                             tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTATCGCTGTCGCTCGACTCGAACCGTGGTCGAAg                                                                                                         >  1:3690757/1‑150 (MQ=14)

                                               tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCCGCAACGCACTCACGCTGGCGCTCGAGTCGAACCATGGTCGAAGCt                                                                                                       >  2:64872/1‑150 (MQ=255)

                                                   aTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACGCACGCGCGCAGGCGCACCAATAGAAACATAGTCGAAGCTtcta                                                                                                   >  1:1206226/1‑149 (MQ=14)

                                                    ttGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAAAGCTCTCTCGCTGCCGCGCGAGTCGAAACTTGGTAGAAGCTTCTCa                                                                                                  >  2:1285302/1‑150 (MQ=14)

                                                    ttGCGGATTAGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                  <  2:1875585/150‑1 (MQ=255)

                                                      gCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCGCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCATCACATc                                                                                                >  1:4894663/1‑150 (MQ=14)

                                                         gATTCGGTTGAGAATTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGGCGCTCGAGGCGAACCTTGGTCGAAGCATCTCATccat                                                                                             >  2:2444705/1‑148 (MQ=11)

                                                           ttCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCACGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCatccatcc                                                                                           >  2:204857/1‑150 (MQ=255)

                                                           ttCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCACACTCGCTGGCGCTCGAGTCGAACCATGGTCGAAGCATCTCAtacaacc                                                                                           >  1:1599891/1‑144 (MQ=14)

                                                              gCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCAGCGGCAAAGCTCTCACGCTAGCGCTAGAGTCGAAACTTGGTCGAAGCTTCTCATCAATCCCCg                                                                                        >  2:2310845/1‑150 (MQ=14)

                                                              gCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCGCGCTGGCGCTCGACGCGAACCTTGGTAGAAGCTTCTCAACCATccccc                                                                                        >  2:6242451/1‑149 (MQ=14)

                                                                ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                      <  1:988626/150‑1 (MQ=255)

                                                                ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                      <  1:6294909/150‑1 (MQ=255)

                                                                    gaattCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                  <  1:3908065/147‑1 (MQ=21)

                                                                       ttCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCACGCTGGCGCTCCAGTCGAACCGTGGTCGAAGCATCTCATCCAACCCCGCAAGGGCAg                                                                               >  2:1352368/1‑150 (MQ=14)

                                                                           gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCTACGCTCTCTCGCTGGCGCTCGAGTCGAAACTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                           >  2:3610373/1‑150 (MQ=255)

                                                                            ggtcTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAatat                                                                          <  1:2444697/147‑1 (MQ=37)

                                                                            ggccTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAatat                                                                          <  1:5765429/147‑1 (MQ=37)

                                                                             gACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACACTAGCTGGCGCACGAGTCGAACCATTGTCGAAGCTTCACAAACATCACCCCAAGGGCAGAataat                                                                         >  1:488185/1‑148 (MQ=11)

                                                                             gACTATTGAATGTTTTGGTGGTGGGGTAAGTATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                         <  1:59494/150‑1 (MQ=11)

                                                                                ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCATCTCATCCTTCCCCCCATGGGCAGAATAATTaa                                                                      >  2:4153302/1‑150 (MQ=255)

                                                                                     aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCg                                                                 <  1:5936591/150‑1 (MQ=39)

                                                                                       aGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGa                                                               <  1:4763702/150‑1 (MQ=39)

                                                                                             ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                         <  2:4034709/150‑1 (MQ=37)

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

                                                                                                   gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAg                                                   >  2:3728890/1‑150 (MQ=255)

                                                                                                    ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGt                                                  >  2:3535999/1‑150 (MQ=255)

                                                                                                     gggAAGGATTACTCAGCGCTGCGCGCTTCGCACTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGtt                                                 <  2:4324171/150‑1 (MQ=255)

                                                                                                       gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCATGGTCGAAGCTTCTCATCCTTCCCCGCAAGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGttaa                                               >  2:682322/1‑148 (MQ=255)

                                                                                                        aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGTGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAg                                                                                                         <  1:5590271/91‑1 (MQ=18)

                                                                                                         ccgatctcTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                             <  1:2705585/143‑1 (MQ=255)

                                                                                                         aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                             <  1:3845058/150‑1 (MQ=255)

                                                                                                          ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCtt                                                                                                      >  2:488254/1‑92 (MQ=33)

                                                                                                                cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtttt                                      <  2:4820353/150‑1 (MQ=255)

                                                                                                                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTTATTGCGGATTCGGTGGGGAGTTCAggggggtcgg                                     >  2:2705585/1‑143 (MQ=11)

                                                                                                                   aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTTGTCGAAGCTTCTCAAACTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGGTGGGAAGTTCAGGGACTTTTGaa                                   >  2:3224907/1‑150 (MQ=14)

                                                                                                                    gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                  <  1:699057/150‑1 (MQ=255)

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

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

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

                                                                                                                    gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCTATCGCTGGCGCTAGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                  <  1:1753730/150‑1 (MQ=255)

                                                                                                                     cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAg                                 <  1:1424292/150‑1 (MQ=255)

                                                                                                                      gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGt                                <  1:4101148/150‑1 (MQ=255)

                                                                                                                        tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                              <  1:4466097/150‑1 (MQ=255)

                                                                                                                        tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                              >  2:3591417/1‑150 (MQ=255)

                                                                                                                        tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                              <  2:5924557/150‑1 (MQ=255)

                                                                                                                        tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCAAGCTGGCGCTCGGGTCGAACCGTGGTCGAAGCATCTCAAACTTCCCCCCCTGGGCAGAAAATTTAATTTCGGGTTCGTTGGGAAGTTGAGGGGAGTTTGGAAGTGa                              >  2:4508443/1‑150 (MQ=14)

                                                                                                                          cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtg                            <  1:6499501/150‑1 (MQ=255)

                                                                                                                           gcgcTTCGCCCTTCGGGTGGTTGCCTGCGGCTCCTCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTTGTAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtgg                           <  2:5260392/150‑1 (MQ=17)

                                                                                                                                 cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                     <  1:2191681/150‑1 (MQ=255)

                                                                                                                                 cGCCCTTCGGGTCGTTGCCTGCGGATACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                     <  1:893497/150‑1 (MQ=255)

                                                                                                                                  gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTATCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCATCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACGTTTGAAAGTGTTGGTgggggg                    >  2:1300315/1‑150 (MQ=17)

                                                                                                                                   cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTACCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                   <  1:5135170/150‑1 (MQ=255)

                                                                                                                                    ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGCCTTTTGAAAGTGATGGTGGTggggg                  >  2:5410457/1‑150 (MQ=255)

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

                                                                                                                                       tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg               <  1:343094/150‑1 (MQ=255)

                                                                                                                                         gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa             <  2:5830419/150‑1 (MQ=255)

                                                                                                                                         agagcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGACGGa             >  2:1739209/5‑150 (MQ=17)

                                                                                                                                          ggTCGTTGCCTGAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGTAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt            <  1:5957978/150‑1 (MQ=255)

                                                                                                                                          gagcGTTGCCGCAGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAt                             <  2:6101046/130‑1 (MQ=14)

                                                                                                                                           gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTACCCGCATGGGCAGAATATTTAATTGAGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt           >  2:3803005/1‑150 (MQ=255)

                                                                                                                                                tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGATCCTTGGTCGTAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGTTAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAc      <  1:5682897/150‑1 (MQ=14)

                                                                                                                                                 gCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAATTCTGAATCGTTGGGAAGTTCAGGGACTGTTGAAATTGGTGGGGGTGGGGGGAGGGTTCGAAcc     >  1:288898/1‑150 (MQ=17)

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

                                                                                                                                                  ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGGTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCt    >  2:3159600/1‑150 (MQ=255)

                                                                                                                                                  ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTAGAACCTTGGTCGAAGCGTCTAATCCTTACCCGCAGGGGCCGAATATTTAATTGCGGATTCGTTGGGGAGTTCAGGGACTTTTGTAAGTGATGGTGGTGGGGGGAGGATTCGAACCt    >  1:4682511/1‑150 (MQ=17)

                                                                                                                                                   cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCtt   <  2:776196/150‑1 (MQ=255)

                                                                                                                                                   cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCATGGTCGAAGCTTCTCAACCTTCCCCCCAGGGGCAGAATATTTAATTGCGGGTTCGTTGGGGAGTTGAGGGACATTTGAAAGTGATGGTGGTGGGGGGAGGATTCGAAccgt   >  1:3049981/1‑148 (MQ=17)

                                                                                                                                                   cTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                              <  1:967317/107‑1 (MQ=255)

                                                                                                                                                    tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                   <  1:5333219/53‑1 (MQ=35)

                                                                                                                                                    tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  >  2:4676176/1‑150 (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: `‑`