BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,502,790	T→G	6.2%	M141R (ATG→AGG)	PP_1310 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,502,790	0	T	G	6.2%	582.3 / 21.0	217	M141R (ATG→AGG)	PP_1310	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (101/97); new base G (11/4); total (116/101)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 223 total aligned reads displayed.

TGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCG  >  NC_002947/1502641‑1502939

tGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGACCGCCGGGGCGATCAt                                                                                                                                                       <  2:406377/150‑1 (MQ=255)

 gcagcaTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATg                                                                                                                                                      <  2:3344292/150‑1 (MQ=255)

 gcagcaTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATg                                                                                                                                                      <  1:2732895/150‑1 (MQ=255)

  cagcaTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATgc                                                                                                                                                     <  2:1333221/150‑1 (MQ=255)

    gcaTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCCTgcgc                                                                                                                                                   >  2:4736270/1‑150 (MQ=255)

     caTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTGATgggcg                                                                                                                                                  >  2:1661948/1‑149 (MQ=255)

     caTGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGcc                                                                                                                                                  <  1:4440431/150‑1 (MQ=255)

          cTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGAcc                                                                                                                                             <  1:920824/150‑1 (MQ=255)

           tGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGGGGGGAGGCGCCGGTccg                                                                                                                                            >  1:1000548/1‑150 (MQ=255)

           tGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGAccg                                                                                                                                            >  2:2186147/1‑150 (MQ=255)

             cGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGAccgcc                                                                                                                                          <  1:1078329/150‑1 (MQ=255)

             cGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGAccgcc                                                                                                                                          <  2:1969124/150‑1 (MQ=255)

             cGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGAccgcc                                                                                                                                          <  1:1770976/150‑1 (MQ=255)

               tGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGa                                                                                                                                        <  2:5245608/150‑1 (MQ=255)

               tGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGGCGCATGAccgccgg                                                                                                                                        <  1:5420797/150‑2 (MQ=255)

               tGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGGCGCATGAccgccgg                                                                                                                                        <  1:4294096/150‑2 (MQ=255)

               tGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGGCGTGCGCCCGGCCGCCGa                                                                                                                                        >  1:247242/1‑150 (MQ=255)

                   aGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGGGGGCGTGCGCCTGAccgccgccatc                                                                                                                                    >  1:1434527/1‑145 (MQ=255)

                     gCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTccc                                                                                                                                  <  1:5462835/150‑1 (MQ=255)

                     gCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAAATccc                                                                                                                                  >  2:5194365/1‑150 (MQ=255)

                     gCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGGCATGCGCCTGACCGCCGAATTccc                                                                                                                                  >  2:3413618/1‑150 (MQ=255)

                       aCCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCTGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGa                                                                                                                                <  2:4926083/150‑1 (MQ=255)

                       aCCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGa                                                                                                                                <  1:4361029/150‑1 (MQ=255)

                       aCCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGa                                                                                                                                <  2:5183316/150‑1 (MQ=255)

                       aCCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGGCGTGCGCCGGACCGCCCCATTCCCGa                                                                                                                                >  1:1280110/1‑150 (MQ=255)

                       aCCCCGGGCTACCGCTGGTGCTGATGGGCTGTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCTGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGa                                                                                                                                <  2:3841436/150‑1 (MQ=255)

                         cccGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGGCCTGGGCCGGACCGCCGCAGTCCCGATg                                                                                                                              >  1:3836370/1‑150 (MQ=255)

                           cGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCATGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGcc                                                                                                                            <  2:302549/150‑1 (MQ=255)

                                 aCCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAGGCGCCAGACCGCCGACATCCCGCTGCcgccgcc                                                                                                                      >  1:1354757/1‑149 (MQ=255)

                                 aCCGCAGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCCTGCGCCTGACCGCCCCATTCCCGATCCCGCcgcg                                                                                                                      >  2:5324837/1‑150 (MQ=255)

                                     cTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCGGACCGCCGACTTCCCGATGCCCCAgcgcgcc                                                                                                                  >  1:4099608/1‑147 (MQ=255)

                                     cTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGGCCTGCGCCGGGCCGCCGAATTCACCAAGCCGCAgcccggc                                                                                                                  >  1:2010562/1‑145 (MQ=255)

                                      tGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGcc                                                                                                                 >  1:4040443/1‑150 (MQ=255)

                                            tGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCAAGCGCCTGACCGCCGAAATCCCCATGCCGCAGCCCTGCCCGCTCa                                                                                                           >  2:4662961/1‑150 (MQ=255)

                                              aTGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCAcc                                                                                                         <  1:4139039/150‑1 (MQ=255)

                                                  gCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtgg                                                                                                     <  1:3891469/150‑1 (MQ=255)

                                                  gCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtgg                                                                                                     >  2:5085155/1‑150 (MQ=255)

                                                   cTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGACCGCCGGGGCGATCATgaaca                                                                                                                                                  >  1:5070136/1‑100 (MQ=38)

                                                     tttCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGGCAAGCCCCAGACCGCCGAAGACCCCCTGCCCCAGCGCAGCCCCCTCACCGACGTGc                                                                                                  >  1:784278/1‑150 (MQ=255)

                                                         gttcgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCgtgg                                                                                                     <  2:1796100/143‑1 (MQ=255)

                                                         gttcgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGGCATGCGCCTGACCGCCGCATTCCCGATGCCGCAGCGCAGCCCGCGCACCGGGGTGCAGcc                                                                                              >  1:3692975/1‑150 (MQ=255)

                                                          ttcgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCCTGCGCCAGACCCCCCACATCCCGATGCCGCCCCGCGGCCCCCGCACCGGGGGGccgccg                                                                                             >  2:174472/1‑150 (MQ=255)

                                                           tcgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCgg                                                                                            <  1:3231160/150‑1 (MQ=255)

                                                            cgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGa                                                                                           <  2:3759343/150‑1 (MQ=255)

                                                            cgGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGGCTTGCGCCTGACCGCCGACATCCCGATGCCGCAGCGCAGCGCGCGCACCCTGGTGCAGCCggc                                                                                           >  1:66201/1‑149 (MQ=255)

                                                             ggtggtTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGAc                                                                                          <  2:3588358/150‑1 (MQ=255)

                                                              gtggtTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGGCATGCGCCTGACCGCCGAATTACCCAAGCCGCACCGCAGCCCGCACACCGACGTGCAGCCGGACg                                                                                         >  1:4059115/1‑150 (MQ=255)

                                                                  tttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAAATCCCCATGCCGCAGCGCTGCCCGCTCACCGAGGTGCAGCCGGACGCCGa                                                                                     >  1:1293768/1‑150 (MQ=255)

                                                                   ttCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGGGGTCATGCGCCTGACCGCCGAAATCCCCATGCCGCACCGCTGCCCGCTCACCGAGGTGCAGCCGGACGCCgcg                                                                                    >  1:3476897/1‑148 (MQ=255)

                                                                   gtCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCGCACCGTGGTGCAGCCGGACGCCGAt                                                                                    <  2:4815817/149‑1 (MQ=255)

                                                                      gTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGACTTCCCGAAGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGCTGaa                                                                                 >  1:290238/1‑150 (MQ=255)

                                                                      gTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGaa                                                                                 <  1:2180469/150‑1 (MQ=255)

                                                                      gTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGaa                                                                                 <  1:379248/150‑1 (MQ=255)

                                                                          ccaccaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtcg                                                                             >  2:1463340/1‑150 (MQ=255)

                                                                           caccaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtcgt                                                                            <  1:1841479/150‑1 (MQ=255)

                                                                            accaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAgtcgtc                                                                           >  2:1287398/1‑150 (MQ=255)

                                                                            accaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAAATCCCCATGCCGCAGCGCTGCCCGCGCACCCGGGGGCAGCCCGACGCCGAGGGAgtcgtc                                                                           >  1:4901371/1‑150 (MQ=255)

                                                                             ccaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGtcgtcg                                                                          <  1:5292511/150‑1 (MQ=255)

                                                                             ccaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGtcgtcg                                                                          >  2:1983654/1‑150 (MQ=255)

                                                                             ccaGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCGGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCCACGTGGTGCAGCCGGCAGCCGATGGAGtagtcg                                                                          >  2:632177/1‑150 (MQ=255)

                                                                                       cGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCt                                                                <  1:36046/150‑1 (MQ=255)

                                                                                        gggCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGCTCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCACTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTg                                                               <  1:4424244/150‑1 (MQ=255)

                                                                                        gggCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGGCGTGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGcgg                                                               >  1:3059890/1‑148 (MQ=255)

                                                                                           cggcTTGAAACTGCCGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCTTGCGCGTGACCTCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCTGCCGGTCGCCGATGAAGTCGTCGCGCCGCAGCTGGtt                                                            <  2:4430166/147‑1 (MQ=255)

                                                                                             ccTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTa                                                          <  2:4984793/150‑1 (MQ=255)

                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTAc                                                         >  1:4672874/1‑150 (MQ=255)

                                                                                              cTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTAc                                                         <  1:3376928/150‑1 (MQ=255)

                                                                                               ttGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACg                                                        >  2:852328/1‑150 (MQ=255)

                                                                                               ttGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAAATACCGATGCCGCAGCGCTGCCCCCTCACCGAGGTGCAGCCGGACGCCGAAGAACTCGTCGCGCCGCAGCTGGTTTACg                                                        >  1:5620576/1‑150 (MQ=255)

                                                                                               ttGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCCCCCAAATCCCGATGCCGCAGCGCGGCCCGCGCACCCTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCGGGTTTaag                                                        >  1:2202506/1‑148 (MQ=255)

                                                                                               ttGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGGCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCCCTCACCCACGTGCAGCCGGACGCCCCTGAAGTCGTCGCGCCGCAGCTGGTTTAca                                                        >  1:4984082/1‑149 (MQ=255)

                                                                                                 gAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGaa                                                      <  1:1865817/150‑1 (MQ=255)

                                                                                                 gAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGAAGCCGCAGCGCTGCCCGCGCACCCGGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCCGCTGGTTTACGaa                                                      >  1:293681/1‑150 (MQ=255)

                                                                                                           gCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCCCCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGAcc                                            >  2:5382671/1‑149 (MQ=255)

                                                                                                           gCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACa                                            <  1:5231904/150‑1 (MQ=255)

                                                                                                            cGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGACTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAg                                           <  2:3745522/150‑1 (MQ=255)

                                                                                                             gTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGc                                          <  1:2342683/150‑1 (MQ=255)

                                                                                                             gTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGc                                          <  1:1809261/150‑1 (MQ=255)

                                                                                                             gTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCGGCCCGCGCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACCGc                                          >  2:153343/1‑150 (MQ=255)

                                                                                                                   ttGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGGGGTGCAGCCGGACGCCCATGAACTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTcg                                    >  1:1665637/1‑150 (MQ=255)

                                                                                                                      cagcaCCTGTTCATGATCGCCCCGGCGGTCATGCGCCAGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCCGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGt                                 >  2:1999301/1‑150 (MQ=255)

                                                                                                                       agcaCCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTc                                <  2:5163498/150‑1 (MQ=255)

                                                                                                                         caCCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTccc                              <  1:4062987/150‑1 (MQ=255)

                                                                                                                         caCCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTccc                              <  2:4777172/150‑1 (MQ=255)

                                                                                                                            cTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCAt                           <  2:2741708/150‑1 (MQ=255)

                                                                                                                                    gATCGCCCCGGCGGGCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCGCACCGTGGGGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCGGTCGCGTCCCCAAGAgctgct                   >  2:1870540/1‑150 (MQ=255)

                                                                                                                                     aTCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGctgctg                  <  1:4575245/150‑1 (MQ=255)

                                                                                                                                       cGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGaa                <  1:4562010/150‑1 (MQ=255)

                                                                                                                                        gCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGaaa               <  2:2474762/150‑1 (MQ=255)

                                                                                                                                          ttcggcggTCAGGCGCATGACCGCCGGATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGt             >  2:2831529/3‑150 (MQ=255)

                                                                                                                                          ttcggcggTCAGGCGCATGACCGCCGGATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTACTGAAAGt             >  2:319575/3‑150 (MQ=255)

                                                                                                                                          cccGGCGGTCATGCGCCTGACCGCCGAATTCCCGATTCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGAGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGt             >  2:5518035/1‑150 (MQ=255)

                                                                                                                                            cggcggTCATGCGCCTGCCCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTgg           <  1:2809210/150‑1 (MQ=255)

                                                                                                                                            cggcggTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTgg           >  2:3854535/1‑150 (MQ=255)

                                                                                                                                              gcggTCAGGCGCATGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCCCCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGGTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCa         >  2:4294143/1‑150 (MQ=255)

                                                                                                                                              gcggTCAGGCGCATGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCa         >  2:557844/1‑150 (MQ=255)

                                                                                                                                               cggTCAGGCGCATGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTa                                                          >  2:2156727/1‑100 (MQ=255)

                                                                                                                                                 gTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGaa      <  2:3476707/150‑1 (MQ=255)

                                                                                                                                                    aTGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGc   >  2:4360968/1‑150 (MQ=255)

                                                                                                                                                     tGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCCGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCg  <  1:3712242/150‑1 (MQ=255)

TGCAGCATGGCTGCGTGCCAGGCACCCCGGGCTACCGCTGGTGCTGATGGGCTTTTCGTTCGGTGGTTTCGTTGCCACCAGCCTGGCCGGGCGCCTTGAAACTGCAGGCGTCGAGTTGCAGCACCTGTTCATGATCGCCCCGGCGGTCATGCGCCTGACCGCCGAATTCCCGATGCCGCAGCGCTGCCCGCTCACCGTGGTGCAGCCGGACGCCGATGAAGTCGTCGCGCCGCAGCTGGTTTACGAATGGTCCGACAGCCTGTCGCGTCCCCATGAGCTGCTGAAAGTGGCAGAATGCG  >  NC_002947/1502641‑1502939

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