BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	450,921	T→C	6.6%	noncoding (115/127 nt)	PP_mr05 →	RsmY ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	450,921	0	T	C	6.6%	661.2 / 24.1	212	noncoding (115/127 nt)	PP_mr05	RsmY ncRNA
Reads supporting (aligned to +/- strand): ref base T (73/122); new base C (3/11); total (77/135)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

Only 100 of 214 total aligned reads displayed.

TGTTCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCAGGC  >  NC_002947/450772‑451070

   tCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcg                                                                                                                                                    <  2:2866868/150‑1 (MQ=255)

     tAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggc                                                                                                                                                  >  2:599946/1‑150 (MQ=255)

      aaGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggcg                                                                                                                                                 <  2:2856129/150‑1 (MQ=255)

      aaGATAGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggcg                                                                                                                                                 >  2:4133245/1‑150 (MQ=255)

       aGATTGGATCCACTGTGAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggcgg                                                                                                                                                <  2:3036895/150‑1 (MQ=255)

        gATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCggg                                                                                                                                               <  2:1544472/150‑1 (MQ=255)

            ggATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGttt                                                                                                                                           <  2:1921676/150‑1 (MQ=255)

            ggATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGttt                                                                                                                                           <  1:1879205/150‑1 (MQ=255)

             gATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGtttt                                                                                                                                          <  2:4471317/150‑1 (MQ=255)

              aTCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGTTTTc                                                                                                                                         <  2:174219/150‑1 (MQ=255)

                ccACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCtt                                                                                                                                       <  1:1750350/150‑1 (MQ=255)

                 cACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                                      <  2:2899716/150‑1 (MQ=255)

                 cACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                                      >  2:1317493/1‑150 (MQ=255)

                  aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtg                                                                                                                                     <  1:1040228/150‑1 (MQ=255)

                  aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTAGAAAACCCCGCCGAAGCGGGGTTTTccaga                                                                                                                                     <  2:4063268/150‑5 (MQ=255)

                  aCTGTTAATCTTTACCACATGGACGTAGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACTAAGGACCACCTGACAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCACGCTTCGGCGGGGTTTTCTTtg                                                                                                                                     >  1:2785603/1‑150 (MQ=255)

                    tGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcg                                                                                                                                   >  2:1239064/1‑150 (MQ=255)

                    tGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGTTTTccagact                                                                                                                                   >  2:2991655/1‑144 (MQ=255)

                     gtgaaTCCTTACCACAGGGACGTCGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGACACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgt                                                                                                                                  <  2:829137/147‑1 (MQ=255)

                       tAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgtgc                                                                                                                                <  1:2022657/150‑1 (MQ=255)

                             ttACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGttt                                                                                                                          >  1:2364217/1‑150 (MQ=255)

                                 cacaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                                      <  1:3050118/150‑1 (MQ=255)

                                  acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                                     <  2:151301/150‑1 (MQ=255)

                                    aTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGc                                                                                                                   <  1:4353991/150‑1 (MQ=255)

                                     tGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGcc                                                                                                                  <  1:238878/150‑1 (MQ=255)

                                      ggACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCa                                                                                                                 <  2:2431072/150‑1 (MQ=255)

                                      ggACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCa                                                                                                                 <  1:3723595/150‑1 (MQ=255)

                                       gACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAg                                                                                                                <  1:522717/150‑1 (MQ=255)

                                       gACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAg                                                                                                                <  2:4391815/150‑1 (MQ=255)

                                             gcgcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATc                                                                                                          <  1:977785/150‑1 (MQ=255)

                                               gcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCa                                                                                                        >  1:2052579/1‑150 (MQ=255)

                                                      gcgcAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGa                                                                                                 <  2:3959539/150‑1 (MQ=255)

                                                         cAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGc                                                                                              >  2:3718366/1‑150 (MQ=255)

                                                          aaaGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCt                                                                                             <  2:2876723/150‑1 (MQ=255)

                                                          aaaGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCt                                                                                             >  1:4334646/1‑150 (MQ=255)

                                                          aaaGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCGTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCt                                                                                             >  2:3026530/1‑150 (MQ=255)

                                                              gACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCa                                                                                         >  2:1815999/1‑150 (MQ=255)

                                                                cAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATc                                                                                       <  2:3108584/150‑1 (MQ=255)

                                                                 aGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGTCAGCACATCCAGCGGTGATGCTCCCATCa                                                                                      <  1:4075161/150‑1 (MQ=255)

                                                                      aaGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGa                                                                                 <  1:4170062/150‑1 (MQ=255)

                                                                      aaGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGa                                                                                 <  1:812277/150‑1 (MQ=255)

                                                                      aaGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGa                                                                                 <  1:3194733/150‑1 (MQ=255)

                                                                       aGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAg                                                                                >  2:114127/1‑150 (MQ=255)

                                                                          acgacgacgAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgt                                                                             <  2:1751689/150‑1 (MQ=255)

                                                                            gacgacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtcc                                                                           >  1:3960255/1‑150 (MQ=255)

                                                                               gacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtccgtc                                                                        <  1:65000/150‑1 (MQ=255)

                                                                               gacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtccgtc                                                                        <  1:3517892/150‑1 (MQ=255)

                                                                                 cgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGtctc                                                                      <  2:2254457/150‑1 (MQ=255)

                                                                                  gaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTGTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTctct                                                                     >  2:3432110/1‑150 (MQ=255)

                                                                                  gaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGCTGTTCTTTGCGTGCGTGTTTTACAGGCAAGCACATCCAGAGGAGATGCTCCCATCAGGCGAGCGTCAGTctct                                                                     >  2:4472991/1‑150 (MQ=255)

                                                                                    aGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtc                                                                   <  1:4040739/150‑1 (MQ=255)

                                                                                     ggACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcc                                                                  >  1:3114228/1‑150 (MQ=255)

                                                                                     ggACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcc                                                                  <  1:4230895/150‑1 (MQ=255)

                                                                                      gACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcct                                                                 <  1:142496/150‑1 (MQ=255)

                                                                                         caccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCa                                                              >  1:1113714/1‑150 (MQ=255)

                                                                                         caccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCa                                                              >  1:2424206/1‑150 (MQ=255)

                                                                                         caccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCa                                                              <  1:708140/150‑1 (MQ=255)

                                                                                             tGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCaa                                                          <  1:1440708/150‑1 (MQ=255)

                                                                                              gCCAGGATGGCGGGGCGATACGGATATCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAAc                                                         <  1:3765135/150‑1 (MQ=255)

                                                                                                 aGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTg                                                      <  1:2297473/150‑1 (MQ=255)

                                                                                                   gATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTgcg                                                    <  1:3239323/150‑1 (MQ=255)

                                                                                                   gATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTgcg                                                    <  1:3621512/150‑1 (MQ=255)

                                                                                                   gATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTgcg                                                    <  1:1492106/150‑1 (MQ=255)

                                                                                                         ggggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                                      <  2:854251/62‑1 (MQ=255)

                                                                                                         ggggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGcagcag                                              >  1:1184457/1‑150 (MQ=255)

                                                                                                          gggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCagcagc                                             <  2:3457762/150‑1 (MQ=255)

                                                                                                                tACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCcaca                                       <  1:1552068/150‑1 (MQ=255)

                                                                                                                 aCGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACag                                      >  1:258025/1‑150 (MQ=255)

                                                                                                                   ggATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACagag                                    >  1:4467877/1‑150 (MQ=255)

                                                                                                                   ggATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACagag                                    <  2:1058509/150‑1 (MQ=255)

                                                                                                                     aTGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGc                                  <  1:3851843/150‑1 (MQ=255)

                                                                                                                        tcgcaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctg                               <  2:2194744/147‑1 (MQ=255)

                                                                                                                         cacaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGTGTTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgc                              <  2:2559663/150‑1 (MQ=255)

                                                                                                                         cacaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgc                              >  2:730311/1‑150 (MQ=255)

                                                                                                                          acaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgct                             >  2:2172247/1‑150 (MQ=255)

                                                                                                                            aGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTgg                           <  1:1827560/150‑1 (MQ=255)

                                                                                                                            aGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTgg                           >  1:2201312/1‑150 (MQ=255)

                                                                                                                             gggAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgc                          >  2:1526801/1‑150 (MQ=255)

                                                                                                                              ggAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgcg                         <  1:732340/150‑1 (MQ=255)

                                                                                                                               gAAACAGTCTGGAACACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCCGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCACCGTGCGCAGCAGCTCCCCAGAGGCctgctgttcg                        <  2:4104288/150‑5 (MQ=255)

                                                                                                                               gAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgcgc                        >  2:3941643/1‑150 (MQ=255)

                                                                                                                               gAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgcgc                        >  2:2026133/1‑150 (MQ=255)

                                                                                                                                 aaCAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTg                      <  2:1431779/150‑1 (MQ=255)

                                                                                                                                  aCAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg                     <  1:579046/150‑1 (MQ=255)

                                                                                                                                      tCTGGAAATCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCAt                 >  2:283350/1‑150 (MQ=255)

                                                                                                                                      tCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCAt                 >  1:68735/1‑150 (MQ=255)

                                                                                                                                       cTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATc                <  2:694389/150‑1 (MQ=255)

                                                                                                                                        tGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcg               >  2:3577928/1‑150 (MQ=255)

                                                                                                                                        tGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcg               <  2:2364217/150‑1 (MQ=255)

                                                                                                                                         agAAAACCCCGCCGAAGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcgc              >  1:3788138/2‑150 (MQ=255)

                                                                                                                                          gAAAACCCCGCCGAAGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcgcg             >  1:3450453/1‑150 (MQ=255)

                                                                                                                                           aaaaCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGa            <  2:3155729/150‑1 (MQ=255)

                                                                                                                                             aaCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGaaa          <  2:1387278/150‑1 (MQ=255)

                                                                                                                                                cccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATAc       <  1:2649546/150‑1 (MQ=255)

                                                                                                                                                 ccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGaaa          >  1:2652350/1‑146 (MQ=255)

                                                                                                                                                 ccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATAcc      >  1:2848079/1‑150 (MQ=255)

                                                                                                                                                 ccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATAcc      >  1:2542478/1‑150 (MQ=255)

                                                                                                                                                 ccGCGTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGaaa          <  2:2652350/146‑1 (MQ=255)

                                                                                                                                                  cGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGt                                                       >  2:4221811/1‑100 (MQ=255)

                                                                                                                                                     ttCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCAGGc  <  1:2306060/150‑1 (MQ=255)

TGTTCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCAGGC  >  NC_002947/450772‑451070

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