Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	990,583	A→C	5.6%	D155A (GAC→GCC)	PP_0856 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	990,583	0	A	C	5.6%	225.1 / 3.7	71	D155A (GAC→GCC)	PP_0856	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (39/28);  new base C (2/2);  total (41/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 = 8.16e-01

CCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGT  >  NC_002947/990440‑990730                                                                                                                                                 |                                                                                                                                                    ccTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATc                                                                                                                                                >  2:229805/1‑149 (MQ=255)  cTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCAGACAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCg                                                                                                                                               >  1:127200/1‑149 (MQ=255)                   ccTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATg                                                                                                                              >  3:763/1‑149 (MQ=255)                    cTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTAGTGGTGCAGACCCAGGCCGATCAGGCGATCGTCCTgg                                                                                                                                 >  5:342702/1‑145 (MQ=255)                    cTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTAGTGGTGCAGACCCAGGCCGATCAGGCGATCGTCCTgg                                                                                                                                 <  6:342702/145‑1 (MQ=255)                       aaGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATgcc                                                                                                                            <  5:7403/147‑1 (MQ=255)                         gggTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCa                                                                                                                        <  5:16167/149‑1 (MQ=255)                          ggTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCt                                                                                                                                   >  8:165528/1‑137 (MQ=255)                          ggTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCt                                                                                                                                   <  7:165528/137‑1 (MQ=255)                            tGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTg                                                                                                                     >  2:335483/1‑149 (MQ=255)                                          tGGACTCCAGCACCGGTGAGGAGCGTTGGCGGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTg                                                                                                        >  1:226146/1‑149 (MQ=255)                                          tGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCAACCGCCGTTg                                                                                                        >  1:344132/1‑148 (MQ=255)                                            gACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTgg                                                                                                       <  1:358031/147‑1 (MQ=255)                                            gACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGa                                                                                                      <  2:226146/148‑1 (MQ=255)                                            gACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCAACCGCCGTTg                                                                                                        <  2:344132/146‑1 (MQ=255)                                              cTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCAGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATc                                                                                                    >  5:113975/1‑148 (MQ=255)                                               tCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTa                                                                                                  >  2:266113/1‑149 (MQ=255)                                                 cAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGACAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTAc                                                                                                 <  2:127200/148‑1 (MQ=255)                                                      ccGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAAc                                                                                           <  2:79132/149‑1 (MQ=255)                                                       cGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACt                                                                                          >  2:191800/1‑149 (MQ=255)                                                              gagCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATc                                                                                                                                                <  1:369502/88‑1 (MQ=255)                                                              gagCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATc                                                                                                                                                >  2:369502/1‑88 (MQ=255)                                                              gagCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgccggc                                                                                   >  7:262190/1‑149 (MQ=255)                                                              gagCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgccggc                                                                                   >  5:356830/1‑149 (MQ=255)                                                               agCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgccggcc                                                                                  >  5:151394/1‑149 (MQ=255)                                                                 cGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGt                                                                                >  2:153026/1‑149 (MQ=255)                                                                      gc‑gcTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgc                                                                                       >  1:352788/1‑137 (MQ=255)                                                                      gc‑gcTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgc                                                                                       <  2:352788/137‑1 (MQ=255)                                                                         gcTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACg                                                                         <  1:153026/149‑1 (MQ=255)                                                                         gcTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACg                                                                         >  2:273696/1‑149 (MQ=255)                                                                           tCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGTGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCt                                                                       >  6:84210/1‑149 (MQ=255)                                                                               gcgTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGt                                                                   >  3:213684/1‑149 (MQ=255)                                                                               gcgTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGt                                                                   <  1:335483/149‑1 (MQ=255)                                                                                    aCCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGCCGATCAGGCGATCGTCCTgg                                                                                                                                 <  5:169186/82‑1 (MQ=37)                                                                                    aCCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGCCGATCAGGCGATCGTCCTgg                                                                                                                                 >  6:169186/1‑82 (MQ=37)                                                                                     ccAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgccgg                                                                                    <  1:205275/126‑1 (MQ=255)                                                                                     ccAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgccgg                                                                                    >  2:205275/1‑126 (MQ=255)                                                                                      cAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCg                                                            <  4:213684/149‑1 (MQ=255)                                                                                       aGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCgg                                                           >  6:105018/1‑149 (MQ=255)                                                                                            aGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGa                                                                           <  2:217069/128‑1 (MQ=255)                                                                                            aGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGa                                                                           >  1:217069/1‑128 (MQ=255)                                                                                                   gCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCg                                                    >  6:217834/1‑144 (MQ=255)                                                                                                   gCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCg                                                    <  5:217834/144‑1 (MQ=255)                                                                                                   gCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGc                                               <  5:105018/149‑1 (MQ=255)                                                                                                    ccccACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCCCCGGTGCGCCGATCGc                                               >  3:70492/1‑148 (MQ=255)                                                                                                      ccACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGccac                                            >  6:294618/1‑149 (MQ=255)                                                                                                         ccTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCaa                                         <  4:763/149‑1 (MQ=255)                                                                                                               aacaacGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCt                                   >  2:369519/1‑149 (MQ=255)                                                                                                                acaacGGTGACGTGGTGGAGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTg                                  >  7:20277/1‑149 (MQ=255)                                                                                                                   acGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCg                               <  6:356830/149‑1 (MQ=255)                                                                                                                     ggTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGt                             >  2:108952/1‑149 (MQ=255)                                                                                                                       tGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTg                           >  3:23820/1‑149 (MQ=255)                                                                                                                         aCGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGc                          <  8:262190/148‑1 (MQ=255)                                                                                                                              gtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCTCCGTTGGATCTACGAATACTCGCCGGCCGTGCTGACGGTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCt                    <  6:113975/149‑1 (MQ=255)                                                                                                                              gtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCt                    >  8:85392/1‑149 (MQ=255)                                                                                                                              gtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGCCACCGGTGCACCGATCGCCACCAACCGCCTGGCGGTTACCGGCCt                    <  1:266113/149‑1 (MQ=255)                                                                                                                                  tggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTc                 >  3:266481/1‑148 (MQ=255)                                                                                                                                    gtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTAc              >  8:293897/1‑149 (MQ=255)                                                                                                                                      gCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCAGGCGGTTggtggc                      >  4:134369/1‑134 (MQ=255)                                                                                                                                      gCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCAGGCGGTTggtggc                      <  3:134369/139‑6 (MQ=255)                                                                                                                                      gCAGAAACAGGACGATCGCCTGATCGGCCTGGATGCAGCCACTGGCGAACGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCAACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCg            >  8:304815/1‑149 (MQ=255)                                                                                                                                          aCCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGCCCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGTGTGGCACCGGTgc                                                        >  4:29163/1‑101 (MQ=255)                                                                                                                                          aCCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGTGTGGCACCGGTgc                                                        <  3:29163/101‑1 (MQ=255)                                                                                                                                            ccAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAgg      >  3:44274/1‑149 (MQ=255)                                                                                                                                            ccAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAgg      >  8:155465/1‑149 (MQ=255)                                                                                                                                              aGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCgg           >  4:82941/1‑142 (MQ=255)                                                                                                                                              aGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCgg           <  3:82941/142‑1 (MQ=255)                                                                                                                                              aGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTa    <  1:369519/149‑1 (MQ=255)                                                                                                                                                gACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCgg                       >  2:296965/1‑128 (MQ=255)                                                                                                                                                gACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCgg                       <  1:296965/128‑1 (MQ=255)                                                                                                                                                gACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGt  <  1:229805/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGC‑GCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGT  >  NC_002947/990440‑990730

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