Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	336,125	+T	coding (210/237 nt)	PP_0278 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	336,124	1	.	T	98.4%	207.2 / ‑3.3	64	coding (209/237 nt)	PP_0278	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (32/31);  total (33/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.81e-01

ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGA‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGA  >  NC_002947/335980‑336263                                                                                                                                                   |                                                                                                                                            aCCCATCGCTCCTTCTCCTACCTGTAGAGCCCTGCGTCCCCGTCGTCGCGGCCTACGAGCATTTCTCGATCCTGCTGGGCTACATCAAGGACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGA‑ttac                                                                                                                                         <  8:957/149‑4 (MQ=255)   ccATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGA‑ttacct                                                                                                                                       >  2:368129/1‑144 (MQ=255)   ccATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGA‑ttacct                                                                                                                                       >  1:287586/1‑144 (MQ=255)         cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGa                                                                                                                                  <  7:337512/149‑1 (MQ=255)         cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGA‑ttacct                                                                                                                                       <  3:169572/143‑6 (MQ=255)         cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGA‑ttacct                                                                                                                                       >  4:169572/1‑138 (MQ=255)             aTCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATg                                                                                                                              >  1:80977/1‑149 (MQ=255)                  cTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGcc                                                                                                                           >  4:171213/1‑147 (MQ=255)                   tACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAg                                                                                                                        >  7:220162/1‑149 (MQ=255)                   tACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAg                                                                                                                        >  2:25939/1‑149 (MQ=255)                     cgcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATg                                                                                                                              >  1:301615/1‑141 (MQ=255)                     cgcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATg                                                                                                                              <  2:301615/141‑1 (MQ=255)                     cgcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCTAgg                                                                                                                       <  2:80977/148‑1 (MQ=255)                       cgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGgcgc                                                                                                                    <  5:383386/149‑1 (MQ=255)                           gagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAgg                                                                                                                       <  7:342010/142‑1 (MQ=255)                           gagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAgg                                                                                                                       >  8:342010/1‑142 (MQ=255)                           gagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCtgat                                                                                                                <  2:358956/149‑1 (MQ=255)                           gagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCtga                                                                                                                 <  2:328438/148‑1 (MQ=255)                             gCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCtgatg                                                                                                               >  3:360523/1‑148 (MQ=255)                              cccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGaa                                                                                                             <  1:368129/149‑1 (MQ=255)                                cGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACg                                                                                                           <  3:162091/149‑1 (MQ=255)                                       cccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGTGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCg                                                                                                     <  5:304801/148‑1 (MQ=255)                                           tCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATc                                                                                                <  2:287586/149‑1 (MQ=255)                                              aCGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGc                                                                                              >  2:258374/1‑148 (MQ=255)                                              aCGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACCATATCGAGATCGc                                                                                              >  7:368037/1‑148 (MQ=255)                                               cgcgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGcc                                                                                             >  4:128358/1‑148 (MQ=255)                                                gcgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCaa                                                                                           <  3:171213/149‑1 (MQ=255)                                                           gCATGTCTCGATCCTGATGTGCTACACCAAGAACCTGGTGTGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGCCAGACCAT‑TACCTGTGCAACATGGAAAAGGCGCTCATGAACTAAATATAGAACTCCAAGTGTAAggag                                                                                 >  2:32647/1‑146 (MQ=255)                                                                  tcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                          >  7:188488/1‑148 (MQ=255)                                                                   cGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCa                                                                        >  4:169552/1‑149 (MQ=255)                                                                    gATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgcg                                                                                 <  3:305883/139‑1 (MQ=255)                                                                    gATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgcg                                                                                 >  4:305883/1‑139 (MQ=255)                                                                    gATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGc                                                                         <  1:258374/147‑1 (MQ=255)                                                                     aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAtt                                                                      <  4:360523/149‑1 (MQ=255)                                                                         tGATGGGCTACATCAAGCACCTGCTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAcg                                                                  >  3:17341/1‑149 (MQ=255)                                                                             gggCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGt                                                              >  2:179744/1‑149 (MQ=255)                                                                                   cATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtg                                                        >  6:52785/1‑149 (MQ=255)                                                                                    aTCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgt                                                       <  5:319615/149‑1 (MQ=255)                                                                                     tCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgtg                                                      >  1:30243/1‑149 (MQ=255)                                                                                     tCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgtg                                                      >  7:336387/1‑149 (MQ=255)                                                                                      cAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                     >  5:15205/1‑149 (MQ=255)                                                                                        aGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCt                                                                                                                   >  6:309221/1‑85 (MQ=255)                                                                                        aGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCt                                                                                                                   <  5:309221/85‑1 (MQ=255)                                                                                        aGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                     >  5:107425/1‑147 (MQ=255)                                                                                            ccTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGTGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCgg                                               <  2:303111/149‑1 (MQ=255)                                                                                              tGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCgggc                                             >  3:378292/1‑147 (MQ=255)                                                                                               ggTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGCCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGccc                                            >  5:46467/1‑148 (MQ=255)                                                                                                gTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCtt                                           <  4:148136/149‑1 (MQ=255)                                                                                                            gATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGCGCGGCCTTGCGTCGCGAATg                               <  8:188488/149‑1 (MQ=255)                                                                                                              tatGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATggg                             <  3:196683/149‑1 (MQ=255)                                                                                                              tatGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATggg                             <  1:179744/149‑1 (MQ=255)                                                                                                                tGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCt                           <  6:15205/149‑1 (MQ=255)                                                                                                                     gacgacCATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCaaa                      <  6:107425/149‑1 (MQ=255)                                                                                                                       cgcacATAAATTCCTTGGGTCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACCATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAagc                    <  8:368037/145‑1 (MQ=255)                                                                                                                          ccATAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGTGCTGCAAagcagc                 >  1:193222/1‑149 (MQ=255)                                                                                                                            aTAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTGGCGAATGGGCGGCAAAGGAGccc               >  5:360668/1‑149 (MQ=255)                                                                                                                             tAAATTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGcccc              >  3:71231/1‑149 (MQ=255)                                                                                                                               aaTTCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCgg            <  8:220162/149‑1 (MQ=255)                                                                                                                                 ttCCTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGc           >  1:268475/1‑148 (MQ=255)                                                                                                                                    cTGGGGGCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTAGCGAATGGGCTGCAAAGCAGCCCCGGCGAtt       >  4:182851/1‑149 (MQ=255)                                                                                                                                       ggggCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCTAATGGGCTGCAAAGCAGCCCCGGCGATTTCa    <  5:52785/149‑1 (MQ=255)                                                                                                                                       ggggCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTc     <  3:169552/148‑1 (MQ=255)                                                                                                                                       ggggCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCa    <  4:71231/149‑1 (MQ=255)                                                                                                                                        gggCAG‑CCGATTACCTGTGCGACATGGGCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAttt      >  4:234356/1‑146 (MQ=255)                                                                                                                                         ggCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCa                        >  6:208634/1‑127 (MQ=255)                                                                                                                                         ggCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCa                        <  5:208634/127‑1 (MQ=255)                                                                                                                                         ggCAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCaga  <  1:32647/149‑1 (MQ=255)                                                                                                                                           cAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCa    >  5:301538/1‑145 (MQ=255)                                                                                                                                           cAG‑CCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCa    >  4:137741/1‑145 (MQ=255)                                                                                                                                                   |                                                                                                                                            ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAG‑CCGA‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGA  >  NC_002947/335980‑336263

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