Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,568,696	T→G	7.4%	E284D (GAA→GAC)	PP_2252 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,568,696	0	T	G	7.4%	210.5 / 7.3	68	E284D (GAA→GAC)	PP_2252	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/26);  new base G (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 = 9.19e-01

CACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTC  >  NC_002947/2568552‑2568842                                                                                                                                                  |                                                                                                                                                    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                                                                                            >  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                                                                                                                                    >  5:313423/1‑67 (MQ=255)                                                                                                aCTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGggtg                                                                                                                                    <  6:313423/67‑1 (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                                                                           >  4:96694/1‑89 (MQ=255)                                                                                                                                  gATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGccc                                                                           <  3:96694/89‑1 (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  <  6:309058/148‑1 (MQ=255)                                                                                                                                                 gttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTc  <  2:192870/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTT‑GTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCG‑GCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTC  >  NC_002947/2568552‑2568842

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