Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,101,079	T→C	6.0%	V127A (GTC→GCC)	PP_2710 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,101,079	0	T	C	6.0%	196.1 / 2.4	66	V127A (GTC→GCC)	PP_2710	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base C (2/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

GGAGGAATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAGCACAT  >  NC_002947/3100933‑3101226                                                                                                                                                   |                                                                                                                                                     ggaggaATACGGCAACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTCGGAGGCCATGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCg                                                                                                                                                    >  7:103087/1‑149 (MQ=255)  gaggaATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCgg                                                                                                                                                   <  4:241053/149‑1 (MQ=255)  gaggaATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCg                                                                                                                                                    <  3:171963/148‑1 (MQ=255)   aggaATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGTAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGACCggg                                                                                                                                                  >  4:64946/1‑149 (MQ=255)   aggaATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCggg                                                                                                                                                  <  3:75297/149‑1 (MQ=255)   aggaATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCggg                                                                                                                                                  <  2:262444/149‑1 (MQ=255)     gaATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGgtt                                                                                                                                                >  1:34224/1‑149 (MQ=255)      aaTACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGgttg                                                                                                                                               <  3:61430/149‑1 (MQ=255)      aaTACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGACCGGgttg                                                                                                                                               <  7:317789/149‑1 (MQ=255)       aTACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCATGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGgttgt                                                                                                                                              <  8:103087/149‑1 (MQ=255)          cGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGTAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGACCGGgttgcca                                                                                                                                            <  3:64946/148‑4 (MQ=255)            gCCACGAGCAGTGGATCCGCAACGAGATCAGAGCGTGCGGCGGTGACTGCGAGGCCGTGCGCGATGGCCAGCCGGTGGAATCCATCGCGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGAGGCAAGCCGGTCCGGTGGATCgg                                                                                                                                          >  8:80612/1‑148 (MQ=255)                      gTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACg                                                                                                                               >  7:115655/1‑149 (MQ=255)                            ccTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTc                                                                                                                          >  7:160537/1‑148 (MQ=255)                            ccTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTc                                                                                                                          >  4:339308/1‑148 (MQ=255)                              tCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCCGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑g                                                                                                                        >  8:335448/1‑147 (MQ=255)                                aaCGACATCAGAGCTTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACTACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑Gaa                                                                                                                      <  8:175806/147‑1 (MQ=255)                                aaCGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAg                                                                                                                     <  4:309689/148‑1 (MQ=255)                                aaCGACATCAGAGCGTGCGGCGGGGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCTATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACGCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAg                                                                                                                     <  6:164606/148‑1 (MQ=255)                                 aCGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAgg                                                                                                                    >  5:310366/1‑148 (MQ=255)                                    aCATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCAcc                                                                                                                <  4:116258/149‑1 (MQ=255)                                     cATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCa                                                                                                               <  1:374305/149‑1 (MQ=255)                                          gagCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATc                                                                                                          <  1:94248/149‑1 (MQ=255)                                             cGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGcc                                                                                                       >  6:210559/1‑149 (MQ=255)                                                     ggTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCAc                                                                                               <  6:74814/149‑1 (MQ=255)                                                              gaccatcagctcgatggattgcaccggctggccatcgcgcTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAggcggc                                                                                      >  2:330386‑M2/41‑149 (MQ=255)                                                              gaccatcagctcgatggattgcaccggctggccatcgcgcTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAggcggc                                                                                      >  2:325018‑M2/41‑149 (MQ=255)                                                              gAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTAATGGTCGCCTACCTGTACGACCTGATTGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAggcggc                                                                                      >  1:347035/1‑149 (MQ=255)                                                                  ccGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCGTGGCCACCCAGGCggcgggc                                                                                  <  8:266870/149‑1 (MQ=255)                                                                   cGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCggcgggcg                                                                                 >  7:84515/1‑149 (MQ=255)                                                                     agctcgatggattgcaccggctggccatcgcgcTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAg                                                                                                                     >  2:224864‑M2/34‑111 (MQ=255)                                                                     agctcgatggattgcaccggctggccatcgcgcTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAg                                                                                                                     <  1:224864‑M2/78‑1 (MQ=255)                                                                        tcgatggattgcaccggctggccatcgagcTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCCCCCAGGCGGCGGGCGTTATc                                                                            >  6:261373‑M2/31‑149 (MQ=255)                                                                        gcgATGGCCAGCCGGTGCAATCCATCGAGCCGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATc                                                                            <  7:335448/149‑1 (MQ=255)                                                                         cgATGGCCAGCCGGTGCAATCCATCGTGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATc                                                                            >  1:201763/1‑148 (MQ=255)                                                                         cgATGGCCAGCCGGTGCAATCCATCGGGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATc                                                                            >  1:195718/1‑148 (MQ=255)                                                                          gATGGCCAGCCGGTGGAATCCATCGAGCTGTTGGTCGCCTACCTGTACGACCTGATCGGCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATGGCCCTGTCCACCCAGGCGGCGGGCGTTATCCa                                                                          <  7:80612/149‑1 (MQ=255)                                                                              gCCAGCCGGTGCAATCCATCGAGCTAATGGTCGCCTACCTGTACGACCTGATTGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGaa                                                                       <  2:347035/148‑1 (MQ=255)                                                                                cAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAaca                                                                     >  2:168784/1‑148 (MQ=255)                                                                                aaGCCGGTGCAATCAATCGAGCTGATTGAAGCCTAACTGTACGACATGATCGACCGTTGCAACACTGTAGGGTTGTTAGGCATGGTAAACGTGCTC‑GAAGTCACAAGCATCGCACTGGCCACCCAGGCGGCGGGCGTTATCCAGAaca                                                                     >  8:244223/2‑148 (MQ=255)                                                                                    cGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTg                                                                <  2:34224/149‑1 (MQ=255)                                                                                      gTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGc                                                              >  6:217199/1‑149 (MQ=255)                                                                                         cAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCcct                                                           <  3:339308/149‑1 (MQ=255)                                                                                                tCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGAc                                                    <  1:168784/149‑1 (MQ=255)                                                                                                tCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGACCGTGGCAACCAGGTCGGGTTGTTCGGCATTGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCTGGCGTTATCCAGAACACATTGGCCCTGCCTGAc                                                    <  7:244223/149‑1 (MQ=255)                                                                                                   aGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAg                                                 <  8:115655/149‑1 (MQ=255)                                                                                                     cTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGc                                               <  5:217199‑M2/148‑1 (MQ=255)                                                                                                       gATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCCGAACACATTGGCCCTGCCTGACCAGGCGt                                             >  3:4428/1‑149 (MQ=255)                                                                                                        aTGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCcctgcct                                                       <  1:263364/138‑1 (MQ=255)                                                                                                        aTGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCcctgcct                                                       >  2:263364/1‑138 (MQ=255)                                                                                                        aTGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGGTCAGGcag                                             <  6:310366/148‑3 (MQ=255)                                                                                                           gTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAg                                         <  5:222696/149‑1 (MQ=255)                                                                                                               ccTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTAc                                     >  1:128111/1‑149 (MQ=255)                                                                                                               ccTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTAc                                     <  6:220350/149‑1 (MQ=255)                                                                                                                     tGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGc                               >  3:363751/1‑149 (MQ=255)                                                                                                                          gACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCa                          >  6:271459/1‑149 (MQ=255)                                                                                                                          gACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCa                          >  5:356289/1‑149 (MQ=255)                                                                                                                           aCCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCAc                         >  8:14733/1‑149 (MQ=255)                                                                                                                                     ccGTGGCAACCCGGTGGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGGGTTCAGCTACCTCAGCTCGCCCGGTGcgcgcg               >  8:242859/1‑149 (MQ=255)                                                                                                                                     ccGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTc                >  3:324939/1‑148 (MQ=255)                                                                                                                                       cgaacAACCCGACCGGGTTGTTCGGCATGGTCAACGTGCTCGGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCATGCCTGACCAGGCGTTCAGCTACCTCAGATCGCACGGTGCGCTCGa              >  1:259850/5‑149 (MQ=255)                                                                                                                                          gCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAg          <  5:210559/149‑1 (MQ=255)                                                                                                                                            aaCCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAg          >  3:249703/1‑147 (MQ=255)                                                                                                                                             aCCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAg       >  2:291873/1‑149 (MQ=255)                                                                                                                                              cccGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCAAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTAGCACGGTGCGCTCGACCAGGAg       >  7:345212/1‑148 (MQ=255)                                                                                                                                                  gTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAGCACAt  >  5:238770/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GGAGGAATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTC‑GAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAGCACAT  >  NC_002947/3100933‑3101226

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