Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,568,704	T→A	7.2%	M282L (ATG→TTG)	PP_2252 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,568,704	0	T	A	7.2%	202.2 / 5.0	69	M282L (ATG→TTG)	PP_2252	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (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 = 8.65e-01

GGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCTGGCCT  >  NC_002947/2568562‑2568850                                                                                                                                                |                                                                                                                                                    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                                                                                                                               <  7:28743/142‑1 (MQ=255)                        cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATgg                                                                                                                               >  8:28743/1‑142 (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                                                                                                    >  5:139080/1‑148 (MQ=255)                                            aaTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGc                                                                                                    >  1:309084/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                                  <  3:160087/149‑1 (MQ=255)                                                                                                             gCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTgg                                  <  4:145476/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                                                                                >  8:40935/1‑86 (MQ=38)                                                                                                                              ccGCAGAT‑GGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTg                                                                                <  7:40935/86‑1 (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          <  6:309058/148‑1 (MQ=255)                                                                                                                                       gttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTc          <  2:192870/148‑1 (MQ=255)                                                                                                                                           gACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCTg      <  3:335469/148‑1 (MQ=255)                                                                                                                                              cATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCTACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCTGGCCt  >  1:129712/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    GGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCTGGCCT  >  NC_002947/2568562‑2568850

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