Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,568,701	C→A	7.2%	V283F (GTC→TTC)	PP_2252 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,568,701	0	C	A	7.2%	212.7 / 5.3	69	V283F (GTC→TTC)	PP_2252	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/29);  new base A (3/2);  total (38/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 = 6.82e-01

CGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCTG  >  NC_002947/2568555‑2568846                                                                                                                                                    |                                                                                                                                                   cGGCCGCGGCGCACTGGCCAGCATGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCa                                                                                                                                                  >  2:241613/1‑149 (MQ=255)        ggcgCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCgg                                                                                                                                           <  2:297943/149‑1 (MQ=255)        ggcgCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGAT‑GGTCGAACAACTGCgg                                                                                                                                           <  1:352857/149‑1 (MQ=255)                ggCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTGAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCAc                                                                                                                                   <  7:133889/149‑1 (MQ=255)                 gCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACg                                                                                                                                  >  5:342508/1‑149 (MQ=255)                  ccAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACgg                                                                                                                                 >  2:48923/1‑149 (MQ=255)                         gCTGTCCTCAACCTGCCCCGACATGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCTTCATCTCGTGCCATTCCTGGTTCTACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATggg                                                                                                                          <  1:291697/149‑1 (MQ=255)                         gCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATggg                                                                                                                          >  4:301009/1‑149 (MQ=255)                           tGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGt                                                                                                                         >  4:143491/1‑148 (MQ=255)                             tCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGtcat                                                                                                                      >  3:323817/1‑149 (MQ=255)                             tCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGtcat                                                                                                                      <  6:342508/149‑1 (MQ=255)                             tCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCCGGTGGTCACGGTCATGGGtcat                                                                                                                      >  4:160087/1‑149 (MQ=255)                               cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATgg                                                                                                                           >  8:28743/1‑142 (MQ=255)                               cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATgg                                                                                                                           <  7:28743/142‑1 (MQ=255)                               cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGtcatca                                                                                                                    >  3:25448/1‑149 (MQ=255)                                   aCCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAcccc                                                                                                                <  6:106467/149‑1 (MQ=255)                                   aCCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAcccc                                                                                                                <  3:301009/149‑1 (MQ=255)                                    ccTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCg                                                                                                               >  3:145476/1‑149 (MQ=255)                                      tGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGt                                                                                                              >  4:49526/1‑148 (MQ=255)                                         cccGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCAGGGGTCATCACCCCGTAAcc                                                                                                          >  2:232552/1‑149 (MQ=255)                                             aCAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCt                                                                                                     <  4:25448/149‑1 (MQ=255)                                              cAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCACCGCCCCGTGCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTg                                                                                                    >  5:366763/1‑149 (MQ=255)                                              cAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTg                                                                                                    <  2:292092/149‑1 (MQ=255)                                                ggTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGGTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCg                                                                                                  <  6:373049/149‑1 (MQ=255)                                                   aaTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGc                                                                                                >  1:309084/1‑148 (MQ=255)                                                   aaTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGc                                                                                                >  5:139080/1‑148 (MQ=255)                                                    aTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGCCCATCTGCGGGTGGTCACGGTCATGGGTCATCCCCCCGTACCCG‑GCTGCGGCa                                                                                               >  4:44940/1‑148 (MQ=255)                                                         cccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGa                                                                                          >  2:100861/1‑148 (MQ=255)                                                          ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGaaa                                                                                        <  1:48923/149‑1 (MQ=255)                                                           cccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAat                                                                                       <  1:241613/149‑1 (MQ=255)                                                            ccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAat                                                                                       >  7:168816/1‑148 (MQ=255)                                                                   ggCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGccc                                                                               >  1:72191/1‑149 (MQ=255)                                                                       ccGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGc                                                                           >  2:22409/1‑149 (MQ=255)                                                                        cGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGcc                                                                          >  4:349371/1‑149 (MQ=255)                                                                               gCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGc                                                                   >  5:51828/1‑149 (MQ=255)                                                                                      tcCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGc                                                            >  5:270223/1‑149 (MQ=255)                                                                                        cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTg                                                             <  3:90810/146‑1 (MQ=255)                                                                                        cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGc                                                            <  6:51828/147‑1 (MQ=255)                                                                                         aTCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGa                                                         >  1:173973/1‑149 (MQ=255)                                                                                            aaCTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATg                                                       <  2:309084/148‑1 (MQ=255)                                                                                             aCTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGggtg                                                                                                                                        <  6:313423/67‑1 (MQ=255)                                                                                             aCTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGggtg                                                                                                                                        >  5:313423/1‑67 (MQ=255)                                                                                                    ccATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGTGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACgg                                                 <  1:100861/147‑1 (MQ=255)                                                                                                       ttCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGTGGCATCGGCAGTgg                                           <  4:323817/149‑1 (MQ=255)                                                                                                             gTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATcc                                     >  5:365802/1‑149 (MQ=255)                                                                                                                cAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCcgc                                  >  6:298133/1‑149 (MQ=255)                                                                                                                   cGCCACGTTCAGGATCACCCGCAGTTGGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGt                                >  1:192870/1‑149 (MQ=255)                                                                                                                    gCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTgg                              <  4:145476/149‑1 (MQ=255)                                                                                                                    gCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTgg                              <  3:160087/149‑1 (MQ=255)                                                                                                                     ccACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTggg                             >  4:29154/1‑149 (MQ=255)                                                                                                                      cccGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTgggg                            <  6:366763/147‑1 (MQ=255)                                                                                                                         gTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGc                           >  5:274200/1‑147 (MQ=255)                                                                                                                               gATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGccc                                                                               <  3:96694/89‑1 (MQ=255)                                                                                                                               gATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGccc                                                                               >  4:96694/1‑89 (MQ=255)                                                                                                                               gATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTgg                   >  4:153006/1‑149 (MQ=255)                                                                                                                                aTCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGCTCCCGCGTGGGGCCAGGGTGgc                  >  2:115359/1‑149 (MQ=255)                                                                                                                                 tCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcg                 >  4:335469/1‑149 (MQ=255)                                                                                                                                   aCCCGCAGAT‑GGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCAGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcgc                >  2:108416/1‑148 (MQ=255)                                                                                                                                    cccGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGTGATCCCGCGTGGGGCCAGGGTGGCGCCa              <  6:139080/149‑1 (MQ=255)                                                                                                                                    cccGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGcc               >  5:309058/1‑148 (MQ=255)                                                                                                                                     ccGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCaa             <  2:170096/149‑1 (MQ=255)                                                                                                                                     ccGCAGAT‑GGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTg                                                                            <  7:40935/86‑1 (MQ=38)                                                                                                                                     ccGCAGAT‑GGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTg                                                                            >  8:40935/1‑86 (MQ=38)                                                                                                                                     ccGCAGAT‑GGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGcc               >  8:33222/1‑147 (MQ=255)                                                                                                                                       gCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGc           <  3:49526/149‑1 (MQ=255)                                                                                                                                           tt‑gttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCtttt       <  1:22409/149‑1 (MQ=255)                                                                                                                                              gttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTc      <  2:192870/148‑1 (MQ=255)                                                                                                                                              gttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTc      <  6:309058/148‑1 (MQ=255)                                                                                                                                                  gACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCTg  <  3:335469/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCTG  >  NC_002947/2568555‑2568846

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