BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,986,918	G→A	5.9%	S60S (AGC→AGT)	ycaC‑II ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,986,918	0	G	A	5.9%	664.1 / 17.6	202	S60S (AGC→AGT)	ycaC‑II	putative hydrolase
Reads supporting (aligned to +/- strand): ref base G (86/103); new base A (6/6); total (92/110)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 206 total aligned reads displayed.

TCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGCGTCGTCTTT  >  NC_002947/5986770‑5987067

   gTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtgg                                                                                                                                                   <  1:1114095/150‑1 (MQ=255)

     tgcttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtg                                                                                                                                                 <  2:5102406/150‑1 (MQ=255)

        ttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGgtggtg                                                                                                                                              >  2:2083388/1‑150 (MQ=255)

         tgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                                                                             <  1:773346/150‑1 (MQ=255)

         tcgtgCCGGTGGCCTTCACTGCCTTGACGAAGTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTggtggtca                                                                                                                                             <  2:1225634/146‑3 (MQ=2)

                 gTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTc                                                                                                                                     <  2:1157140/150‑1 (MQ=255)

                     ccTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCa                                                                                                                                 <  1:37140/150‑1 (MQ=255)

                     ccTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCa                                                                                                                                 <  2:3650042/150‑1 (MQ=255)

                     ccTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCa                                                                                                                                 <  1:4781582/150‑1 (MQ=255)

                      cttcactGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTGGTGGTCAGAACAGTCGGCAg                                                                                                                                >  2:2000625/7‑150 (MQ=2)

                      cTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                                                                                                                                >  1:558135/1‑150 (MQ=255)

                        tGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGc                                                                                                                              <  1:4365405/150‑1 (MQ=255)

                         gATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGcc                                                                                                                             >  2:2441149/1‑150 (MQ=255)

                         gATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGAAGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGcc                                                                                                                             >  1:615106/1‑150 (MQ=255)

                            tGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTGGTGGTCAGAACAGTCGGCAgcttgaa                                                                                                                          <  1:4486982/150‑7 (MQ=2)

                            tGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCaaa                                                                                                                          <  2:3235303/150‑1 (MQ=255)

                               cTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGaa                                                                                                                       >  2:4322015/1‑150 (MQ=255)

                               cTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGaa                                                                                                                       >  1:2125481/1‑150 (MQ=255)

                               cTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGaa                                                                                                                       >  1:328101/1‑150 (MQ=255)

                               cTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGaa                                                                                                                       >  2:4322104/1‑150 (MQ=255)

                                ttGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAAc                                                                                                                      >  2:1648870/1‑150 (MQ=255)

                                ttGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAAc                                                                                                                      >  1:2614272/1‑150 (MQ=255)

                                  gtcgtaGTCTTCGTTGTCACAGGCGTTGATCTGGCCGGGACGGGCTATGTACGGCGCTTCCTGGAACATCTCTTTCAGCTGGGGAACCAGCTGGCCATTCGGGCCTTTCTCTAAGCTGGTGGTGAGGATTGTCGGCAGGCCAAAGAAAtt                                                                                                                    <  2:2728906/145‑1 (MQ=255)

                                  gACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACtt                                                                                                                    <  2:214587/150‑1 (MQ=255)

                                  gACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACtt                                                                                                                    <  1:1648870/150‑1 (MQ=255)

                                   aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                                                                                                   <  1:2877689/150‑1 (MQ=255)

                                   aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                                                                                                   <  2:2493216/150‑1 (MQ=255)

                                   aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                                                                                                   <  2:1394052/150‑1 (MQ=255)

                                   aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                                                                                                   <  1:621865/150‑1 (MQ=255)

                                     gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc                                                                                                                 >  2:2067071/1‑150 (MQ=255)

                                     gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc                                                                                                                 >  1:1145974/1‑150 (MQ=255)

                                      aaGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGcc                                                                                                                <  1:3135239/150‑1 (MQ=255)

                                        gTCTTCGTGGGCCCCGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTGTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAg                                                                                                              <  1:5043262/150‑1 (MQ=255)

                                         tCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAgg                                                                                                             >  1:3628679/1‑150 (MQ=255)

                                         tCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAgg                                                                                                             >  2:4505916/1‑150 (MQ=255)

                                              gTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGcc                                                                                                        >  1:205030/1‑150 (MQ=255)

                                                    ccAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAgac                                                                                                  <  1:2476973/150‑3 (MQ=255)

                                                     cAGGCGTTGAGCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGcc                                                                                                 >  2:1979738/1‑150 (MQ=255)

                                                        gCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGc                                                                                              <  2:4413128/150‑1 (MQ=255)

                                                        gCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGc                                                                                              >  1:1427800/1‑150 (MQ=255)

                                                          gTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCAc                                                                                            >  2:344018/1‑150 (MQ=255)

                                                              aTCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACgttg                                                                                        >  1:2730480/1‑150 (MQ=255)

                                                               tCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGTCCAGCACgttgt                                                                                       >  2:3691014/1‑150 (MQ=255)

                                                                 tGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTc                                                                                     >  1:4289839/1‑150 (MQ=255)

                                                                 tGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTc                                                                                     >  1:5107273/1‑150 (MQ=255)

                                                                 tGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTc                                                                                     >  1:1709553/1‑150 (MQ=255)

                                                                   gCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCtt                                                                                   <  2:3770914/150‑1 (MQ=255)

                                                                   gCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCtt                                                                                   <  2:2331535/150‑1 (MQ=255)

                                                                    ccGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTg                                                                                  <  2:2950690/150‑1 (MQ=255)

                                                                    ccGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATAGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTg                                                                                  <  2:3339171/150‑1 (MQ=255)

                                                                     cGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGa                                                                                 <  2:34929/150‑1 (MQ=255)

                                                                     cGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGa                                                                                 >  2:1762841/1‑150 (MQ=255)

                                                                     cGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGa                                                                                 <  2:4695719/150‑1 (MQ=255)

                                                                     cGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCCTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGa                                                                                 >  1:4539899/1‑150 (MQ=255)

                                                                            ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                                                                          <  1:1997437/150‑1 (MQ=255)

                                                                            ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                                                                          <  1:2266979/150‑1 (MQ=255)

                                                                            ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                                                                          <  2:1201335/150‑1 (MQ=255)

                                                                             gCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTg                                                                         <  2:848909/150‑1 (MQ=255)

                                                                             gCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTg                                                                         <  1:2401672/150‑1 (MQ=255)

                                                                             gCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTg                                                                         <  1:1919584/150‑1 (MQ=255)

                                                                                aTGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGc                                                                      <  2:1544656/150‑1 (MQ=255)

                                                                                aTGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGc                                                                      <  1:4613171/150‑1 (MQ=255)

                                                                                 tGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCg                                                                     >  2:1050842/1‑150 (MQ=255)

                                                                                   tACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGaa                                                                   <  1:1513163/150‑1 (MQ=255)

                                                                                       gcgcATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAAt                                                               >  2:77401/1‑150 (MQ=255)

                                                                                          cATCCGGGAACATCTCTTTCCGCTCGTGCACCAGCGGGCCATTCGGGCCTTGCTCGAATCTGGTGGTGAGGATGGTCGGCATGCCAAAGAACTTGGCCAGGTCTCCAAAGGCCAGCACGTTTTTCTTGAATTCGTTGGGCGAAAAATcgt                                                            <  2:2862337/150‑3 (MQ=25)

                                                                                           aTCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCTGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTg                                                           <  2:1133208/150‑1 (MQ=255)

                                                                                           aTCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTg                                                           <  1:2603258/150‑1 (MQ=255)

                                                                                           aTCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTg                                                           >  1:2981448/1‑150 (MQ=255)

                                                                                           aTCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTg                                                           <  1:316012/150‑1 (MQ=255)

                                                                                            tCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGc                                                          <  1:1589822/150‑1 (MQ=255)

                                                                                              cGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCAc                                                        >  2:1864985/1‑150 (MQ=255)

                                                                                                ggAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCa                                                      <  2:1217530/150‑1 (MQ=255)

                                                                                                      tctcTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCgaga                                                <  1:3980801/150‑1 (MQ=255)

                                                                                                       ctctTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGAt                                               >  2:2329441/1‑150 (MQ=255)

                                                                                                        tctTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATc                                              >  1:1279401/1‑150 (MQ=255)

                                                                                                          tttCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAg                                            >  1:498834/1‑150 (MQ=255)

                                                                                                           ttCGGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGa                                           <  2:48273/150‑1 (MQ=255)

                                                                                                           ttCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                                                                          <  1:3048665/119‑1 (MQ=255)

                                                                                                           ttCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGCTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                                                                          >  2:3048665/1‑119 (MQ=255)

                                                                                                            tCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGAc                                          >  1:1213247/1‑150 (MQ=255)

                                                                                                              aGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCg                                        <  1:2329441/150‑1 (MQ=255)

                                                                                                                cTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGt                                      >  1:2661742/1‑150 (MQ=255)

                                                                                                                   gggCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCtg                                   >  2:1796213/1‑150 (MQ=255)

                                                                                                                   gggCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCtg                                   <  2:4758533/150‑1 (MQ=255)

                                                                                                                     gCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCtggt                                 >  1:4261202/1‑150 (MQ=255)

                                                                                                                      cACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCtggtg                                <  1:4703405/150‑1 (MQ=255)

                                                                                                                           gCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGa                           <  2:1098511/150‑1 (MQ=255)

                                                                                                                            cGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGAc                          <  1:2083388/150‑1 (MQ=255)

                                                                                                                                        gCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGcaccac              <  2:5050045/150‑1 (MQ=255)

                                                                                                                                          cTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACgg            >  2:2390061/1‑150 (MQ=255)

                                                                                                                                          cTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACgg            >  1:1400096/1‑150 (MQ=255)

                                                                                                                                           ttGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGc           <  2:558135/150‑1 (MQ=255)

                                                                                                                                            tGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGcg          >  2:1031449/1‑150 (MQ=255)

                                                                                                                                              cTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGcgtc        >  1:3431971/1‑150 (MQ=255)

                                                                                                                                               tCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGcgtcg       >  1:944668/1‑150 (MQ=255)

                                                                                                                                               tCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGcgtcg       >  2:629369/1‑150 (MQ=255)

                                                                                                                                                 gAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGCgtcgtc     >  1:4085619/1‑150 (MQ=255)

                                                                                                                                                    gCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGCGTCGTCttt  <  2:1806135/150‑1 (MQ=255)

                                                                                                                                                    gCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGCGTCGTCttt  <  1:3667815/150‑1 (MQ=255)

TCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGCGTCGTCTTT  >  NC_002947/5986770‑5987067

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