Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,993,334	A→G	7.5%	L374P (CTG→CCG)	PP_2617 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,993,334	0	A	G	7.5%	178.3 / 2.3	63	L374P (CTG→CCG)	PP_2617	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base G (2/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCGCCTTGCTCATGCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCCGGCGGTCTTCCAGCGTGGTCGTCATCACGCCGAATTGT  >  NC_002947/2993189‑2993480                                                                                                                                                  |                                                                                                                                                   aTTCGCCTTGCTCATGCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGa                                                                                                                                                >  1:1080818/1‑150 (MQ=255) aTTCGCCTTGCTCATGCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGa                                                                                                                                                <  2:1382218/150‑1 (MQ=255)     gCCTTGCTCATGCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCa                                                                                                                                            >  2:829641/1‑150 (MQ=255)       cTTGCTCATGCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAg                                                                                                                                          >  1:1260609/1‑150 (MQ=255)           cTCATGCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAg                                                                                                                                      <  2:544380/150‑1 (MQ=255)            tCATGCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAgg                                                                                                                                     >  2:175872/1‑150 (MQ=255)              aTGCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGTTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGt                                                                                                                                         >  2:163292/1‑144 (MQ=255)              aTGCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGt                                                                                                                                         <  1:163292/144‑1 (MQ=255)              aTGCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTg                                                                                                                                   >  2:1886496/1‑150 (MQ=255)                gCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGAt                                                                                                                                 <  2:1080818/150‑1 (MQ=255)                gCACGCCCAGGCCGGGTCACGTCGCTGTGGCACGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGAt                                                                                                                                 >  1:989250/1‑150 (MQ=37)                 cACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATc                                                                                                                                <  1:1886496/150‑1 (MQ=255)                   cGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCtt                                                                                                                              <  1:1252701/150‑1 (MQ=255)                    gCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTg                                                                                                                             <  2:1249677/150‑1 (MQ=255)                    gCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTg                                                                                                                             >  2:631458/1‑150 (MQ=255)                      ccAGGCCGGGTCGCGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGc                                                                                                                           <  2:888772/150‑1 (MQ=255)                           ccGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCt                                                                                                                      <  1:829641/150‑1 (MQ=255)                           ccGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCt                                                                                                                      <  2:1671788/150‑1 (MQ=255)                               gTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGAt                                                                                                                  <  1:1530171/150‑1 (MQ=255)                                        gTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTg                                                                                                         <  1:1877695/150‑1 (MQ=255)                                             gCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCa                                                                                                    >  1:1259497/1‑150 (MQ=255)                                               tCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAAc                                                                                                  >  1:869951/1‑150 (MQ=255)                                                    gTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCtg                                                                                             <  2:1796961/150‑1 (MQ=255)                                                     tAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTgg                                                                                                        >  1:1805246/1‑138 (MQ=255)                                                     tAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTgg                                                                                                        <  2:1805246/138‑1 (MQ=255)                                                         tGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGAc                                                                                                                                               >  2:1411613/1‑95 (MQ=255)                                                         tGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGAc                                                                                                                                               <  1:1411613/95‑1 (MQ=255)                                                         tGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCAGCTCCAGGTCATTCAGTGACTCGGGTCGGAGAATGGTCGCCAGCACGATTGACATCTGGCGCAACGCCTgtggtg                                                                                        >  1:1374867/1‑150 (MQ=19)                                                             aCCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCAt                                                                                    <  1:1596578/150‑1 (MQ=255)                                                                cGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCt                                                                                 >  2:1558500/1‑150 (MQ=255)                                                                 gCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTg                                                                                <  1:1012353/150‑1 (MQ=255)                                                                     ggAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGgcgc                                                                            <  1:1242742/150‑1 (MQ=255)                                                                          ggCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACg                                                                       >  2:1701160/1‑150 (MQ=255)                                                                              ggTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTgtggtg                                                                                        <  2:503865/129‑1 (MQ=255)                                                                              ggTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTgtggtg                                                                                        >  1:503865/1‑129 (MQ=255)                                                                                        cgcCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTTTGGATGAAc                                                                                                              <  2:275381/97‑1 (MQ=255)                                                                                        cgcCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTTTGGATGAAc                                                                                                              >  1:275381/1‑97 (MQ=255)                                                                                            aagggTAATGGCAGTCTAGGGATGTACGGTATTTGCAGTGACAACAGTTCATTCCGCGATGCAAGGTAAGTGAACTGGACTCTCAGTAGGAACATCTGGCGCAACGCCTGTGATGGCATGCTGGCGCGCACGGCGAGGGTGAAGTCGGCt                                                     >  2:846437/3‑150 (MQ=255)                                                                                             cGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTgtggtg                                                                                        >  1:1784180/1‑114 (MQ=255)                                                                                             cGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTgtggtg                                                                                        <  2:1784180/114‑1 (MQ=255)                                                                                               ggTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCGGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAgg                                                                 <  2:1240605/135‑1 (MQ=255)                                                                                               ggTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCGGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAgg                                                                 <  2:1240654/135‑1 (MQ=255)                                                                                               ggTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCGGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAgg                                                                 >  1:1240654/1‑135 (MQ=255)                                                                                               ggTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCGGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAgg                                                                 >  1:1240605/1‑135 (MQ=255)                                                                                                gTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTtcgt                                                 >  1:158156/1‑150 (MQ=255)                                                                                                 tAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTtcgtc                                                <  1:1701160/150‑1 (MQ=255)                                                                                                     ggCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGAta                                            >  2:228541/1‑149 (MQ=255)                                                                                                     ggCAGCAAAGGGATGCCCGGCAATTGCAGCGACACCAGGTCATTCGGCCTGGCTGTCTAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATc                                            <  2:989250/150‑1 (MQ=25)                                                                                                        agcaAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGagcagc                                         >  2:1629942/1‑150 (MQ=255)                                                                                                        agcaAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGc                                         >  2:1742734/1‑150 (MQ=255)                                                                                                           aaaGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGCTGGACGGCGGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCcgg                                      >  1:230907/1‑150 (MQ=255)                                                                                                             aGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCcggcg                                    <  2:1259497/150‑1 (MQ=255)                                                                                                                gATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACGCAAGTGGGGTGATCTTGGCTTTCTGGATGAACGTCTGGCGCAACGCCTGTGGTGGCATGCGGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGAGCAGCCGGCGGTc                                 <  2:564068/150‑1 (MQ=255)                                                                                                                      cGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCTGCCATCGTGCTTCCAg                           >  2:1789248/1‑150 (MQ=37)                                                                                                                       ggCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCCGGCGGTCTTCCAGc                          <  1:631458/150‑1 (MQ=255)                                                                                                                            ttGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCCGGCGGTCTTCCAGCGTGgt                     <  1:1249571/150‑1 (MQ=255)                                                                                                                             tGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCCGGCGGTCTTCCAGCGTGgtc                    <  1:1742734/150‑1 (MQ=255)                                                                                                                             ggcggCGTCACCAGGTCATTCAGTGACTCAAGTGCGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCCGGCGGTCTTCCAGCGTGgtc                    <  2:1260609/146‑1 (MQ=255)                                                                                                                              gCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCCGGCGGTCTTCCAGCGTGgtcg                   >  1:359688/1‑150 (MQ=255)                                                                                                                                       ccAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCCGGCGGTCTTCCAGCGTGGTCGTCATCACGc          <  1:232382/150‑1 (MQ=255)                                                                                                                                           gTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCCGGCGGTCTTCCAGCGTGGTCGTCATCACGCCGaa      <  1:1558500/150‑1 (MQ=255)                                                                                                                                              aTTCAGTGACTCAAGAGAGGGGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCCGGCGGTCTTCCAGCGTGGTCGTCATCACGCCGAATTg   >  2:1633823/1‑150 (MQ=255)                                                                                                                                               ttCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCCGGCGGTCTTCCAGCGTGGTCGTCATCACGCCGAATTGt  >  1:204044/1‑150 (MQ=255)                                                                                                                                                  aGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCcggcg                                    >  1:133311/1‑113 (MQ=255)                                                                                                                                                  aGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCcggcg                                    <  2:133311/113‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   ATTCGCCTTGCTCATGCACGCCCAGGCCGGGTCACGTCGGTCCAGCTCGAAGTAGCTGAAACCCGCATGGAAGGGCAGGTCGCGGGGCGCCACGGGTAATGGCAGCAAAGGGATGCCCGGCAATTGCAGCGGCACCAGGTCATTCAGTGACTCAAGTGAGGTGATCTTGGCTTTCTGGATGAACATCTGGCGCAACGCCTGTGGTGGCATGCTGGCGCGCACGGCGAGGATGAAGTCGGCTTCGTCGATCAGCCGGCGGTCTTCCAGCGTGGTCGTCATCACGCCGAATTGT  >  NC_002947/2993189‑2993480

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