BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,122	T→G	10.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	10.4%	309.5 / 10.3	142	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 (64/55); new base G (14/4); total (82/60)
Fisher's exact test for biased strand distribution p-value = 7.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286973‑1287271

gCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCt                                                                                                                                                       >  1:3750758/1‑150 (MQ=255)

  cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGc                                                                                                                                                     >  2:1914808/1‑150 (MQ=255)

           gTCGAAGCTTCTCATCCTTCCCGGCTTGGGATTAATTTTTGTTTGCGTATTCGCTTTAGTGTTCTGGTACTTTTGAATGTTTTGGTGGTGGGGGATGGTTTACTCAGCGCTGCGGGCTTCGCCCTTGGGGTGGTTGCCGCTGGCAACgac                                                                                                                                            <  2:89105/150‑3 (MQ=2)

           gTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCatagca                                                                                                                                            >  2:3863113/1‑145 (MQ=255)

           gTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGct                                                                                                                                            >  1:3913579/1‑150 (MQ=255)

           gTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGca                                                                                                                                            >  1:1900710/1‑149 (MQ=255)

               aaGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGTGGCCTGCGGCAGCGctctct                                                                                                                                        >  1:2605372/1‑150 (MQ=255)

                 gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCg                                                                                                                                      >  1:3549176/1‑150 (MQ=255)

                  cTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAACCGCTCTCGc                                                                                                                                     >  2:3251748/1‑150 (MQ=255)

                   ttCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCt                                                                                                                                    <  1:58803/150‑1 (MQ=255)

                   ttCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCt                                                                                                                                    <  1:2411376/150‑1 (MQ=37)

                   ttCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCt                                                                                                                                    >  2:258579/1‑150 (MQ=255)

                    tctcATCCCTCCCCGCAAGTGCATAATATATGCTGGCTGTTTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                                                                                                                   <  2:32073/150‑1 (MQ=14)

                       cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcg                                                                                                                                <  2:4085772/150‑1 (MQ=255)

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

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCTCCCTTCGGGTCGTTTCATGCGGCAACGCTATCTCGCTGGCGATCGa                                                                                                                           >  2:2540443/1‑150 (MQ=14)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCGGCCGCAACGCTCTCTCGCTGGCGCACGa                                                                                                                           >  2:1506173/1‑150 (MQ=14)

                             tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCGGGCGCTCGAg                                                                                                                          >  1:500480/1‑150 (MQ=255)

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

                               cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTTCGCTCATCTCCCATATGGTCATAGCGTGGTGCAACGCACTCTAGGTGTCGATAGAGta                                                                                                                        >  2:3097017/1‑149 (MQ=14)

                               cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTTCGCGCATCTCCCTTCTGGTCGTTGCATGCGGCAACGCACTCTCGCTGACGCTCGAGTc                                                                                                                        >  2:3097034/1‑150 (MQ=14)

                                 cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGCGCTGAGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAACGCTCTCTCGCTGGCGCGCGAGTCGa                                                                                                                      >  2:4145539/1‑150 (MQ=255)

                                       ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTTAAAGTGATGGTGTTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg                                                                                                                <  1:1157548/150‑1 (MQ=255)

                                        gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                                                                                               >  1:1013877/1‑150 (MQ=255)

                                        gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                                                                                               >  2:2737999/1‑150 (MQ=255)

                                            aaTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGa                                                                                                           >  1:3705183/1‑150 (MQ=255)

                                            aaTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCATGGTCGa                                                                                                           >  2:1550959/1‑150 (MQ=255)

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

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

                                                tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCt                                                                                                       >  2:542861/1‑150 (MQ=255)

                                                tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAATCt                                                                                                       >  2:2224575/1‑150 (MQ=255)

                                                 ttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCGGGCGCGCGACTCGAAACGTGGTCGGAGGtt                                                                                                      >  2:1810251/1‑150 (MQ=14)

                                                     ttGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCATTGTCGAAGCTTCTCa                                                                                                  >  2:2613071/1‑150 (MQ=255)

                                                     ttGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGGCGCGCGAGTCGAACCGTGGTCGAAGCTTCTCa                                                                                                  >  1:4046415/1‑150 (MQ=255)

                                                     ttGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCCGCAACGCTATAACGCGGGCGCGCGAGTCGAACCGTTGGCGAAACTTCTCa                                                                                                  >  1:101250/1‑150 (MQ=14)

                                                     ttGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCGTCGCCCTTCGGGTCGTTGCCAGCGGCCACCCTCTCTCGCTGGCGCTCGAGTAGAACCTTGGTCGTAGCTTCTCa                                                                                                  >  1:4046545/1‑150 (MQ=14)

                                                      tGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTTGCCGGCCGCCACCCACACACCCAGGCGCGCCAGTCCAACCAGGGTCGAAGCAtctcca                                                                                                 >  1:332812/1‑148 (MQ=14)

                                                       gCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                >  2:2144739/1‑150 (MQ=255)

                                                       gCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGAGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                >  2:2144909/1‑150 (MQ=255)

                                                       gCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                <  2:1650660/150‑1 (MQ=255)

                                                                cTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCAGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCATCCCCga                                                                                       >  1:3919707/1‑149 (MQ=255)

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

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCAAGGTCGAAGCTTGACAACCATCCCCCCa                                                                                      >  2:3496880/1‑150 (MQ=14)

                                                                    agagTTCAGGGACTTTTGAATGTGATGGTGGTGGTGGAAGTATGACTCCGCGCTGCGCGCTTGGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                   <  1:2472891/150‑1 (MQ=14)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGGCCAAGCTACTCATCCTTCCCCGCATGGGcc                                                                                >  2:2404426/1‑149 (MQ=14)

                                                                          cAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                             <  1:1814007/150‑1 (MQ=40)

                                                                          cAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                             <  1:1813980/150‑1 (MQ=40)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCATCTCATCCATCCCCGCCAGGGCAGAat                                                                            >  2:4007184/1‑150 (MQ=255)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACACTAGCAGGCGCGCGAGTCGAACATTGGTCGAATCTTCTCAAACTTACCCGCAAGGGCAGaaa                                                                            >  1:1188514/1‑149 (MQ=14)

                                                                                cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcg                                                                                                                                <  1:179302/93‑1 (MQ=33)

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

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCAGTCGCACGACTCGAACCATGGTCGAAGCTTCTAAAACATACCCCCATGGGCAGAaaataaaa                                                                      >  1:230299/1‑146 (MQ=11)

                                                                                      aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCg                                                                 <  1:2177713/150‑1 (MQ=21)

                                                                                      aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCTGGCGCTAGAGTCGAAACTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCACAATATTTAATTGCg                                                                 >  2:1680945/1‑150 (MQ=255)

                                                                                           gagggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                           <  1:3017470/116‑1 (MQ=17)

                                                                                           gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                           >  2:3017470/1‑119 (MQ=32)

                                                                                               gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTg                                                        >  2:2329320/1‑150 (MQ=38)

                                                                                               gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAAAGCACAAACGCTGGCGCTCGAGTAGAAACTTGGTCGAAACTTCTCATCCATACCCCCATGGGCCGAATATGTAAATGAGGCttagatg                                                        >  2:2806/1‑145 (MQ=11)

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

                                                                                                tggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGACTGCGGCAACGCTCTCTCGCTGGCCCTCGAGTCGAACCTTGGTAGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATGTAAATGCGGAGGCGTTgg                                                       >  2:3067419/1‑150 (MQ=11)

                                                                                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                  <  1:2521407/105‑1 (MQ=32)

                                                                                                   tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTTTCCTGCGGCAACGCACTCTCGCAGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCATCCCCCCATGGGCAGAATATTTAATTTCGTAATACTTGGGaa                                                    >  2:3623468/1‑150 (MQ=11)

                                                                                                      ccgAACGAGTACTCCGCGCTCCGCGCTCCTCCCTTCTGCTCGTTGAGTGCGGCGAGGCTATATCTGTGCAGGTCGGCTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGtt                                                 <  1:1894260/148‑1 (MQ=9)

                                                                                                       ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCACACGCTGGCGCTCGAGTCGAACCGTGGTCGAATCGTCAAATACATACCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTc                                                >  2:1443299/1‑150 (MQ=14)

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

                                                                                                          aGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCTGCAACGCTCTCTCGCTGGCGATCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                             <  2:3750758/150‑1 (MQ=255)

                                                                                                          aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                             <  2:2135617/150‑1 (MQ=255)

                                                                                                              ttACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACt                                         <  1:717320/150‑1 (MQ=255)

                                                                                                               tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtt                                        <  1:258579/150‑1 (MQ=255)

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

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

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGGTGGGGAGTTCAGGGACTTTTg                                     >  2:2383527/1‑150 (MQ=11)

                                                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                   >  2:4154845/1‑150 (MQ=255)

                                                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                   <  1:2144739/150‑1 (MQ=255)

                                                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCACGAGTCGAACCTTGGTCGAAGCTTCTCAACCGTACCCGCAGGGGCACAAAATTTAATAGCGGCTTCGATTGGACGATCAGGGAATTTTGaa                                   >  2:3985723/1‑150 (MQ=14)

                                                                                                                       gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAAGATTTAATTGCGGATTCGTTGGGAAGGTCAGGGACTGTTGATAGt                                >  2:3588088/1‑150 (MQ=17)

                                                                                                                           cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                     <  1:4119341/61‑1 (MQ=35)

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

                                                                                                                                   gCCCGTCGGGTCGTGGCCGGCGGCAACGCTCTCTCGCTGTCGCTCGAGACGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                    <  2:2955350/150‑1 (MQ=14)

                                                                                                                                     ccTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAAACTTCCCCCCCTGGGCCGAATATTTAATTGCGGTGTCGGTGGGAAGTTGAGGGCCGTTTGTAAGTGAgggggggggggg                  >  2:1076327/1‑150 (MQ=14)

                                                                                                                                      cTTCGGGTCGTTGGCTGCGGCATAGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGGTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTTGTTAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGa                 <  2:3373551/150‑1 (MQ=11)

                                                                                                                                      cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGa                 >  2:677278/1‑150 (MQ=255)

                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                     >  2:1616173/1‑65 (MQ=35)

                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                     <  1:1616173/65‑1 (MQ=35)

                                                                                                                                         gagagcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg              >  2:1163888/6‑150 (MQ=255)

                                                                                                                                           ggTCGTTGCCTGCGGCAACCCTCTATCGCTGGCGCTCGAGTCGAACCATGGTCGAAGCTTCTCATACATACCCCCATGGGCAGAATATTTAATTGCGGATTAGTTTGGAAGTTCAGGGACGTTTGAAAGTGATTGTGGTGGGGGAAGGAt            >  2:868276/1‑150 (MQ=14)

                                                                                                                                            ctCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt           <  2:2554380/149‑1 (MQ=255)

                                                                                                                                             tCGTTGCCTGCGGCACCGCTCTCTCGGTGCCGCTCGAGTGTTACCTTTGTCGACGCTTGTGCTCCTCCGCCGCTTGGGCAGAATATTTTATTGCGGATTCGTTGGGTATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc          <  1:3386041/150‑1 (MQ=11)

                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGa                                                                                                           <  1:3401878/53‑1 (MQ=35)

                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCACGAGTCGAACCTTGGTCGATGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGCAGTTCAGGGCCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc          <  1:2425050/150‑1 (MQ=14)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa        <  1:961167/150‑1 (MQ=255)

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

                                                                                                                                                 tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTTGTAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAc      <  1:936771/150‑1 (MQ=255)

                                                                                                                                                  gCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGACGGATTCGAAcc     >  2:3245041/1‑150 (MQ=255)

                                                                                                                                                   ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTATCATCCATACCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGGAGGATTAGAACCt    >  1:2782739/1‑150 (MQ=17)

                                                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGACCCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCtt   <  2:332812/150‑1 (MQ=255)

                                                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTAATTGCGGATTCGTTGGGACGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGATTCGAACCtt   >  2:1705502/1‑150 (MQ=255)

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

                                                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCCCATGGGCAGAATATTTAATTTCGGATTCGTTGGGTAGTTCAGGGACTTTTGTAAGTGATTGTGGTGGGGGAAGGATTAGAACCtt   >  1:2668335/1‑150 (MQ=17)

                                                                                                                                                     tGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                         <  2:3524412/95‑1 (MQ=18)

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