Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,568,693	T→A	7.3%	Q285H (CAA→CAT)	PP_2252 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,568,693	0	T	A	7.3%	196.5 / 5.0	68	Q285H (CAA→CAT)	PP_2252	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/26);  new base A (3/2);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

GCTGCTCCACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTT  >  NC_002947/2568545‑2568841                                                                                                                                                     |                                                                                                                                                       gctgctCCACCGGCCGCGGCGCACTCGCCAGCAAGCTGTCCTCAACCTCCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAgt                                                                                                                                                        >  4:90810/1‑149 (MQ=255)     ctCCAAAGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGAAATAGACCCCCGCCATGGCCCCCCTCTGGCGGTTCTCCATCAACTACTGACATTCCTGGTTCAACGCCACCATCAGGATCACCAGCAgtacat                                                                                                                                                    >  2:291697/1‑145 (MQ=255)        cgccGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGa                                                                                                                                                <  3:93880/147‑1 (MQ=255)        cACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGGTCAGGATCACCCGCACTT‑GTTCGa                                                                                                                                                >  8:133889/1‑149 (MQ=255)        cACCGGCAGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGa                                                                                                                                                >  5:106467/1‑149 (MQ=255)           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                                                                                           >  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                         <  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)                                                                                                                                                     |                                                                                                                                                       GCTGCTCCACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTT  >  NC_002947/2568545‑2568841

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