BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,983,700	C→G	5.7%	F45L (TTC→TTG)	cmk →	cytidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,983,700	0	C	G	5.7%	1205.9 / 44.1	362	F45L (TTC→TTG)	cmk	cytidylate kinase
Reads supporting (aligned to +/- strand): ref base C (169/172); new base G (4/17); total (173/189)
Fisher's exact test for biased strand distribution p-value = 6.70e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 368 total aligned reads displayed.

AGAAGAGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCCGTAC  >  NC_002947/1983551‑1983849

  aagaGGGCAAGTCGTGAATTCAATGACACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCaa                                                                                                                                                     <  1:7158360/150‑1 (MQ=255)

   agagGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGTAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAAt                                                                                                                                                    <  1:5662664/150‑1 (MQ=255)

   agagGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAAt                                                                                                                                                    <  1:1252395/150‑1 (MQ=255)

       ggCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCa                                                                                                                                                <  1:4713505/150‑1 (MQ=255)

       ggCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCa                                                                                                                                                <  2:1625485/150‑1 (MQ=255)

           aGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaa                                                                                                                                            <  1:8604013/150‑1 (MQ=255)

            gTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaag                                                                                                                                           <  2:210527/150‑2 (MQ=255)

               gTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCAc                                                                                                                                        <  2:1837822/150‑1 (MQ=255)

                tGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACg                                                                                                                                       <  1:4256959/150‑1 (MQ=255)

                tGAATTCAAGGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACg                                                                                                                                       <  1:8689182/150‑1 (MQ=255)

                  aaTTCTATGCCACGGGTTATCACCCTCGACGGACCATGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGAGTGCCCGCGAGCTGGGCTGGAAGCTGCGGGACTCCTGGGCGCTGTATCGCTTGCTGGCATTCTATGCCAGAAACCACGGc                                                                                                                                     <  2:5531391/150‑1 (MQ=255)

                  aaTTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGTTTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGc                                                                                                                                     <  1:2574468/150‑1 (MQ=255)

                       aaTGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGa                                                                                                                                >  1:7588138/1‑150 (MQ=255)

                       aaTGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGa                                                                                                                                >  1:138354/1‑150 (MQ=255)

                       aaTGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGa                                                                                                                                >  1:712277/1‑150 (MQ=255)

                            cACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCTGCAACCACGGCGTCGACCTGa                                                                                                                           >  2:1523000/1‑150 (MQ=255)

                             aCCGGTTACCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTTATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTTTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGAc                                                                                                                          <  1:2512029/150‑1 (MQ=255)

                              ccGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGAcc                                                                                                                         <  2:84758/150‑1 (MQ=255)

                                     tCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAag                                                                                                                  >  1:4922268/1‑150 (MQ=255)

                                     tCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTAACGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAag                                                                                                                  <  2:7776509/150‑1 (MQ=255)

                                       aCCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAagag                                                                                                                >  1:878759/1‑150 (MQ=255)

                                        ccATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTTTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgc                                                                                                               <  1:6822460/150‑1 (MQ=255)

                                        ccATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgc                                                                                                               <  2:234545/150‑1 (MQ=255)

                                        ccATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgc                                                                                                               <  1:7296016/150‑1 (MQ=255)

                                         cATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgct                                                                                                              >  2:3206457/1‑150 (MQ=255)

                                          aTCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGTATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgctg                                                                                                             <  1:4437220/150‑1 (MQ=255)

                                           tCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgctgc                                                                                                            <  2:4825651/150‑1 (MQ=255)

                                              aCGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGTTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGa                                                                                                         >  2:8790192/1‑150 (MQ=255)

                                                  aCCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCaa                                                                                                     >  1:8149514/1‑150 (MQ=255)

                                                    cAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGc                                                                                                   >  2:2728274/1‑150 (MQ=255)

                                                     aaGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCt                                                                                                  >  2:1442672/1‑150 (MQ=255)

                                                       gCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTg                                                                                                >  2:2443409/1‑150 (MQ=255)

                                                        cGGCTCGGGCAAGGGCACGGTCGCCGGGCTGTTGGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgc                                                                                               <  1:457466/150‑1 (MQ=255)

                                                         ggCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgcc                                                                                              <  2:7009266/150‑1 (MQ=255)

                                                           cTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgccgc                                                                                            <  2:2505217/150‑1 (MQ=255)

                                                              gggCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGTGCTGCTGACCAAGCTTGCCGCCCa                                                                                         >  1:3808706/1‑150 (MQ=255)

                                                               ggCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCAt                                                                                        >  2:6430813/1‑150 (MQ=255)

                                                                gCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATc                                                                                       >  2:4331380/1‑150 (MQ=255)

                                                                gCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATc                                                                                       >  1:7914802/1‑150 (MQ=255)

                                                                    gggCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGa                                                                                   <  1:3493211/150‑1 (MQ=255)

                                                                    gggCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGa                                                                                   >  2:5484564/1‑150 (MQ=255)

                                                                        aCGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtgtg                                                                               >  2:861121/1‑150 (MQ=255)

                                                                        aCGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtgtg                                                                               <  1:6544445/150‑1 (MQ=255)

                                                                        aCGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtgtg                                                                               <  2:4549775/150‑1 (MQ=255)

                                                                          ggTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCa                                                                             >  2:5032425/1‑150 (MQ=255)

                                                                           gTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAg                                                                            >  1:5364552/1‑150 (MQ=255)

                                                                             cGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGtt                                                                          <  2:1183370/150‑1 (MQ=255)

                                                                                 gggCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATc                                                                      <  2:8125332/150‑1 (MQ=255)

                                                                                 gggCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATc                                                                      <  1:8462450/150‑1 (MQ=255)

                                                                                 gggCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATc                                                                      <  2:3942695/150‑1 (MQ=255)

                                                                                    cTGATGGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCGGTATCGCTTGCGGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCt                                                                   <  1:3205957/150‑1 (MQ=255)

                                                                                      gATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGc                                                                 >  1:520130/1‑150 (MQ=255)

                                                                                       aTTGCCCGCGTGCTGGGGTGGATGCTGGTGGACTCCGGCGCGGTGTATCGGTTGCTGGCTTTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGcc                                                                <  1:4066807/150‑1 (MQ=255)

                                                                                         tGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTTTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGa                                                              <  1:1199316/150‑1 (MQ=255)

                                                                                            ccGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGcc                                                           >  1:2877310/1‑150 (MQ=255)

                                                                                            ccGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGcc                                                           <  2:3478539/150‑1 (MQ=255)

                                                                                              gcgAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCg                                                         >  1:3070633/1‑150 (MQ=255)

                                                                                                   ctgggctggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAg                                                    <  1:6049916/150‑1 (MQ=255)

                                                                                                      ggctggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCATGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTa                                                 >  1:3602549/1‑150 (MQ=255)

                                                                                                       gctggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTac                                                <  2:5952637/150‑1 (MQ=255)

                                                                                                         tggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTACGCTacaa                                              >  2:7105180/1‑150 (MQ=255)

                                                                                                         tggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTacaa                                              >  2:5257295/1‑150 (MQ=255)

                                                                                                          ggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGATAAGCTacaac                                             <  2:3708801/150‑1 (MQ=255)

                                                                                                           gAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTacaaca                                            <  2:3271842/150‑1 (MQ=255)

                                                                                                           gAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTacaaca                                            >  1:1786115/1‑150 (MQ=255)

                                                                                                           gAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTacaaca                                            <  2:7544986/150‑1 (MQ=255)

                                                                                                           gAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGAATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTacaaca                                            <  1:4789592/150‑1 (MQ=255)

                                                                                                             aGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGa                                          <  1:2426063/150‑1 (MQ=255)

                                                                                                             aGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGa                                          <  2:5368254/150‑1 (MQ=255)

                                                                                                                   tgGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATcc                                    >  1:1120157/1‑150 (MQ=255)

                                                                                                                    ggACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCAAGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCAACTGGATGTGCAGTACATCGCAGCCGAGCCCGGTAAGCTACAACAGATCATCCt                                   >  1:3660302/1‑150 (MQ=255)

                                                                                                                     gACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCtt                                  >  1:6297265/1‑150 (MQ=255)

                                                                                                                       cTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGa                                <  2:2977552/150‑1 (MQ=255)

                                                                                                                       cTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGa                                <  2:6509056/150‑1 (MQ=255)

                                                                                                                        tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAg                               >  2:786945/1‑150 (MQ=255)

                                                                                                                        tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAg                               >  2:4436358/1‑150 (MQ=255)

                                                                                                                        tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAg                               >  2:7473169/1‑150 (MQ=255)

                                                                                                                        tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAg                               >  1:1940983/1‑150 (MQ=255)

                                                                                                                        tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAg                               <  1:5572105/150‑1 (MQ=255)

                                                                                                                        tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAg                               >  2:287546/1‑150 (MQ=255)

                                                                                                                         ccGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAgg                              <  1:6155798/150‑1 (MQ=255)

                                                                                                                         ccGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAgg                              <  2:8244961/150‑1 (MQ=255)

                                                                                                                         ccGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAgg                              <  1:7011993/150‑1 (MQ=255)

                                                                                                                         ccGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAgg                              <  2:6101963/150‑1 (MQ=255)

                                                                                                                          cGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATACGTGAggg                             >  2:4120262/1‑150 (MQ=255)

                                                                                                                           ggCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGt                            >  1:3865906/1‑150 (MQ=255)

                                                                                                                           ggCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGt                            >  2:7994329/1‑150 (MQ=255)

                                                                                                                              gcgcTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAg                         <  2:7440861/150‑1 (MQ=255)

                                                                                                                                  tGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACg                     >  2:2020525/1‑150 (MQ=255)

                                                                                                                                    tATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTc                   <  2:6831998/150‑1 (MQ=255)

                                                                                                                                     aTCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCa                  <  2:4922268/150‑1 (MQ=255)

                                                                                                                                     aTCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCa                  <  2:5281801/150‑1 (MQ=255)

                                                                                                                                       cGCTTGCTGGCATTCACTGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGGAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGc                <  1:1681829/150‑1 (MQ=255)

                                                                                                                                          ttGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAAc             <  1:947962/150‑1 (MQ=255)

                                                                                                                                          ttGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCGTGAGGGTGAGGACGTCAGCAAc             >  2:5094870/1‑150 (MQ=255)

                                                                                                                                             cTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGACGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTc          >  2:7151005/1‑150 (MQ=255)

                                                                                                                                              tGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCa         >  1:811949/1‑150 (MQ=255)

                                                                                                                                               ggCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCAt        <  1:6590690/150‑1 (MQ=255)

                                                                                                                                                  aTTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCCg     >  1:2088586/1‑150 (MQ=255)

                                                                                                                                                  aTTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCCg     >  2:8724499/1‑150 (MQ=255)

AGAAGAGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCCGTAC  >  NC_002947/1983551‑1983849

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