Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,568,693	T→A	8.6%	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	8.6%	201.5 / 6.1	69	Q285H (CAA→CAT)	PP_2252	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/23);  new base A (3/3);  total (43/26)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CTCCACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGC  >  NC_002947/2568549‑2568837                                                                                                                                                  |                                                                                                                                                 ctCCACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑Agttgtt                                                                                                                                              <  2:223492/149‑1 (MQ=255)    cACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGa                                                                                                                                           >  7:2433/1‑149 (MQ=255)    cACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGa                                                                                                                                           >  2:190468/1‑149 (MQ=255)              ggcgCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGATGGTCGAACAACTGCgg                                                                                                                                 <  2:338305/149‑1 (MQ=255)               gcgcACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCggg                                                                                                                                <  8:270670/149‑1 (MQ=255)                  cACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGggtgg                                                                                                                             <  6:286324/149‑1 (MQ=255)                   aCTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGggtggt                                                                                                                            >  7:45571/1‑149 (MQ=255)                      ggCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCAc                                                                                                                         >  8:253762/1‑149 (MQ=255)                        ccAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACgg                                                                                                                       <  8:2433/149‑1 (MQ=255)                         cAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGt                                                                                                                      >  2:263250/1‑149 (MQ=255)                                   tCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtca                                                                                                             >  3:156598/1‑148 (MQ=255)                                   tCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtc                                                                                                              <  1:190468/147‑1 (MQ=255)                                     cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtcatc                                                                                                           <  8:45571/148‑1 (MQ=255)                                      tCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAc                                                                                                         >  7:105816/1‑149 (MQ=255)                                      tCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAc                                                                                                         >  3:112050/1‑149 (MQ=255)                                        aaCCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAcc                                                                                                        <  1:263250/148‑1 (MQ=255)                                          ccTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCg                                                                                                     <  3:150365/149‑1 (MQ=255)                                           cTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCg                                                                                                     >  5:304940/1‑148 (MQ=255)                                               cccGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAAcc                                                                                                >  6:291128/1‑149 (MQ=255)                                                       gTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCgg                                                                                        <  6:304940/149‑1 (MQ=255)                                                         aaTAGACCCCCGCCATGGCCCCGCTCTGGCGGGTCTCCATCAACTCCTGCCATTCCTGGTTCACCGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGGCATCACCCCGTAACCGGCTGCGGc                                                                                       >  4:31673/1‑148 (MQ=255)                                                          aTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGc                                                                                       <  5:149481/147‑1 (MQ=255)                                                          aTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGc                                                                                       >  6:149481/1‑147 (MQ=255)                                                            agccccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCCCCCCGTAACCGGCTGCGGCAATc                                                                                   >  5:301961/4‑149 (MQ=255)                                                            aGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATc                                                                                   >  2:166556/1‑149 (MQ=255)                                                              aCCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCg                                                                                  >  3:34306/1‑148 (MQ=255)                                                               cccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGa                                                                                 >  7:76536/1‑148 (MQ=255)                                                                ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACg                                                                                                                        <  8:199928/108‑1 (MQ=255)                                                                ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACg                                                                                                                        >  7:194697/1‑108 (MQ=255)                                                                ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACg                                                                                                                        <  8:194697/108‑1 (MQ=255)                                                                ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACg                                                                                                                        >  7:199928/1‑108 (MQ=255)                                                                ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGaaa                                                                               <  7:253762/149‑1 (MQ=255)                                                                ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGaaa                                                                               >  8:214345/1‑149 (MQ=255)                                                                       aTGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGc                                                                        >  2:245184/1‑149 (MQ=255)                                                                            cccGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTg                                                                   >  4:136358/1‑149 (MQ=255)                                                                                     gCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCtt                                                            <  7:214345/147‑1 (MQ=255)                                                                                      cGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGcc                                                         <  8:76536/149‑1 (MQ=255)                                                                                           ctcCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTg                                                    >  2:187864/1‑149 (MQ=255)                                                                                              cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGcc                                                  >  1:266078/1‑149 (MQ=255)                                                                                              cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGcc                                                  >  5:301811/1‑148 (MQ=255)                                                                                              cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGcc                                                  >  7:178100/1‑148 (MQ=255)                                                                                                     tCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCAc                                          >  3:149152/1‑149 (MQ=255)                                                                                                       cTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCAc                                          <  4:156598/147‑1 (MQ=255)                                                                                                       cTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACgg                                        >  2:82704/1‑149 (MQ=255)                                                                                                        tGCCATTCCTGGTTCAACGACACGTTCAGGATAACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACTGTCATGGGTCATCACCCCGTAAACGGCAGCGGCAATCGAACTCATCGCCCCTGCCACCTTGCACACTGCCGATGCCAAGGc                                       >  4:104497/1‑149 (MQ=255)                                                                                                          ccATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCAc                                     >  3:275051/1‑149 (MQ=255)                                                                                                                cTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATg                                              >  5:122196/1‑134 (MQ=255)                                                                                                                cTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATg                                              <  6:122196/134‑1 (MQ=255)                                                                                                                cTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGa                               >  7:274057/1‑149 (MQ=255)                                                                                                                  ggTTCAACGCCACGTTCAGGATCACCCGC‑AGTAGTTCGACCATCTGCTGGTGGTCACGGTCATGGGTCATAACCCCGTAACCGCCTACGGCAATCGAAGTCATAGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTTGCCATc                             >  7:190341/1‑149 (MQ=255)                                                                                                                      cAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCcgc                         >  4:77637/1‑149 (MQ=255)                                                                                                                       aaCGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCcgc                         <  4:149152/148‑1 (MQ=255)                                                                                                                        aCGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGt                       <  2:26787/149‑1 (MQ=255)                                                                                                                         cGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTg                      <  6:332830/149‑1 (MQ=255)                                                                                                                               gTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGc                  >  1:184170/1‑147 (MQ=255)                                                                                                                                 tCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCTTGGCGCCacc              >  5:111533/1‑147 (MQ=255)                                                                                                                                    tgaccACCCGC‑AGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTg           <  8:249516/145‑1 (MQ=255)                                                                                                                                     gATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCAcc                                                              >  4:247523/1‑97 (MQ=255)                                                                                                                                     gATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCAcc                                                              <  3:247523/97‑1 (MQ=255)                                                                                                                                     gATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTgg          >  8:186319/1‑149 (MQ=255)                                                                                                                                       tCACGCTC‑AGTTGTTCTACCATCTGCGTTTGGTCACGGTCATGGGTCATCAGCCCTTATCCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACATTGCCCACTGCAGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgc         <  8:190341/148‑1 (MQ=255)                                                                                                                                       tCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcg        <  4:112050/149‑1 (MQ=255)                                                                                                                                        cACCCGC‑AGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgc         <  1:90715/147‑1 (MQ=255)                                                                                                                                        cACCCGC‑AGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgc         >  2:90715/1‑147 (MQ=255)                                                                                                                                         aCCCGC‑AGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcgc       >  1:41304/1‑148 (MQ=255)                                                                                                                                          cccGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGTCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGcc      >  6:147197/1‑148 (MQ=255)                                                                                                                                          cccGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGcc      >  2:149325/1‑148 (MQ=255)                                                                                                                                           ccGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGcc      >  1:106070/1‑147 (MQ=255)                                                                                                                                           ccGC‑AGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCCCTGGCGATCCCGCGGGGGGCCAGGGTGGCGcc      >  4:255232/1‑147 (MQ=255)                                                                                                                                             gC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGc  <  8:105816/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CTCCACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGC  >  NC_002947/2568549‑2568837

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