Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,709,214	C→T	10.2%	A101T (GCG→ACG)	PP_4166 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,709,214	0	C	T	10.2%	186.7 / 7.8	68	A101T (GCG→ACG)	PP_4166	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base T (4/3);  total (37/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 = 1.00e+00

GAATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAGGG  >  NC_002947/4709066‑4709356                                                                                                                                                     |                                                                                                                                               gaATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAgcgc                                                                                                                                                >  7:245555/1‑149 (MQ=255)   atgatgGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCt                                                                                                                                              >  1:243904/1‑149 (MQ=255)    tgatgGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTc                                                                                                                                             <  1:149762/149‑1 (MQ=255)         gCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCggg                                                                                                                                        >  2:297645/1‑149 (MQ=255)         gCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCgg                                                                                                                                         <  8:276802/148‑1 (MQ=255)           ccATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCgg                                                                                                                                         >  7:168510/1‑146 (MQ=255)           ccATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCgg                                                                                                                                         <  8:168510/146‑1 (MQ=255)                   ggCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTg                                                                                                                              <  1:297645/149‑1 (MQ=255)                    gCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGTCTGGCCCGCGAGGCGCTgg                                                                                                                             <  5:43800/149‑1 (MQ=255)                       tCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGTCTGGCCCGCGAGGCGCTggaga                                                                                                                          <  7:132650/149‑4 (MQ=255)                       tCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGTCTGGCCCGCGAGGCGCTggaga                                                                                                                          <  2:124355/149‑4 (MQ=255)                       tCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTcc                                                                                                                          >  6:237317/1‑149 (MQ=255)                       tCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTcc                                                                                                                          <  5:124241/149‑1 (MQ=255)                          cgGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTgg                                                                                                                       <  1:44684/149‑1 (MQ=255)                            gggCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTa                                                                                                                     <  3:170742/149‑1 (MQ=255)                                aGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTc                                                                                                                  >  4:364375/1‑148 (MQ=255)                                 gTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGCGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGCCGCGGGTCCTGGTAGTc                                                                                                                  >  4:308286/1‑147 (MQ=255)                                  ttCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAg                                                                                                               >  2:110469/1‑149 (MQ=255)                                        aaTGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGt                                                                                                              >  5:390084/1‑144 (MQ=255)                                        aaTGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGt                                                                                                              <  6:390084/144‑1 (MQ=255)                                         aTGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAgg                                                                                                        <  8:257851/149‑1 (MQ=255)                                               aTAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGGTGGTCCTGGTAGTCCAGTTGCAGgcggcg                                                                                                   >  6:380188/1‑148 (MQ=255)                                               aTAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTgg                                                                                                                       <  4:282451/128‑1 (MQ=255)                                               aTAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTgg                                                                                                                       >  3:282451/1‑128 (MQ=255)                                               aTAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTc                                                                                                                  <  1:304576/133‑1 (MQ=255)                                               aTAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACAGACGCAGGCAGGCCTGCATGTAGGCAACTTGGCTCTCCAGCGCCTCGAGGGACAGACGCTGGTCCTGGTAGTc                                                                                                                  >  2:304576/1‑133 (MQ=255)                                                       cATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGc                                                                                             >  6:346091/1‑146 (MQ=255)                                                       cATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGc                                                                                             >  3:351965/1‑146 (MQ=255)                                                        aTCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                         >  1:357829/1‑149 (MQ=255)                                                        aTCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                         >  5:46722/1‑149 (MQ=255)                                                          cGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGAc                                                                                        >  3:20394/1‑148 (MQ=255)                                                                  gggATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGttt                                                                               <  1:110469/149‑1 (MQ=255)                                                                   ggATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTg                                                                              <  6:46722/149‑1 (MQ=255)                                                                         ggTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGt                                                                                 <  5:376213/140‑1 (MQ=255)                                                                         ggTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGt                                                                                 >  6:376213/1‑140 (MQ=255)                                                                          gTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                         >  3:92347/1‑131 (MQ=255)                                                                          gTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                         <  4:92347/131‑1 (MQ=255)                                                                          gTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAgg                                                                       >  8:282697/1‑149 (MQ=255)                                                                             aTGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGc                                                                    <  8:245555/149‑1 (MQ=255)                                                                                ccGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCga                                                                  <  2:243904/148‑1 (MQ=255)                                                                                      ggTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                         >  7:94700/1‑119 (MQ=255)                                                                                      ggTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                         <  8:94700/119‑1 (MQ=255)                                                                                        tCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGGATGCAGCAGGTGCGAGACCCgcag                                                         >  6:77525/1‑149 (MQ=255)                                                                                             aCTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGcaggcagg                                                    >  6:241747/1‑149 (MQ=255)                                                                                                  tCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATcc                                               >  2:337097/1‑149 (MQ=255)                                                                                                  tCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATcc                                               >  1:366302/1‑149 (MQ=255)                                                                                                    aGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTg                                             <  7:282697/149‑1 (MQ=255)                                                                                                        ctacatgcaggcctGCCTGCGTGTAGGCCACTTGGCTCTCCAGCGCATCGCCGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGc                                         >  4:192233/7‑149 (MQ=255)                                                                                                        ctacatgcaggcctGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGc                                         >  1:163585/7‑149 (MQ=255)                                                                                                        accgacGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCt                                              >  4:332111/1‑144 (MQ=255)                                                                                                        accgacGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCt                                              <  3:332111/144‑1 (MQ=255)                                                                                                          cgacGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCCTCCTGCAGCg                                        >  6:323042/1‑148 (MQ=255)                                                                                                          cgacGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCg                                        >  3:339911/1‑148 (MQ=255)                                                                                                                 ggcaggcCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAgg                                <  5:346091/149‑1 (MQ=255)                                                                                                                       ccTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGc                          >  3:124619/1‑149 (MQ=255)                                                                                                                                aGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgt                 >  2:146774/1‑149 (MQ=255)                                                                                                                                 ggCCACTTTGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGTTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgtt                <  1:134905/149‑1 (MQ=255)                                                                                                                                 ggCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgtt                <  5:237317/149‑1 (MQ=255)                                                                                                                                 ggCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgtt                <  4:84007/149‑1 (MQ=255)                                                                                                                                 ggCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgtt                <  4:351965/149‑1 (MQ=255)                                                                                                                                  gCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgttg               >  8:87009/1‑149 (MQ=255)                                                                                                                                    cACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGcc             >  7:9968/1‑149 (MQ=255)                                                                                                                                     aCTTGGGTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCt            <  1:268894/149‑1 (MQ=255)                                                                                                                                         ggCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCg          <  4:20394/147‑1 (MQ=255)                                                                                                                                         ccaggaccAGCGTCTGGCCCGCGAGGCGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCg          >  8:206140/7‑147 (MQ=255)                                                                                                                                           aggaccAGCGTCTGGCCCGCGAGGCGCTGGTCCTGGTAGTCCAGATGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGCCCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGGTCGCATTGTTGCCTCgctc       >  8:169009/5‑146 (MQ=255)                                                                                                                                               ccAGCGTCTGGCCCGCGAGGCGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAggg  >  4:201110/1‑149 (MQ=255)                                                                                                                                               ccAGCGTCTGGCCCGCGAGGCGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAggg  >  1:310520/1‑149 (MQ=255)                                                                                                                                                cAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAgg   <  2:163585/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               GAATGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAGGG  >  NC_002947/4709066‑4709356

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