Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	874,702	C→A	7.2%	A94D (GCC→GAC) ‡	PP_0759 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	874,702	0	C	A	7.2%	221.1 / 5.9	69	A94D (GCC→GAC) ‡	PP_0759	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base A (2/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

CATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCAGCTCAAGCTCAC  >  NC_002947/874562‑874843                                                                                                                                             |                                                                                                                                              cATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGGGCGATCAGgcgg                                                                                                                                       <  2:126750/149‑2 (MQ=255) cATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGgcgc                                                                                                                                       <  2:20717/149‑1 (MQ=255) cATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGgcgc                                                                                                                                       <  2:255730/149‑1 (MQ=255)        tcGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGGGCGATCAGgcggt                                                                                                                                      <  6:39452/143‑3 (MQ=255)        tcGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGGGCGATCAGgcggt                                                                                                                                      >  5:39452/1‑141 (MQ=255)             cGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGcc                                                                                                                           >  7:9177/1‑149 (MQ=255)              gATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGccc                                                                                                                          >  1:78246/1‑149 (MQ=255)              gATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGccc                                                                                                                          <  2:78246/149‑1 (MQ=255)                  aGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGAt                                                                                                                      <  1:51059/149‑1 (MQ=255)                   ggTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATg                                                                                                                     >  5:114177/1‑149 (MQ=255)                    gTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGccc                                                                                                                                             <  4:177372/124‑1 (MQ=255)                    gTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGccc                                                                                                                                             >  3:177372/1‑124 (MQ=255)                       aaCCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTTTGATCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGACATTGCCAAAACCCTGATCGCCAAGTCGCCCAAATTACGCCCGGATGCagc                                                                                                                 >  3:196598/1‑149 (MQ=255)                       aaCCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCcgc                                                                                                                 >  5:266547/1‑149 (MQ=255)                        aCCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCcgcg                                                                                                                <  3:148179/149‑1 (MQ=255)                              tGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTgg                                                                                                          <  5:288169/149‑1 (MQ=255)                                tACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTggg                                                                                                         <  6:242395/148‑1 (MQ=255)                                     gcAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTgg                                                                                                    >  5:174766/1‑148 (MQ=255)                                     gcAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGc                                                                                                   <  2:191721/149‑1 (MQ=255)                                       aGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGGGCGATCAGgcggttg                                                                                                                                    >  4:273158/1‑110 (MQ=255)                                       aGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGGGCGATCAGgcggttg                                                                                                                                    <  3:273158/114‑5 (MQ=255)                                       aGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGc                                                                                                   >  4:49379/1‑147 (MQ=255)                                       aGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCaa                                                                                                 <  8:9177/149‑1 (MQ=255)                                       aGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCaa                                                                                                 <  4:49192/149‑1 (MQ=255)                                                agcagACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGt                                                                                         <  4:46064/148‑1 (MQ=255)                                                 gcagACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGa                                                                                       >  4:99921/1‑149 (MQ=255)                                                 gcagACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGa                                                                                       <  3:138910/149‑1 (MQ=255)                                                       ccTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATc                                                                                     >  3:154182/1‑145 (MQ=255)                                                       ccTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCaa                                                                                   >  4:282091/1‑147 (MQ=255)                                                           ggTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGtg                                                                              <  6:174766/148‑1 (MQ=255)                                                           ggTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGtg                                                                              <  3:282091/148‑1 (MQ=255)                                                           ggTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGtg                                                                              <  3:49379/148‑1 (MQ=255)                                                            gTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGtgat                                                                            >  8:72987/1‑149 (MQ=255)                                                                 gCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGAtgaact                                                                       >  4:88385/1‑149 (MQ=255)                                                                    tCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGAtgaactgaa                                                                    >  8:196094/1‑149 (MQ=255)                                                                    tCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGAtgaactgaa                                                                    <  6:266547/149‑1 (MQ=255)                                                                       aGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGc                                                                 >  2:273331/1‑149 (MQ=255)                                                                        gCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGcc                                                                <  6:114177/149‑1 (MQ=255)                                                                               tgtgCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCg                                                         >  7:17713/1‑149 (MQ=255)                                                                               tgtgCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCg                                                         >  3:243556/1‑149 (MQ=255)                                                                                    tCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGAATGCCGCGCAGTGGGACTgg                                                                                                    <  7:102088/101‑1 (MQ=255)                                                                                    tCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGAATGCCGCGCAGTGGGACTgg                                                                                                    >  8:102088/1‑101 (MQ=255)                                                                                      gACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCggg                                                  <  4:243556/149‑1 (MQ=255)                                                                                        cAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCCGCCCGggcg                                                >  7:191322/1‑149 (MQ=255)                                                                                           gTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAaga                                                                                 <  1:203405/113‑1 (MQ=255)                                                                                           gTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAaga                                                                                 >  2:203405/1‑113 (MQ=255)                                                                                            tCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGggcgg                                               >  5:138036/1‑146 (MQ=255)                                                                                            tCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGggcgg                                               <  6:138036/146‑1 (MQ=255)                                                                                               aGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGggcgg                                               <  4:195606/143‑1 (MQ=255)                                                                                               aGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGggcgg                                               >  3:195606/1‑143 (MQ=255)                                                                                               aGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCATGa                                          >  2:137627/1‑148 (MQ=255)                                                                                                gcgcCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCACg                                                    >  4:185916/1‑137 (MQ=255)                                                                                                           aaGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTacac                             >  5:128063/1‑149 (MQ=255)                                                                                                           aaGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGCCGCCCAAATTCCGCACGGATGCCGCGCACTGGGACTGGAAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTacac                             >  8:185995/1‑149 (MQ=255)                                                                                                              cGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCgg                          <  3:88385/149‑1 (MQ=255)                                                                                                              cGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCgg                          >  8:219760/1‑149 (MQ=255)                                                                                                               gtgtGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTac                               <  1:303034/143‑1 (MQ=255)                                                                                                               gtgtGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTac                               >  2:303034/1‑143 (MQ=255)                                                                                                                    cAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTgat                    >  5:131503/1‑149 (MQ=255)                                                                                                                    cAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTga                     <  7:105085/148‑1 (MQ=255)                                                                                                                       gCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCCACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGc                                                       >  5:279995/1‑111 (MQ=255)                                                                                                                       gCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGc                                                       <  6:279995/111‑1 (MQ=255)                                                                                                                       gCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGa                                          >  6:259880/1‑124 (MQ=255)                                                                                                                       gCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGa                                          <  5:259880/124‑1 (MQ=255)                                                                                                                       gCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTgatcg                  >  3:90799/1‑148 (MQ=255)                                                                                                                         cATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTgatcgat                >  4:87702/1‑148 (MQ=255)                                                                                                                          aTTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCATGATCATCGTCTACACCGGGCTGATCGATCa              <  1:137627/149‑1 (MQ=255)                                                                                                                              ccAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCagctc          <  4:154182/149‑1 (MQ=255)                                                                                                                                 aCCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCagctcaag       <  1:273331/149‑1 (MQ=255)                                                                                                                                  ccGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCagctcaagc      >  3:159230/1‑149 (MQ=255)                                                                                                                                    gCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCagctcaagct     >  8:131415/1‑148 (MQ=255)                                                                                                                                             cccAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCAGCTCAAGCTCAc  <  4:228361/142‑1 (MQ=255)                                                                                                                                             cccAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCAGCTCAAGCTCAc  >  3:228361/1‑142 (MQ=255)                                                                                                                                             |                                                                                                                                              CATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCAGCTCAAGCTCAC  >  NC_002947/874562‑874843

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