BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,602,750	T→A	5.1%	K685M (AAG→ATG) ‡	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,602,750	0	T	A	5.1%	798.7 / 22.1	233	K685M (AAG→ATG) ‡	PP_4071	membrane protein
Reads supporting (aligned to +/- strand): ref base T (99/122); new base A (8/4); total (107/126)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 236 total aligned reads displayed.

AAGCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAGCTGCCCC  >  NC_002947/4602603‑4602899

aaGCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGTTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTc                                                                                                                                                     <  2:3408446/150‑1 (MQ=255)

aaGCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTc                                                                                                                                                     >  2:5289723/1‑150 (MQ=255)

aaGCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTc                                                                                                                                                     >  1:2404184/1‑150 (MQ=255)

 aGCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCa                                                                                                                                                    <  1:3992356/150‑1 (MQ=255)

 aGCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCa                                                                                                                                                    >  2:3607902/1‑150 (MQ=255)

  gCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAg                                                                                                                                                   <  2:3048920/150‑1 (MQ=255)

   ccTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAgg                                                                                                                                                  <  1:5045351/150‑1 (MQ=255)

   ccTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAgg                                                                                                                                                  <  2:168763/150‑1 (MQ=255)

    cTTTGGCAGTGAGCATTGGCGCCAGCAAATCACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGACCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGc                                                                                                                                                 <  2:76852/150‑1 (MQ=255)

     tttGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGcc                                                                                                                                                >  2:444379/1‑150 (MQ=255)

     tttGGCAGTGAGCATTGGCGCCAGAAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGATCATGGGCCTGCTGCTTCAGGcc                                                                                                                                                >  2:1365411/1‑150 (MQ=255)

         gCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCa                                                                                                                                            <  1:4477367/150‑1 (MQ=255)

         gCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCCGGGCCTGAAGCAGCAGGcccagg                                                                                                                                            <  1:2073288/150‑4 (MQ=255)

            gTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATac                                                                                                                                         >  2:4608697/1‑150 (MQ=255)

            gTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATac                                                                                                                                         >  1:3850864/1‑150 (MQ=255)

            gTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATac                                                                                                                                         <  1:3099829/150‑1 (MQ=255)

                     ggCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtgt                                                                                                                                <  1:3615164/150‑1 (MQ=255)

                      gcgcCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtgtg                                                                                                                               <  2:4913131/150‑1 (MQ=255)

                      gcgcCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtgtg                                                                                                                               >  1:5471372/1‑150 (MQ=255)

                      gcgcCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtgtg                                                                                                                               <  1:1756615/150‑1 (MQ=255)

                       cgcCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTgg                                                                                                                              <  2:5464018/150‑1 (MQ=255)

                             cAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACgg                                                                                                                        >  1:5309707/1‑150 (MQ=255)

                              aaaaCACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGc                                                                                                                       >  1:3165455/1‑150 (MQ=255)

                              aaaaCACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGc                                                                                                                       >  1:4843342/1‑150 (MQ=255)

                              aaaaCACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGc                                                                                                                       >  1:25035/1‑150 (MQ=255)

                                  cacCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCt                                                                                                                   <  1:2550560/150‑1 (MQ=255)

                                          tCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAAt                                                                                                           >  1:101228/1‑150 (MQ=255)

                                          gcGCCCTGGCTCGAGCCGGAAGACAGCTGGTTAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAAt                                                                                                           <  2:3344387/149‑1 (MQ=255)

                                            gCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGt                                                                                                         <  2:2728748/150‑1 (MQ=255)

                                             cccTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTc                                                                                                        >  1:2949258/1‑150 (MQ=255)

                                             cccTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTc                                                                                                        >  1:5258833/1‑150 (MQ=255)

                                               cTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTctg                                                                                                      <  1:1451706/150‑1 (MQ=255)

                                                tGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTctgc                                                                                                     <  1:3665369/150‑1 (MQ=255)

                                                   cTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCtgctgc                                                                                                  >  2:4825078/1‑150 (MQ=255)

                                                    tCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCACTGTCTGCTGCa                                                                                                 >  2:3303971/1‑150 (MQ=255)

                                                    tCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCa                                                                                                 >  1:4440666/1‑150 (MQ=255)

                                                      gAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACg                                                                                               >  2:602088/1‑150 (MQ=255)

                                                         ccGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTgg                                                                                            <  1:565243/150‑1 (MQ=255)

                                                          cGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGc                                                                                           >  2:5118711/1‑150 (MQ=255)

                                                           ggAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCt                                                                                          >  1:2832230/1‑150 (MQ=255)

                                                           ggAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCt                                                                                          >  1:909581/1‑150 (MQ=255)

                                                            gAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTg                                                                                         >  2:3215267/1‑150 (MQ=255)

                                                                 cAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACg                                                                                    <  2:3165455/150‑1 (MQ=255)

                                                                  aGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCtt                                                                                                                  >  2:4166788/1‑119 (MQ=255)

                                                                  aGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGc                                                                                   >  1:2846738/1‑150 (MQ=255)

                                                                   gCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCt                                                                                  <  1:3407345/150‑1 (MQ=255)

                                                                    cTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTg                                                                                 <  2:3901351/150‑1 (MQ=255)

                                                                     tGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTCTGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCCACGGCACCGTACGCAATGTCTGCTGCACGTGGCTGACACGCTgg                                                                                >  2:4855339/1‑150 (MQ=255)

                                                                     tGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGtt                                                                                                                                  <  2:4255595/100‑1 (MQ=255)

                                                                     tGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTgg                                                                                <  1:2268902/150‑1 (MQ=255)

                                                                     tGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTgg                                                                                <  1:1577558/150‑1 (MQ=255)

                                                                     tGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTgg                                                                                <  1:3299129/150‑1 (MQ=255)

                                                                     tGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTgg                                                                                >  1:1929850/1‑150 (MQ=255)

                                                                      ggTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTgcc                                                                               <  1:3030897/150‑3 (MQ=255)

                                                                             gACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTg                                                                        <  2:2832366/150‑1 (MQ=255)

                                                                              aCATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGc                                                                       >  1:4283302/1‑150 (MQ=255)

                                                                              aCATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGc                                                                       <  1:643059/150‑1 (MQ=255)

                                                                               cATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCg                                                                      <  1:2156061/150‑1 (MQ=255)

                                                                                aTCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACGGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCgg                                                                     <  2:765203/150‑1 (MQ=255)

                                                                                 tCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGc                                                                    >  2:4744098/1‑150 (MQ=255)

                                                                                 tCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGc                                                                    >  1:179794/1‑150 (MQ=255)

                                                                                   gAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAg                                                                  >  1:1135258/1‑150 (MQ=255)

                                                                                   gAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAg                                                                  >  2:672212/1‑150 (MQ=255)

                                                                                    aaCGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAgg                                                                 >  1:2452857/1‑150 (MQ=255)

                                                                                       gCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCg                                                              <  1:3507864/150‑1 (MQ=255)

                                                                                        cAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCgg                                                             <  2:1151529/150‑1 (MQ=255)

                                                                                         aGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGt                                                            <  1:4394001/150‑1 (MQ=255)

                                                                                           gCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTAt                                                          <  1:3200692/150‑1 (MQ=255)

                                                                                            ccccACCTGATGACGCAGGTCCAGTCCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATc                                                         <  2:260948/150‑1 (MQ=255)

                                                                                             cccACCTGATGACACAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATcc                                                        >  2:634949/1‑150 (MQ=255)

                                                                                                  cTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTa                                                   >  1:1341782/1‑150 (MQ=255)

                                                                                                    gatgaCGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAgg                                                 <  2:243881/150‑1 (MQ=255)

                                                                                                    gatgaCGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAgg                                                 >  1:1300681/1‑150 (MQ=255)

                                                                                                          gCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTCGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAgcg                                           <  1:225963/150‑1 (MQ=255)

                                                                                                                 cAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACtt                                    >  1:4610708/1‑150 (MQ=255)

                                                                                                                  aGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTc                                   <  1:4080896/150‑1 (MQ=255)

                                                                                                                   gACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCa                                  <  1:4188802/150‑1 (MQ=255)

                                                                                                                   gACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCa                                  <  2:3850864/150‑1 (MQ=255)

                                                                                                                      cGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGt                               <  1:634874/150‑1 (MQ=255)

                                                                                                                      cGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGt                               >  1:3422003/1‑150 (MQ=255)

                                                                                                                         tATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGcc                            <  2:514324/150‑1 (MQ=255)

                                                                                                                         tATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGcc                            >  1:5515352/1‑150 (MQ=255)

                                                                                                                           tGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGc                          >  2:1500583/1‑150 (MQ=255)

                                                                                                                             cAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCAt                        <  2:2404184/150‑1 (MQ=255)

                                                                                                                             cAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCAt                        >  1:560529/1‑150 (MQ=255)

                                                                                                                                cTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAg                     <  1:709784/150‑1 (MQ=255)

                                                                                                                                    tCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGtct                 >  2:473932/1‑150 (MQ=255)

                                                                                                                                    tCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGtc                  >  2:2506281/1‑149 (MQ=255)

                                                                                                                                     tgcgggCCTGAAGCAGCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGtctc                >  2:1912574/4‑150 (MQ=255)

                                                                                                                                     ccTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGtctc                <  1:3607902/150‑1 (MQ=255)

                                                                                                                                     ccTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGtctc                <  2:208993/150‑1 (MQ=255)

                                                                                                                                     ccTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGtctc                <  1:262705/150‑1 (MQ=255)

                                                                                                                                     ccTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGtctc                >  2:2203839/1‑150 (MQ=255)

                                                                                                                                       cgggCCTGAAGCAGCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGt              >  2:3511062/2‑150 (MQ=255)

                                                                                                                                        gggCCTGAAGCAGCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGtt             >  2:5504750/1‑150 (MQ=255)

                                                                                                                                        gggCCTGAAGCAGCAGGCCCACATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGtt             >  1:1074506/1‑150 (MQ=255)

                                                                                                                                           ccTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCa          <  1:889636/150‑1 (MQ=255)

                                                                                                                                            ctgctgCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAg         <  1:2758117/150‑1 (MQ=255)

                                                                                                                                                tgctTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAGCTGc     >  2:2953711/1‑150 (MQ=255)

                                                                                                                                                   ttCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAGCTGcccc  <  1:586860/150‑1 (MQ=255)

AAGCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAGCTGCCCC  >  NC_002947/4602603‑4602899

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