BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,875,613	T→C	6.8%	V143A (GTG→GCG) ‡	cigR →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,875,613	0	T	C	6.8%	514.4 / 17.9	163	V143A (GTG→GCG) ‡	cigR	membrane protein
Reads supporting (aligned to +/- strand): ref base T (69/82); new base C (7/4); total (77/86)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 166 total aligned reads displayed.

GCCTCCAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCG  >  NC_002947/2875464‑2875762

gCCTCCAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTAGGAGTGGATGCAAGCGGGCGCAGACCTGAGTCTGGTCGCCGt                                                                                                                                                       >  2:4133891/1‑150 (MQ=255)

 ccTCCAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTg                                                                                                                                                      <  2:423171/150‑1 (MQ=255)

 ccTCCAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTg                                                                                                                                                      <  2:786581/150‑1 (MQ=255)

     cAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGa                                                                                                                                                  >  2:3682429/1‑150 (MQ=255)

     cAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGa                                                                                                                                                  >  1:1389374/1‑150 (MQ=255)

     cAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGa                                                                                                                                                  <  2:2224525/150‑1 (MQ=255)

             cAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATCCCGGTCGCCACGGCGACCAGAatc                                                                                                                                          <  1:957280/150‑1 (MQ=255)

              agaagaACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGatca                                                                                                                                         >  1:803451/1‑150 (MQ=255)

               gaagaaCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATCCCGGTCGCCACGGCGACCAGAatcag                                                                                                                                        <  2:298116/150‑2 (MQ=255)

               gaagaaCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATCCCGGTCGCCACGGCGACCAGAatcag                                                                                                                                        <  2:4001653/150‑2 (MQ=255)

                 agaaCCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCt                                                                                                                                      <  2:763669/150‑1 (MQ=255)

                      cTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAg                                                                                                                                 <  2:3339205/150‑1 (MQ=255)

                      cTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAg                                                                                                                                 >  2:2354257/1‑150 (MQ=255)

                        ggCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGgt                                                                                                                               >  1:1847681/1‑150 (MQ=255)

                        ggCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGgt                                                                                                                               >  1:1847702/1‑150 (MQ=255)

                         gCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGgtg                                                                                                                              <  1:1276019/150‑1 (MQ=255)

                                gCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACg                                                                                                                       >  1:2927844/1‑150 (MQ=255)

                                gCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACg                                                                                                                       >  1:1829409/1‑150 (MQ=255)

                                 cAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACgg                                                                                                                      <  2:2565052/150‑1 (MQ=255)

                                  aaGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGTCGACCGGGATCATCTATGAGGTGTTGAACggg                                                                                                                     <  2:1579082/150‑1 (MQ=255)

                                  aaGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACggg                                                                                                                     <  2:1866922/150‑1 (MQ=255)

                                   aGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACgggg                                                                                                                    <  1:3704411/150‑1 (MQ=255)

                                       ccTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCAtt                                                                                                                >  1:2688461/1‑150 (MQ=255)

                                          gCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGa                                                                                                             <  1:38085/150‑1 (MQ=255)

                                          gCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGa                                                                                                             <  1:1054819/150‑1 (MQ=255)

                                          gCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGa                                                                                                             <  2:996397/150‑1 (MQ=255)

                                          gCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGa                                                                                                             <  1:2575726/150‑1 (MQ=255)

                                          gCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGa                                                                                                             <  1:1417987/150‑1 (MQ=255)

                                             aCCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACAACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTa                                                                                                          >  1:654062/1‑150 (MQ=255)

                                              ccTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTaa                                                                                                         <  2:2296134/150‑1 (MQ=255)

                                              ccTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTaa                                                                                                         <  2:100592/150‑1 (MQ=255)

                                              ccTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTaa                                                                                                         <  1:3954376/150‑1 (MQ=255)

                                               cTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAg                                                                                                        <  1:760729/150‑1 (MQ=255)

                                                tGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAgg                                                                                                       <  1:3609877/150‑1 (MQ=255)

                                                 ggTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGc                                                                                                      <  2:572626/150‑1 (MQ=255)

                                                 ggTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGc                                                                                                      <  1:404360/150‑1 (MQ=255)

                                                  gTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGcc                                                                                                     >  2:2307675/1‑150 (MQ=255)

                                                       gCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAgg                                                                                                >  2:2551936/1‑150 (MQ=255)

                                                       gCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAgg                                                                                                >  1:4025931/1‑150 (MQ=255)

                                                       gCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGCATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAgg                                                                                                >  2:1761711/1‑150 (MQ=255)

                                                          aagaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGc                                                                                             <  1:3671812/150‑1 (MQ=255)

                                                          aagaagCTGGATGGTCGTCTGTTGGGGCAGCGGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGc                                                                                             <  1:2714330/150‑1 (MQ=255)

                                                            gaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGTGCGCAGACCTGATTCTGGTCGCCTTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGa                                                                                           <  2:32655/150‑1 (MQ=255)

                                                            gaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGa                                                                                           <  1:4235419/150‑1 (MQ=255)

                                                            gaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGa                                                                                           <  2:1736321/150‑1 (MQ=255)

                                                            gaagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGa                                                                                           <  1:1281538/150‑1 (MQ=255)

                                                             aagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAg                                                                                          <  2:1764393/150‑1 (MQ=255)

                                                             aagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAg                                                                                          <  2:1777145/150‑1 (MQ=255)

                                                             aagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAg                                                                                          <  2:1245463/150‑1 (MQ=255)

                                                             aagCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGGGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATTAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAg                                                                                          <  1:549278/150‑1 (MQ=255)

                                                                   gATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTcc                                                                                    >  1:170053/1‑150 (MQ=255)

                                                                      ggTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCt                                                                                 <  1:2749540/150‑1 (MQ=255)

                                                                       gtcgtcTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAATGCCGCAGGTGCGAGCAGTCCGCTc                                                                                <  2:3791877/150‑1 (MQ=255)

                                                                           tcTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCGGACTGCTCgcac                                                                            <  2:2853996/150‑4 (MQ=255)

                                                                           tcTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTc                                                                            >  2:2494714/1‑150 (MQ=255)

                                                                            cTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAg                                                                                          >  1:3236741/1‑135 (MQ=255)

                                                                            cTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAg                                                                                          <  2:3236741/135‑1 (MQ=255)

                                                                              gTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCt                                                                         >  2:1607559/1‑150 (MQ=255)

                                                                              gTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCt                                                                         >  2:1607596/1‑150 (MQ=255)

                                                                                     cAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGgcg                                                                  <  2:664397/150‑1 (MQ=255)

                                                                                      aGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGgcgc                                                                 >  2:4035960/1‑150 (MQ=255)

                                                                                       gCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCt                                                                >  1:4338733/1‑150 (MQ=255)

                                                                                        cggCCTCACTACGACGGGTATGAGTGGCTGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTg                                                               <  1:163633/148‑1 (MQ=255)

                                                                                            cTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGt                                                           >  1:3091282/1‑150 (MQ=255)

                                                                                             tCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTg                                                          >  2:2286532/1‑150 (MQ=255)

                                                                                                  acgacgGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGc                                                     <  2:505638/150‑1 (MQ=255)

                                                                                                  acgacgGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGc                                                     <  1:2261093/150‑1 (MQ=255)

                                                                                                      cgGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATg                                                 <  2:855947/150‑1 (MQ=255)

                                                                                                      cgGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATg                                                 <  2:3483198/150‑1 (MQ=255)

                                                                                                         gTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGc                                              >  2:2008106/1‑150 (MQ=255)

                                                                                                             gAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTgcg                                          >  1:3119109/1‑150 (MQ=255)

                                                                                                               gTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCgca                                        >  2:2623932/1‑150 (MQ=255)

                                                                                                                 ggATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGTTGGGCTgctgcgc                                      <  2:2305021/150‑6 (MQ=255)

                                                                                                                  gATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCgcagca                                     >  2:543980/1‑150 (MQ=255)

                                                                                                                    tGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCagcagc                                   <  2:3904250/150‑1 (MQ=255)

                                                                                                                       aaGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCTGCGCAGCCCa                                <  1:1735345/150‑1 (MQ=255)

                                                                                                                       aaGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCa                                >  2:3871321/1‑150 (MQ=255)

                                                                                                                         gCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAAc                              >  1:639535/1‑150 (MQ=255)

                                                                                                                         gCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCGATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCGCGGTCCCTGCAGGCGCTGGCGTGCCTGCGCTGGGCTGCGCAGCAGcccgcc                              >  1:1068062/1‑147 (MQ=255)

                                                                                                                           gggCGCAGACCTTATTCTGGTCGCGGTGGCGACGGGTATCTTTTATGAGGTGTTGAACGGGGCTTTCGACTAAGGCCGCAGTTGCGAGCAGTCCGCTGGGTCCCTGCTGGCGCTGGGGTGCCTGCTTTGGGCTGCGCAGCAGCCCAACgg                            <  2:1236319/150‑1 (MQ=255)

                                                                                                                           gggCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTTTTGAACGGGGCATTCGACTAAGGCCGCTGGTGCGAGCTGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACgg                            <  2:4018752/150‑1 (MQ=255)

                                                                                                                            ggCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGTTGGGCTGCTGCGCAGCCCATCgca                           <  1:2044725/150‑3 (MQ=255)

                                                                                                                               tcatagatGATCCCGGTCGCCACGGCGACCAGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTa                        >  1:2602281/8‑150 (MQ=255)

                                                                                                                               gcAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGc                                              >  2:1677739/1‑128 (MQ=255)

                                                                                                                               gcAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTa                        >  1:477562/1‑150 (MQ=255)

                                                                                                                                 aGACCTTTTTCTGGTCGCGGTGGCGACGGGTATAATCTTTGTGGTTTTTAAGGGGGATTTCTACTAAGGCCGCTGGTGCGAGCAGTCCGCTCGGTCCGTGCTGGCGGTGGGGTGCCTGCTTTGGGCTGCGCAGCAGCCCAACGGCCTATc                      <  2:3584179/150‑1 (MQ=255)

                                                                                                                                 aGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATc                      >  1:213453/1‑150 (MQ=255)

                                                                                                                                 aGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATc                      >  1:2317929/1‑150 (MQ=255)

                                                                                                                                     cTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACg                  >  1:2808182/1‑150 (MQ=255)

                                                                                                                                     cTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACg                  >  1:3097020/1‑150 (MQ=255)

                                                                                                                                         ttCTGGTCGCCGTGGCGACCGGGGTCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGc              >  1:1983563/1‑150 (MQ=255)

                                                                                                                                         tcccggTCGCCACGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCCGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGc              >  1:2837089/5‑150 (MQ=255)

                                                                                                                                         tcccggTCGCCACGGCGACCAGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGc              >  2:921121/5‑150 (MQ=255)

                                                                                                                                         tcccggTCGCCACGGCGACCAGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGc              >  2:921172/5‑150 (MQ=255)

                                                                                                                                           ccggTCGCCACGGCGACCAGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCa                                                                                                  >  2:2723246/3‑64 (MQ=37)

                                                                                                                                            tGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGt           >  2:1762523/1‑150 (MQ=255)

                                                                                                                                            cggTCGCCACGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGt           >  1:1293362/2‑150 (MQ=255)

                                                                                                                                                    gTGGCGCCCGGTATCATCTTTGAGGGGTAGAACGGGGCTTTCCCCGTAGTCCGCAGGTGCGAGCACGCCGCCCGGTCCCTGCTGGCGCGGGGCTGCCTGCGTTGGGCTGCGCTGCAGCCCCACGGCCTATCAACGCTGCCGTAACACcgc   <  2:2955417/150‑1 (MQ=255)

                                                                                                                                                    gTGGCGACCGGGTTCATCTTTGAGGTGTTGACCGGGGCATTGGACTACGGCCGCCGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGGGTGCCTGGGATGGGGTGCGCAGCTGCCCAACGGCCTATCAACGCTGCCGTAACACcgc   <  2:2955418/150‑1 (MQ=255)

                                                                                                                                                     tGGCGACCGTTATAATCTATGTGTTTTTGAAGGGGGAATTCTACTAAGGCCGCTGGTGCGAGAAGTCCGCTGGGTCCCTGATGGGGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACcgcg  <  2:2688461/150‑1 (MQ=255)

GCCTCCAGGCATTCAGAAGAACCTGGCCCGCGGCAAGCCCCTGCCACCTGGTATCGCCAAGAAGCTGGATGGTCGTCTGTTGGGGCAGCTGCCTCACTACGACGGGTATGAGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCG  >  NC_002947/2875464‑2875762

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