Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,133,521	A→G	8.0%	A36A (GCA→GCG)	PP_0995 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,133,521	0	A	G	8.0%	169.3 / 4.5	62	A36A (GCA→GCG)	PP_0995	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base G (2/3);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAGCGCTATCCAGGCTGATTGGAGGAGCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGTTGATCGGGCTGTGGTTGACGAGCCGCTTGGCCAGGCCAG  >  NC_002947/1133375‑1133666                                                                                                                                                   |                                                                                                                                                  tGAAAGCGCTATCCAGGCTGATTGGAGGAGCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCgcag                                                                                                                                                  >  8:60836/1‑148 (MQ=255)   aaaGCGCTATCCAGGCTGATTGGAGGAGCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCAGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCAc                                                                                                                                               >  8:137059/1‑149 (MQ=255)     aGCGCTATCCAGGCTGATTGGAGGAGCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCActg                                                                                                                                             <  8:16032/149‑1 (MQ=255)     aGCGCTATCCAGGCTGATTGGAGGAGCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCActg                                                                                                                                             >  7:291053/1‑149 (MQ=255)            tCCAGGCTGATTGGAGGAGCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGa                                                                                                                                      <  7:60836/149‑1 (MQ=255)              cAGGCTGATTGGAGGAGCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACt                                                                                                                                    >  2:218068/1‑149 (MQ=255)                   tGATTGGAGGAGCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTAc                                                                                                                               >  1:300564/1‑149 (MQ=255)                        ggaggaGCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCtt                                                                                                                           >  1:157466/1‑148 (MQ=255)                         gaggagCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTa                                                                                                                         >  2:71358/1‑149 (MQ=255)                         gaggagCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTa                                                                                                                         <  3:108370/149‑1 (MQ=255)                           ggagCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACg                                                                                                                       <  2:147830/149‑1 (MQ=255)                           ggagCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACg                                                                                                                       <  4:78114/149‑1 (MQ=255)                            gagCCACCAACCATGAAGCCTGCCGATGTGCTTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGACGCGCCACGCCGTTGGAAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACg                                                                                                                       <  3:207307/148‑1 (MQ=255)                            gagCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTCGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTAACgg                                                                                                                      >  1:31804/1‑149 (MQ=255)                              gCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTa                                                                                                                         >  1:30700/1‑144 (MQ=255)                              gCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTa                                                                                                                         <  2:30700/144‑1 (MQ=255)                               ccAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGgtgt                                                                                                                   <  5:185780/149‑1 (MQ=255)                               ccAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTAACCTTTACGgtg                                                                                                                    >  2:124373/1‑148 (MQ=255)                                  acaacCATGAAGACTGACGATTTGATTGCATTGATGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGAAGGCAGTTGCTTCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACg                                                                                                                >  7:230921/1‑149 (MQ=255)                                    aacCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACg                                                                                                                >  8:41917/1‑147 (MQ=255)                                     acCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAg                                                                                                             <  2:157466/149‑1 (MQ=255)                                      ccATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCgcagca                                                                                                                                                <  3:72117/113‑1 (MQ=255)                                      ccATGAAGACAGACGATTAGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGACCCGCGCAACGCCGTACGCAAGCCGAAGGCAGTGGCAGCGCTGGGCCGGCTGACCgcagca                                                                                                                                                >  4:72117/1‑113 (MQ=255)                                       cATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCGGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGt                                                                                                            >  4:91478/1‑148 (MQ=255)                                       cATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCGGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGt                                                                                                            >  4:271847/1‑148 (MQ=255)                                        aTGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGGGAAGGGACAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGAGGCTGCGCTGAGCGGGCTGACCGCAGCACAGCTCCTGACTATTACCCTTTACGGTGTACGCAGTGa                                                                                                          >  2:311186/1‑149 (MQ=255)                                         tGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGat                                                                                                         <  8:136707/149‑1 (MQ=255)                                         tGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGa                                                                                                          >  5:215169/1‑148 (MQ=255)                                           aaGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGata                                                                                                        <  2:300564/148‑1 (MQ=255)                                            aGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGatat                                                                                                       <  8:291053/148‑1 (MQ=255)                                             gACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATca                                                                                                     >  6:2564/1‑149 (MQ=255)                                                tGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATca                                                                                                     >  8:128113/1‑146 (MQ=255)                                                     aTTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACCGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTc                                                                                             <  6:241809/149‑1 (MQ=255)                                                       ttGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCGTTACGGTGTACGCAGTGATATCACAGAGGTCg                                                                                            >  3:305478/1‑148 (MQ=255)                                                         gATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGTGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGca                                                                                          >  7:281601/1‑148 (MQ=255)                                                         gATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGca                                                                                          >  1:4335/1‑148 (MQ=255)                                                         gATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGca                                                                                          <  4:91996/148‑1 (MQ=255)                                                                   cTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCgtt                                                                               <  7:128113/149‑1 (MQ=255)                                                                             gCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCGGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATcaca                                                                                                   <  8:174450/119‑1 (MQ=255)                                                                               gAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCa                                                                    >  8:107791/1‑148 (MQ=255)                                                                                  gggCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGTGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGt                                                                <  8:281601/149‑1 (MQ=255)                                                                                     ccAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGc                                                             <  5:237051/149‑1 (MQ=255)                                                                                         tGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCttt                                                         <  5:51339/149‑1 (MQ=255)                                                                                           ccGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATAACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTc                                                        <  1:71358/148‑1 (MQ=255)                                                                                            cgcgCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCa                                                      <  1:218068/149‑1 (MQ=255)                                                                                             gcgcCACGCCGTTGGCAGGCGGATGGGAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCaa                                                     <  8:230921/149‑1 (MQ=255)                                                                                             gagcCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCaa                                                     >  6:135116/3‑149 (MQ=255)                                                                                              cccttcGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAAc                                                    <  5:272141/144‑1 (MQ=255)                                                                                               gcCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAAc                                                    >  4:188599/1‑148 (MQ=255)                                                                                                ccACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAAc                                                    >  3:117510/1‑147 (MQ=255)                                                                                                ccACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCa                                                  <  4:45864/149‑1 (MQ=255)                                                                                                 cACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCCTTCCAACCAg                                                 >  3:136099/1‑149 (MQ=255)                                                                                                 cACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGAACTGCTGCTGAAAATTACCCTTTACGGTGTACGCAGGGATATCACAGAGGTCGAACAAACCCCGTTGTTATGGGGCAAAGTCGCCTTTCCAACCAg                                                 >  1:33525/1‑149 (MQ=255)                                                                                                          ggCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCGGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGt                                         <  3:91478/148‑1 (MQ=255)                                                                                                          ggCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCGGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGt                                         <  3:271847/148‑1 (MQ=255)                                                                                                           gCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGtt                                        <  4:305478/148‑1 (MQ=255)                                                                                                           gCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGTTg                                       <  4:117510/149‑1 (MQ=255)                                                                                                                 ggATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGTTGATCggg                                 <  7:137059/149‑1 (MQ=255)                                                                                                                 ggATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGTTGATCggg                                 <  6:215169/149‑1 (MQ=255)                                                                                                                       caatGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGTTGATCGGGCTGTgg                           >  5:132188/4‑149 (MQ=255)                                                                                                                              tGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTAGGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGTTGATCGGGCTGTGGTTGAc                      >  6:10786/1‑147 (MQ=255)                                                                                                                              tGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACGAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGTTGATCGGGCTGTGGTTGAc                      >  4:293467/1‑147 (MQ=255)                                                                                                                              tGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGTTGATCGGGCTGTGGTTGAc                      >  1:107166/1‑147 (MQ=255)                                                                                                                                            tGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACAAGCTTGGCGTTGATCGGGCTGTGGTTGACGAGCCGCTTggccagg      >  3:48689/1‑149 (MQ=255)                                                                                                                                                cGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGTTGATCGGGCTGTGGTTGACGAGCCGCTTGgccaggccag  >  6:183712/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TGAAAGCGCTATCCAGGCTGATTGGAGGAGCCAACAACCATGAAGACTGACGATTTGATTGCCTTGCTGGCCGCAGGCGAAGGGCCAGTGCCGCGCCACGCCGTTGGCAGGCGGATGGCAGTGGCTGCGCTGGGCGGGCTGACCGCAGCACTGCTGCTGACTATTACCCTTTACGGTGTACGCAGTGATATCACAGAGGTCGCACAAACCCCGTTGTTCTGGGGCAAAGTCGCCTTTCCAACCAGCTTGGCGTTGATCGGGCTGTGGTTGACGAGCCGCTTGGCCAGGCCAG  >  NC_002947/1133375‑1133666

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