Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,709,217	C→G	10.6%	E100Q (GAG→CAG)	PP_4166 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,709,217	0	C	G	10.6%	194.2 / 11.8	66	E100Q (GAG→CAG)	PP_4166	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/27);  new base G (4/3);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAGGG  >  NC_002947/4709069‑4709356                                                                                                                                                     |                                                                                                                                            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                                                                                                                          <  5:124241/149‑1 (MQ=255)                    tCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTcc                                                                                                                          >  6:237317/1‑149 (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                                                                                             >  3:351965/1‑146 (MQ=255)                                                    cATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGc                                                                                             >  6:346091/1‑146 (MQ=255)                                                     aTCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                         >  5:46722/1‑149 (MQ=255)                                                     aTCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                         >  1:357829/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                                              <  3:332111/144‑1 (MQ=255)                                                                                                     accgacGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCt                                              >  4:332111/1‑144 (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)                                                                                                                                                     |                                                                                                                                            TGATGGCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAGGG  >  NC_002947/4709069‑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: ‑