BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,699,409	G→A	6.2%	S15S (AGC→AGT) ‡	PP_5595 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,699,409	0	G	A	6.2%	616.7 / 18.5	194	S15S (AGC→AGT) ‡	PP_5595	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (87/94); new base A (7/5); total (95/99)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

GCATGTCCTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTCACCACTGCC  >  NC_002947/3699261‑3699557

   tGTCCTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTgg                                                                                                                                                  <  1:2139901/150‑1 (MQ=255)

    gTCCTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGc                                                                                                                                                 >  1:2933566/1‑150 (MQ=255)

       cTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAg                                                                                                                                              >  2:382070/1‑150 (MQ=255)

       cTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAg                                                                                                                                              <  1:2091705/150‑1 (MQ=255)

        ttGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAgg                                                                                                                                             >  1:3116904/1‑150 (MQ=255)

        ttGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAgg                                                                                                                                             <  2:3798882/150‑1 (MQ=255)

        ttGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAgg                                                                                                                                             >  2:4273472/1‑150 (MQ=255)

         tGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCAGCGGCCAggt                                                                                                                                            <  1:4174164/150‑2 (MQ=255)

           cTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAG‑CCTGGCCAGCGGCCAggtctg                                                                                                                                         <  2:606723/150‑5 (MQ=255)

            tttCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTgg                                                                                                                                         >  1:459099/1‑150 (MQ=255)

             ttCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGc                                                                                                                                        >  1:1060445/1‑150 (MQ=255)

             ttCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCAGCGGCCAggtctgg                                                                                                                                        <  2:4505299/150‑6 (MQ=255)

               cGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGGCCTGGCCAGCGGCCAggtctgggc                                                                                                                                      <  1:2808724/150‑8 (MQ=255)

                  tgctgcATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACca                                                                                                                                   <  1:3340486/150‑1 (MQ=255)

                    ctgcATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGTCGCTGGCCAGGCTGGCAACcatc                                                                                                                                 <  2:2703050/150‑1 (MQ=255)

                    ctgcATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACcatc                                                                                                                                 >  2:129611/1‑150 (MQ=255)

                     tgcATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACcatca                                                                                                                                <  2:2207476/150‑1 (MQ=255)

                       cATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTAGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCCCTGGCCAGGCTGGCAACcatcatc                                                                                                                              >  1:2711662/1‑150 (MQ=255)

                         tGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACcatcatcat                                                                                                                            <  1:1712336/150‑1 (MQ=255)

                            ccAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCaa                                                                                                                         <  2:823480/150‑1 (MQ=255)

                               gTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGc                                                                                                                      >  2:3483858/1‑150 (MQ=255)

                                ttGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCa                                                                                                                     >  2:2139200/1‑150 (MQ=255)

                                ttGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCa                                                                                                                     <  2:1623275/150‑1 (MQ=255)

                                    tCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGaa                                                                                                                 <  2:4004451/150‑1 (MQ=255)

                                    tCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGACCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGaa                                                                                                                 >  1:1310213/1‑150 (MQ=255)

                                     cGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAac                                                                                                                <  1:4415589/150‑1 (MQ=255)

                                           ggcggGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCAt                                                                                                          >  1:4169978/1‑150 (MQ=255)

                                           ggcggGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCAt                                                                                                          >  1:4170017/1‑150 (MQ=255)

                                              gggTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCAtttt                                                                                                       >  1:2990687/1‑150 (MQ=255)

                                               ggTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTc                                                                                                      <  2:2323949/150‑1 (MQ=255)

                                                 tATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCa                                                                                                    <  1:1218949/150‑1 (MQ=255)

                                                  aTGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCAc                                                                                                   <  2:4307306/150‑1 (MQ=255)

                                                   tGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCAcc                                                                                                  >  1:3957649/1‑150 (MQ=255)

                                                   tGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCAcc                                                                                                  >  2:1797463/1‑150 (MQ=255)

                                                    gTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCt                                                                                                 >  2:1919709/1‑150 (MQ=255)

                                                       gctgtagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCa                                                                                              >  2:3218640/1‑150 (MQ=255)

                                                       gctgtagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCa                                                                                              >  1:1282174/1‑150 (MQ=255)

                                                        ctgtagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAg                                                                                             <  2:866085/150‑1 (MQ=255)

                                                        ctgtagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAg                                                                                             >  2:2828215/1‑150 (MQ=255)

                                                         tgtagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAgg                                                                                            >  2:927788/1‑150 (MQ=255)

                                                          gtagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAggg                                                                                           >  2:3593945/1‑150 (MQ=255)

                                                           tagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGc                                                                                          >  2:301237/1‑150 (MQ=255)

                                                                gtaACAGGTGTGCCGCTCCACCAGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGACCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAg                                                                                     >  2:4179179/1‑150 (MQ=255)

                                                                 taACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAgg                                                                                    <  1:1169633/150‑1 (MQ=255)

                                                                  aaCAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGc                                                                                   <  1:3329719/150‑1 (MQ=255)

                                                                      ggTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTc                                                                               <  1:4088567/150‑1 (MQ=255)

                                                                      ggTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTc                                                                               <  2:1830310/150‑1 (MQ=255)

                                                                      ggTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTc                                                                               <  2:1385925/150‑1 (MQ=255)

                                                                      ggTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTc                                                                               <  2:3924773/150‑1 (MQ=255)

                                                                               cTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCt                                                                      <  1:3316531/150‑1 (MQ=255)

                                                                                 ccACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGa                                                                    <  1:52278/150‑1 (MQ=255)

                                                                                 ccACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGa                                                                    <  2:3685742/150‑1 (MQ=255)

                                                                                    ccGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCa                                                                 <  1:3369702/150‑1 (MQ=255)

                                                                                    ccGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCa                                                                 <  2:1515932/150‑1 (MQ=255)

                                                                                     cGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAg                                                                <  2:1302673/150‑1 (MQ=255)

                                                                                      ggCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAgg                                                               <  2:372420/150‑1 (MQ=255)

                                                                                         gTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCtt                                                            <  1:1995406/150‑1 (MQ=255)

                                                                                          tCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCATGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTg                                                           >  2:295551/1‑150 (MQ=255)

                                                                                            aCGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTc                                                         <  1:2239564/150‑1 (MQ=255)

                                                                                            aCGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTc                                                         >  2:3490640/1‑150 (MQ=255)

                                                                                              ggCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCCGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTcgc                                                       >  2:4419302/1‑150 (MQ=255)

                                                                                              ggCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACcatcat                                                                                                                               >  1:4294422/1‑78 (MQ=255)

                                                                                              ggCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACcatcat                                                                                                                               <  2:4294712/78‑1 (MQ=255)

                                                                                                  tCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCg                                                   >  1:1873963/1‑150 (MQ=255)

                                                                                                     cTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGgcg                                                >  2:1209089/1‑150 (MQ=255)

                                                                                                      tGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGgcgc                                               >  1:2858982/1‑150 (MQ=255)

                                                                                                        gaggagCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCAt                                             >  2:1360469/1‑150 (MQ=255)

                                                                                                         aggagCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATg                                            <  1:1002787/150‑1 (MQ=255)

                                                                                                          ggagCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATgg                                           <  1:947333/150‑1 (MQ=255)

                                                                                                              cGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACATGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCttt                                       <  2:4080750/150‑1 (MQ=255)

                                                                                                              cGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCttt                                       <  2:2520561/150‑1 (MQ=255)

                                                                                                               gAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTg                                      <  2:1060445/150‑1 (MQ=255)

                                                                                                               gAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTg                                      <  2:1479172/150‑1 (MQ=255)

                                                                                                                   tgtcacaTTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGgcgc                                  <  1:3780491/147‑1 (MQ=255)

                                                                                                                   cAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGgcgc                                  <  1:635448/150‑1 (MQ=255)

                                                                                                                     gCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTGTGGCGCGt                                >  2:2555817/1‑150 (MQ=255)

                                                                                                                         aTTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGttttt                            >  2:1250872/1‑150 (MQ=255)

                                                                                                                          ttGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTc                           >  1:3911185/1‑150 (MQ=255)

                                                                                                                            gCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCa                         >  2:3605866/1‑150 (MQ=255)

                                                                                                                            gCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCa                         <  2:353862/150‑1 (MQ=255)

                                                                                                                            gCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCa                         <  1:485083/150‑1 (MQ=255)

                                                                                                                              aaGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAgg                       >  2:2171066/1‑150 (MQ=255)

                                                                                                                              aaGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAgg                       >  1:2632430/1‑150 (MQ=255)

                                                                                                                              aaGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAgg                       >  2:2008402/1‑150 (MQ=255)

                                                                                                                              aaGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTGTGGCGCGTTTTTCCAgg                       >  1:680860/1‑150 (MQ=255)

                                                                                                                                       cccAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCaggtag              <  2:1864433/150‑1 (MQ=255)

                                                                                                                                       cccAGACCTGGCCGCTGGCCAGGCGGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCCGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCaggtag              <  2:4117154/150‑1 (MQ=255)

                                                                                                                                         ccagccTGGCCAGCGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCaggtaggt            >  1:1199129/5‑150 (MQ=255)

                                                                                                                                         cAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCaggtaggt            <  2:224470/150‑1 (MQ=255)

                                                                                                                                          cagccTGGCCAGCGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTc           >  1:3648898/4‑150 (MQ=255)

                                                                                                                                          cagccTGGCCAGCGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTc           >  1:1569178/4‑150 (MQ=255)

                                                                                                                                          cagccTGGCCAGCGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTc           >  1:1371552/4‑150 (MQ=255)

                                                                                                                                            gccTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTcac         >  2:3672028/2‑150 (MQ=255)

                                                                                                                                            aCCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTcac         <  2:1282174/150‑1 (MQ=255)

                                                                                                                                              cTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTcacca       >  2:3972713/1‑150 (MQ=255)

                                                                                                                                               tGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTcaccac      >  2:3927860/1‑150 (MQ=255)

                                                                                                                                               tGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTcaccac      >  2:647678/1‑150 (MQ=255)

                                                                                                                                                  ccGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTCACCACTGc   <  1:4331800/150‑1 (MQ=255)

                                                                                                                                                  ccGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTCACCACTGc   >  2:1125451/1‑150 (MQ=255)

                                                                                                                                                  ccGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTCACCACTGc   <  1:2323368/150‑1 (MQ=255)

GCATGTCCTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTCACCACTGCC  >  NC_002947/3699261‑3699557

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