Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,212,999	T→C	5.5%	T242T (ACA→ACG)	PP_3695 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,212,999	0	T	C	5.5%	222.4 / 2.3	73	T242T (ACA→ACG)	PP_3695	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base C (2/2);  total (36/37)
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.25e-01

GACATCCCCCCAGTAAGCAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGTATCGTCGTACGCATCCAGAAAGTCGCAGAGCCCTTGC  >  NC_002947/4212852‑4213140                                                                                                                                                     |                                                                                                                                              gACATCCCCCCAGTAAGCAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCtt                                                                                                                                               <  8:268603/148‑1 (MQ=255)     tCCCCCCAGTAAGCAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTccc                                                                                                                                          <  5:38028/149‑1 (MQ=255)     tCCCCCCAGTAAGCAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTccc                                                                                                                                          <  3:200144/149‑1 (MQ=255)     tCCCCCCAGTAAGCAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTcc                                                                                                                                           >  8:137191/1‑148 (MQ=255)      ccccccAGTAAGCAATGGACTTGCCATGGAGTTGGGCGTATCCGAACCAAGCGATATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGAGCAAGATATGAAGCCACTTTTGA‑TTCTCGACTATCAGAAGCCTTCTcc                                                                                                                                           >  6:35988/1‑147 (MQ=255)      ccccccAGTAAGCAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTcc                                                                                                                                           >  2:281454/1‑147 (MQ=255)       cccccAGTAAGCAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAg                                                                                                                                        <  8:304089/149‑1 (MQ=255)        ccccAGTAAGCAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTc                                                                                                                                            <  6:114520/144‑1 (MQ=255)        ccccAGTAAGCAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTc                                                                                                                                            >  5:114520/1‑144 (MQ=255)             gTAAGCAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGc                                                                                                                                  <  4:187624/149‑1 (MQ=255)                  cAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGcc                                                                                                                             <  2:169620/149‑1 (MQ=255)                   aaTGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGccc                                                                                                                            >  2:340529/1‑149 (MQ=255)                   aaTGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGccc                                                                                                                            >  4:328471/1‑149 (MQ=255)                        aCTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAg                                                                                                                                        >  1:319501/1‑132 (MQ=255)                        aCTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAg                                                                                                                                        <  2:319501/132‑1 (MQ=255)                         cTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCaa                                                                                                                      <  8:46347/149‑1 (MQ=255)                         cTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCaa                                                                                                                      <  7:137191/149‑1 (MQ=255)                          ttGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAg                                                                                                                     >  5:160761/1‑149 (MQ=255)                             ccATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCt                                                                                                                   >  4:342210/1‑148 (MQ=255)                              cATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCt                                                                                                                   >  8:349985/1‑147 (MQ=255)                                tGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCa                                                                                                                 >  5:52292/1‑147 (MQ=255)                                    gTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTc                                                                                                           <  6:128908/149‑1 (MQ=255)                                       gggCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGt                                                                                                         <  7:329351/148‑1 (MQ=255)                                        ggCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTg                                                                                                        >  7:32365/1‑148 (MQ=255)                                         gCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGACATCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGa                                                                                                       <  1:292106/149‑1 (MQ=255)                                          cGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGcc                                                                                                                        >  5:349269/1‑130 (MQ=255)                                          cGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGcc                                                                                                                        <  6:349269/130‑1 (MQ=255)                                                 gAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGAc                                                                                              >  3:294522/1‑149 (MQ=255)                                                     cAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGttt                                                                                          <  6:160761/149‑1 (MQ=255)                                                      aaGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGttt                                                                                          >  2:138116/1‑148 (MQ=255)                                                      aaGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGttt                                                                                          >  3:190223/1‑148 (MQ=255)                                                      aaGAGAGATCCGAACCACATCCTAAAATGTCTATCGTTGCGGGCTGCTCGCGCTAGAGATGAAGGCAATTTTTA‑TTCTCGACGATAAGAAGAGTTGTCCCAGTTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTa                                                                                         <  6:123232/149‑1 (MQ=255)                                                             aTCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGgca                                                                                  <  6:52292/149‑1 (MQ=255)                                                              tCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGgcag                                                                                 >  3:231464/1‑149 (MQ=255)                                                              tCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGgca                                                                                  <  7:349985/148‑1 (MQ=255)                                                              tCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGACCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGgcag                                                                                 >  8:264907/1‑149 (MQ=255)                                                                  aCCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGa                                                                             >  3:300165/1‑149 (MQ=255)                                                                   cccACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGAc                                                                            <  4:190223/149‑1 (MQ=255)                                                                    acacaTGATAAAATCGCGATCGTTCCGTGCAGCTCGCGCTAGATATGAAGGAACTTTTTA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGAcc                                                                           <  5:35988/148‑1 (MQ=255)                                                                      acaTCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAg                                                                         >  5:231195/1‑149 (MQ=255)                                                                        aTCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTCGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTaa                                                                                                   >  7:15094/1‑121 (MQ=255)                                                                        aTCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTCGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTaa                                                                                                   <  8:9535/121‑1 (MQ=255)                                                                        aTCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTCGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTaa                                                                                                   >  7:9535/1‑121 (MQ=255)                                                                        aTCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTCGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTaa                                                                                                   <  8:15094/121‑1 (MQ=255)                                                                           cTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTcc                                                                    <  3:342210/149‑1 (MQ=255)                                                                           acgatctCGCGATCGTTCCGGTCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTc                                                                     <  5:234/142‑1 (MQ=255)                                                                               aaTCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCGTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCCGCAGACCAGCCTCCAgcg                                                                >  8:254631/1‑149 (MQ=255)                                                                                   gcgATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCATACCAGccgc                                                                     >  6:234/1‑138 (MQ=255)                                                                                    cgATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGgcgc                                                            >  6:228193/1‑148 (MQ=255)                                                                                     gATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCTGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGACTCCAg                                                                  >  5:101768/1‑141 (MQ=255)                                                                                     gATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAg                                                                  <  6:101768/141‑1 (MQ=255)                                                                                      aTCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGCCCAGCCTCCAGCGgcgc                                                            >  8:226351/1‑146 (MQ=255)                                                                                         gTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGaa                                                       <  4:294522/148‑1 (MQ=255)                                                                                          ttCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAg                                                      >  2:192865/1‑148 (MQ=255)                                                                                          ttCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGc                                                     <  7:264907/149‑1 (MQ=255)                                                                                             cGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAATCAc                                                   <  4:300165/148‑1 (MQ=255)                                                                                              gggCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCAcc                                                  >  7:355270/1‑148 (MQ=255)                                                                                               ggCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCa                                                <  1:138116/149‑1 (MQ=255)                                                                                                gCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAg                                               <  1:281454/149‑1 (MQ=255)                                                                                                             gagaTGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGaa                                                       <  5:233856/128‑1 (MQ=255)                                                                                                             gagaTGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGaa                                                       >  6:233856/1‑128 (MQ=255)                                                                                                              agaTGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGTAtcgtc                                 <  4:231464/149‑1 (MQ=255)                                                                                                                      gCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGTATCGTCGTACGCAt                         >  3:335776/1‑149 (MQ=255)                                                                                                                      gCACTATTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGTATCGTCGTACGCAt                         >  2:247746/1‑149 (MQ=255)                                                                                                                                    tcGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTg                                                                                                                          <  6:301478/39‑1 (MQ=255)                                                                                                                                    tcGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTg                                                                                                                          >  5:301478/1‑39 (MQ=255)                                                                                                                                        cgaTCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGTATCGTCGTACGCATCCAGAAAGTCGCAGAGc        <  6:231195/149‑1 (MQ=255)                                                                                                                                        cgaTCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGTATCGTCGTACGCATCCAGAAAGTCGCAGAGc        <  1:192865/149‑1 (MQ=255)                                                                                                                                           tCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGTTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGt                                       <  2:343099/115‑1 (MQ=255)                                                                                                                                           tCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGTTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGt                                       >  1:343099/1‑115 (MQ=255)                                                                                                                                            cAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGTATCGTCGTACGCATCCAGAAAGTCGCAGAGCcctt    >  1:68661/1‑149 (MQ=255)                                                                                                                                             aGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGTATCGTCGTACGCATCCAGAAAGTCGCAGAGCcctt    <  5:228193/148‑1 (MQ=255)                                                                                                                                             aGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGTATCGTCGTACGCATCCAGAAAGTCGCAGAGCcctt    <  8:2839/148‑1 (MQ=255)                                                                                                                                              gAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGTATCGTCGTACGCATCCAGAAAGTCGCAGAGCccttgc  >  5:70878/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                              GACATCCCCCCAGTAAGCAATGGACTTGCCATGGAGTTGGGCGGATCCGAGCCAAGCGAGATCCGACCCACATCCTAAAATCGCGATCGTTCCGGGCAGCTCGCGCAAGAGATGAAGGCACTTTTGA‑TTCTCGACGATCAGAAGCCTTGTCCCAGGTAAGCTTGCCCTGCCAAGCTCAGAGGTCGTGACTAACTGACGTTTAAACGGCAGCAGACCAGCCTCCAGCGGCGCCAGAAGCACCCAGTGACTCGTATCGTCGTACGCATCCAGAAAGTCGCAGAGCCCTTGC  >  NC_002947/4212852‑4213140

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