BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,375,490	T→G	8.3%	R30R (CGA→CGC)	PP_5637 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,375,490	0	T	G	8.3%	330.2 / 19.3	129	R30R (CGA→CGC)	PP_5637	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (65/52); new base G (6/5); total (71/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 131 total aligned reads displayed.

ACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGTATAAAGA  >  NC_002947/4375342‑4375638

acGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTc                                                                                                                                                     >  2:1687573/1‑150 (MQ=255)

      gcgcATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCa                                                                                                                                               <  2:274148/150‑1 (MQ=255)

      gcgcATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCa                                                                                                                                               <  2:395078/150‑1 (MQ=255)

         cATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCagag                                                                                                                                            <  1:1950686/150‑1 (MQ=255)

         cATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGAGTCGCGACCCaggc                                                                                                                                            <  1:751991/150‑3 (MQ=255)

          atatAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCAAGGCGTAGGGTTAAACACCCAGGCCTAGGTCGAGACTCCACCCagatt                                                                                                                                           <  2:90898/150‑3 (MQ=255)

                    ccTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGTCTCGACCCAGAGGCTCGCTGTGa                                                                                                                                 <  2:1184011/150‑1 (MQ=255)

                      ttCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATg                                                                                                                               >  2:2144472/1‑150 (MQ=255)

                            ctctTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGcc                                                                                                                         >  2:817475/1‑150 (MQ=255)

                             tcttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCa                                                                                                                        >  1:2095546/1‑150 (MQ=255)

                             tcttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGAACTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCa                                                                                                                        >  1:1021816/1‑150 (MQ=255)

                             tctgctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCa                                                                                                                        >  2:1521709/1‑150 (MQ=255)

                              cttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCaa                                                                                                                       <  2:2095546/150‑1 (MQ=255)

                               ttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                                      >  1:924028/1‑150 (MQ=255)

                               ttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                                      <  2:450545/150‑1 (MQ=255)

                                 ctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAtgt                                                                                                                    <  1:815152/150‑1 (MQ=255)

                                   aCCCCGTTGATGTTGAGGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTgg                                                                                                                  <  2:1783662/150‑1 (MQ=255)

                                     cccGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGAc                                                                                                                >  1:209848/1‑150 (MQ=255)

                                       cGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGc                                                                                                                                          >  2:182591/1‑122 (MQ=255)

                                        gttgatgttgatGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTc                                                                                                             >  1:1828347/1‑150 (MQ=255)

                                         ttgatgttgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTcc                                                                                                            >  2:69474/1‑150 (MQ=255)

                                         ttgatgttgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTcc                                                                                                            <  1:1521709/150‑1 (MQ=255)

                                           gatgttgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgc                                                                                                          >  1:236096/1‑150 (MQ=255)

                                             tgttgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgc                                                                                                        <  1:2444207/150‑1 (MQ=255)

                                             tgttgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCATAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgc                                                                                                        >  1:2224433/1‑150 (MQ=255)

                                               ttgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCgg                                                                                                      >  1:2213990/1‑150 (MQ=255)

                                                tgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGa                                                                                                     <  2:1828347/150‑1 (MQ=255)

                                                tgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCCACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGa                                                                                                     <  2:2433047/150‑1 (MQ=255)

                                                   tgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAggg                                                                                                  >  1:1152174/1‑150 (MQ=255)

                                                      cGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGACGCTCACTTTGATGGTCGCCAAGGTGGACGTCCGCGCGGAGGGGcg                                                                                               >  2:631829/1‑149 (MQ=255)

                                                       gAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTg                                                                                              <  1:2101144/150‑1 (MQ=255)

                                                         gCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGa                                                                                            <  2:1470341/150‑1 (MQ=255)

                                                         gCTAGGCCATGCCGTTGACCTGCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGTGCGGAGGGGTTGGa                                                                                            <  2:369758/150‑1 (MQ=255)

                                                               ccATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTcc                                                                                      >  2:2070764/1‑150 (MQ=255)

                                                                cATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTATGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGagatcgga                                                                                     >  1:1262448/1‑143 (MQ=255)

                                                                 aTGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTg                                                                                              <  1:2239194/140‑1 (MQ=255)

                                                                 aTGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTg                                                                                              >  2:2239194/1‑140 (MQ=255)

                                                                  tGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTa                                                                                   >  2:2050185/1‑150 (MQ=255)

                                                                  tGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTa                                                                                   >  2:1846848/1‑150 (MQ=255)

                                                                   gCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACACGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAg                                                                                  <  2:1370145/150‑1 (MQ=255)

                                                                       ttGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCgg                                                                              >  1:1712809/1‑150 (MQ=255)

                                                                           ccTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCCATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTg                                                                          >  2:697470/1‑150 (MQ=255)

                                                                            cTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGa                                                                         >  1:1105266/1‑150 (MQ=255)

                                                                            cTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGa                                                                         >  1:1620070/1‑150 (MQ=255)

                                                                            cTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGa                                                                         >  2:2249487/1‑150 (MQ=255)

                                                                             ttCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGATGGGCTGGACTTTCCGTAGGTGGTCTGAg                                                                        >  1:1631273/1‑150 (MQ=255)

                                                                             ttCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAg                                                                        <  1:2249487/150‑1 (MQ=255)

                                                                              tCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                                      <  2:451310/103‑1 (MQ=255)

                                                                              tCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                                      >  1:451310/1‑103 (MQ=255)

                                                                              tCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAgg                                                                       >  1:1793287/1‑150 (MQ=255)

                                                                                   tCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAg                                                                  <  2:924028/150‑1 (MQ=255)

                                                                                    cGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAgg                                                                 <  2:236096/150‑1 (MQ=255)

                                                                                      tGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGcc                                                               >  2:2088039/1‑150 (MQ=255)

                                                                                        cAGCCGGGGGATAATGTCGGGCCCAGCCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCAc                                                             <  2:1072423/150‑1 (MQ=255)

                                                                                             gggggATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACGCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGGCTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCAt                                                        >  2:3662/1‑150 (MQ=255)

                                                                                             gggggATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTGGGCGGTCTGAGGGCAAGGCCACGCCAt                                                        <  1:680495/150‑1 (MQ=255)

                                                                                                ggATAAGGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGGCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCt                                                     <  2:835318/150‑1 (MQ=255)

                                                                                                   tAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTcgc                                                  <  1:1846848/150‑1 (MQ=255)

                                                                                                       gTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCgcggcg                                              >  1:1938299/1‑150 (MQ=255)

                                                                                                        tCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGc                                                                                >  1:1848999/1‑115 (MQ=255)

                                                                                                            gCCCAGGCGTAGGGTTAAACACACCGCAGGTGGGCGCGGATCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGaa                                         <  1:5/150‑1 (MQ=255)

                                                                                                              ccAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACg                                       >  1:1132119/1‑150 (MQ=255)

                                                                                                              ccAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACg                                       >  2:503537/1‑150 (MQ=255)

                                                                                                               cAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGGCTGAGGGCAAAGCAATGCCATGCTCGCGGCGTCGAACga                                      >  1:602/1‑150 (MQ=255)

                                                                                                                aGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACgag                                     >  2:1790703/1‑150 (MQ=255)

                                                                                                                aGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACgag                                     <  2:1938299/150‑1 (MQ=255)

                                                                                                                       gggTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGa                                                                                                                                 <  2:2230236/51‑1 (MQ=255)

                                                                                                                       gggTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGa                                                                                                                                 >  1:2230236/1‑51 (MQ=255)

                                                                                                                         gTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTa                            <  2:209848/150‑1 (MQ=255)

                                                                                                                             aaCACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCt                                                                                               <  2:1742696/79‑1 (MQ=255)

                                                                                                                             aaCACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCt                                                                                               >  1:1742696/1‑79 (MQ=255)

                                                                                                                             aaCACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATg                        >  1:1280111/1‑150 (MQ=255)

                                                                                                                              acacCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGc                       <  1:477584/150‑1 (MQ=255)

                                                                                                                              acacCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGc                       <  2:1599589/150‑1 (MQ=255)

                                                                                                                                acCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTg                                                                                              <  2:222009/77‑1 (MQ=255)

                                                                                                                                acCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTg                                                                                              >  1:220638/1‑77 (MQ=255)

                                                                                                                                acCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCgg                     >  2:1680822/1‑150 (MQ=255)

                                                                                                                                acCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCgg                     >  1:1270146/1‑150 (MQ=255)

                                                                                                                                acCCAGGCCTGGGTCGCGACGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTg                                                                                              <  2:220638/77‑1 (MQ=255)

                                                                                                                                  ccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGTCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCa                   >  1:241914/1‑150 (MQ=255)

                                                                                                                                  ccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCa                   >  1:1529817/1‑150 (MQ=255)

                                                                                                                                     cctcTGGGTCGAGTCGCGACCCAGGGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCa                >  1:915982/4‑150 (MQ=255)

                                                                                                                                      ctcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCaa               >  1:1560051/3‑150 (MQ=255)

                                                                                                                                       tcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAAt              >  1:598361/2‑150 (MQ=255)

                                                                                                                                       tcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAAt              >  1:658929/2‑150 (MQ=255)

                                                                                                                                        cTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCgg                     >  1:2160591/1‑142 (MQ=255)

                                                                                                                                        cTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCgg                     <  2:2160591/142‑1 (MQ=255)

                                                                                                                                        cTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCg                                                                               >  2:1699479/1‑84 (MQ=255)

                                                                                                                                        cTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCg                                                                               <  1:1699479/84‑1 (MQ=255)

                                                                                                                                        cTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATa             <  2:598361/150‑1 (MQ=255)

                                                                                                                                         tGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATAc            <  2:1529817/150‑1 (MQ=255)

                                                                                                                                          gggTCGCCACTCGACCCGGAGGCTCGTTGTGATGGTCGCCAATGTGCACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACt           <  1:631829/150‑1 (MQ=255)

                                                                                                                                           ggTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTg          <  1:1680822/150‑1 (MQ=255)

                                                                                                                                             tCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGta        >  1:3780/1‑150 (MQ=255)

                                                                                                                                               gcgACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGtata      >  1:2226909/1‑150 (MQ=255)

                                                                                                                                                 gACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGTATaaa    <  1:215892/150‑1 (MQ=255)

                                                                                                                                                  aCTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGTCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTg          >  1:1972471/1‑143 (MQ=255)

                                                                                                                                                  aCTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTg          <  2:1972471/143‑1 (MQ=255)

                                                                                                                                                   cTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGTATAAAGa  >  1:2237186/1‑150 (MQ=255)

                                                                                                                                                    tCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAAt              <  2:2280113/137‑1 (MQ=255)

ACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGTATAAAGA  >  NC_002947/4375342‑4375638

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