Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,686,434	A→C	7.7%	D281A (GAC→GCC)	PP_3251 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,686,434	0	A	C	7.7%	225.6 / 8.4	78	D281A (GAC→GCC)	PP_3251	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base C (3/3);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

CCGCCAGCGGCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAG  >  NC_002947/3686287‑3686582                                                                                                                                                      |                                                                                                                                                     ccgccAGCGGCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGAc                                                                                                                                                     <  4:138218/149‑1 (MQ=255)           ccACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGc                                                                                                                                            >  5:369913/1‑148 (MQ=255)             aCGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                         >  7:286755/1‑149 (MQ=255)             aCGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                         <  4:305121/149‑1 (MQ=255)             aCGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                         >  3:255009/1‑149 (MQ=255)                cccgccTGGCAATTCAGGTGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGCCTTGCAGGTCGAACTg                                                                                                                                      <  8:286517/149‑1 (MQ=255)                 ccgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGt                                                                                                                                     >  1:216328/1‑149 (MQ=255)                  cgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGCCTTGCAGGTCGAACTgg                                                                                                                                     >  4:406548/1‑147 (MQ=255)                  cgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGCCTTGCAGGTCGAACTgg                                                                                                                                     <  3:406548/148‑2 (MQ=255)                     cTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCaa                                                                                                                                               >  8:284905/1‑135 (MQ=255)                     cTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTg                                                                                                                                  >  3:244780/1‑148 (MQ=255)                          aaTTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGCCTTGCAGGTCGAACTgg                                                                                                                                     >  7:333440/1‑139 (MQ=255)                          aaTTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGCCTTGCAGGTCGAACTgg                                                                                                                                     <  8:333440/140‑2 (MQ=255)                          aaTTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGccc                                                                                                                            >  1:100903/1‑149 (MQ=255)                          aaTTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGACAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGccc                                                                                                                            >  7:113179/1‑149 (MQ=255)                              cAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTg                                                                                                                         >  4:276207/1‑148 (MQ=255)                               aGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                         >  7:171932/1‑131 (MQ=255)                               aGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                         <  8:171932/131‑1 (MQ=255)                                   gcTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGtt                                                                                                                                    >  6:81843/1‑132 (MQ=255)                                   gcTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGtt                                                                                                                                    <  5:81843/132‑1 (MQ=255)                                    cTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCaa                                                                                                                  <  8:286755/149‑1 (MQ=255)                                            gCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTc                                                                                                           >  1:390469/1‑148 (MQ=255)                                                  ccccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGcc                                                                                                    <  1:271045/149‑1 (MQ=255)                                                  ccccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGc                                                                                                     >  7:52120/1‑148 (MQ=255)                                                  ccccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTCTTTGCAGCCCTTGAACGCAAGCGCGTCTATGGc                                                                                                     >  5:389699/1‑148 (MQ=255)                                                     cgcgCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAAt                                                                                                 <  3:276207/149‑1 (MQ=255)                                                         cATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCTTCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTg                                                                                               <  1:328592/149‑1 (MQ=255)                                                         cATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGat                                                                                             <  6:369913/149‑1 (MQ=255)                                                          aTGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGatg                                                                                            <  8:52120/149‑1 (MQ=255)                                                                   aGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGAc                                                                                   >  5:228470/1‑149 (MQ=255)                                                                       gTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCAtt                                                                               <  6:278556/149‑1 (MQ=255)                                                                        ttCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTg                                                                              <  5:365432/149‑1 (MQ=255)                                                                        ttCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑AGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTg                                                                              >  2:186178/1‑149 (MQ=255)                                                                           tACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCAtt                                                                               <  1:389072/145‑1 (MQ=255)                                                                           tACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCAtt                                                                               >  2:389072/1‑145 (MQ=255)                                                                            aCGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTgccgc                                                                          >  4:277732/1‑149 (MQ=255)                                                                              gCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCt                                                                        >  4:391765/1‑149 (MQ=255)                                                                                 aGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTgca                                                                     <  4:255009/149‑1 (MQ=255)                                                                                    gcgcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTgcagca                                                                  <  6:228470/149‑1 (MQ=255)                                                                                      gcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCagcagc                                                                <  4:37475/149‑1 (MQ=255)                                                                                         cTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCa                                                             <  4:341517/149‑1 (MQ=255)                                                                                             ggCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGa                                                          <  3:391765/148‑1 (MQ=255)                                                                                             ggCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGa                                                          <  4:244780/148‑1 (MQ=255)                                                                                                tACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGatcaa                                                      >  5:191859/1‑149 (MQ=255)                                                                                                  ccGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGatcaatc                                                    <  8:113179/149‑1 (MQ=255)                                                                                                    gCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGAtcaatcaa                                                  <  2:390469/149‑1 (MQ=255)                                                                                                    gCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTCTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGAtcaatcaa                                                  >  8:140338/1‑149 (MQ=255)                                                                                                     ccTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATcaatcaat                                                 <  6:191859/149‑1 (MQ=255)                                                                                                         aaTGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtgg                                               <  6:6788/147‑1 (MQ=255)                                                                                                         aaTGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtgg                                               <  7:232436/147‑1 (MQ=255)                                                                                                           tGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtggctgg                                           <  2:100903/149‑1 (MQ=255)                                                                                                             aCCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATcaatcaat                                                 >  6:106828/1‑141 (MQ=255)                                                                                                             aCCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATcaatcaat                                                 <  5:106828/141‑1 (MQ=255)                                                                                                                ggggTATCGCGACGCCAGAGCTCATC‑‑CGGCAGTTCGACCTGCTATGCGAACTGTTTGCAGCCCTTGAACGCATGCGGGTCTATGTCCAATTGATGTATGCCGACCATTGCCGCCTgctt                                                                    <  6:386589/119‑3 (MQ=255)                                                                                                                ggggTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGTATGCCGACCATTGCCGCCTgcag                                                                    >  5:386589/1‑119 (MQ=255)                                                                                                                ggggTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGaa                                      <  2:216328/149‑1 (MQ=255)                                                                                                                        gccacgccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGgcgc                              >  2:379397/1‑149 (MQ=255)                                                                                                                           acgccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCAt                           >  2:131301/1‑149 (MQ=255)                                                                                                                           acgccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCAt                           >  7:257789/1‑149 (MQ=255)                                                                                                                             gccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtggc                                              <  3:172707/128‑1 (MQ=255)                                                                                                                             gccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGc                         <  2:377880/149‑1 (MQ=255)                                                                                                                              ccaGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGcc                        >  3:77419/1‑149 (MQ=255)                                                                                                                               caGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGccc                       <  6:389699/149‑1 (MQ=255)                                                                                                                                  aCCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCa                    >  4:390379/1‑149 (MQ=255)                                                                                                                                      cATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCg                >  2:345176/1‑149 (MQ=255)                                                                                                                                        tC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCACCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGtt              >  5:368739/1‑149 (MQ=255)                                                                                                                                           caaacAGTTCGCCTCGCAGGTCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTg             >  4:140568/5‑149 (MQ=255)                                                                                                                                             gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCATTTGATGCATGCCGCCCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCTTTCGTTGcc           >  6:115247/1‑149 (MQ=255)                                                                                                                                             gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCagcagc                                                                >  6:374059/1‑96 (MQ=255)                                                                                                                                             gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCagcagc                                                                <  5:374059/96‑1 (MQ=255)                                                                                                                                             gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGcc           >  3:401756/1‑149 (MQ=255)                                                                                                                                              ccAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCg          >  8:348239/1‑149 (MQ=255)                                                                                                                                                 gTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCAt       >  1:30771/1‑149 (MQ=255)                                                                                                                                                  ttCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCAtt      <  3:277732/149‑1 (MQ=255)                                                                                                                                                  ttCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCAt       >  3:86203/1‑148 (MQ=255)                                                                                                                                                  ttCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCAt       >  3:80643/1‑148 (MQ=255)                                                                                                                                                     gACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGc            <  7:45356/140‑1 (MQ=255)                                                                                                                                                     gACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGc            >  8:45356/1‑140 (MQ=255)                                                                                                                                                      aCCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAg  >  2:187764/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCGCCAGCGGCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATC‑‑CGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAG  >  NC_002947/3686287‑3686582

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