Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,116,922	A→C	8.3%	L158R (CTC→CGC)	PP_2733 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,116,922	0	A	C	8.3%	210.1 / 10.4	72	L158R (CTC→CGC)	PP_2733	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/36);  new base C (4/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

GGGGCAGGCGTTGTGAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCCAGCAGTGCGCCGAT  >  NC_002947/3116774‑3117070                                                                                                                                                     |                                                                                                                                                      ggggCAGGCGTTGTGAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGa                                                                                                                                                       <  1:62896/149‑1 (MQ=255)        gCGTTGGGACACCGAGCCTTCTGCCCAGCACGGGTTGCGTGGGGTGGGGCCGGCCGGTCCCCTGGCGCAAGCGCGGGGCACGCTGGTGGGCGCAGGGGCGTTCCTGCTCCGAGCGAGGTCGCCGGTCAGCGCCGCGCCGCGAGgcgcgc                                                                                                                                                <  7:326111/149‑1 (MQ=255)           ttGTGAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCTCACGgcgcggcgc                                                                                                                                             <  5:386695/149‑5 (MQ=255)              tgAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGGGGTGTGTGGCCGGCCGTGCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGTTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCcagcag                                                                                                                                          <  5:272556/149‑1 (MQ=255)                aaaaaGAGCCTTCTGCCCAGCAAGGTTTGCGTGGGGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAgcagg                                                                                                                                         <  8:104830/148‑1 (MQ=255)                  caaGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTGGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAgcaggggc                                                                                                                                      <  1:216925/148‑1 (MQ=255)                     gagCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAgcaggggca‑gg                                                                                                                                  <  1:59861/149‑1 (MQ=255)                         cTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGcc                                                                                                                               >  7:263389/1‑148 (MQ=255)                          gccTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCaa                                                                                                                             <  3:256657/147‑1 (MQ=255)                              gCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGAGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGcc                                                                                                                         >  5:107055/1‑149 (MQ=255)                               cccAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGccc                                                                                                                        <  2:336537/149‑1 (MQ=255)                               cccAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGcc                                                                                                                         >  1:153861/1‑148 (MQ=255)                                  aGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCAt                                                                                                                      >  5:94672/1‑148 (MQ=255)                                  aGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGACA‑GGGCCAAAGCCCATa                                                                                                                     >  3:266844/1‑149 (MQ=255)                                                 gtgtgtGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAg                                                                                                        >  1:70089/1‑147 (MQ=255)                                                      tgGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCg                                                                                                 <  7:144302/149‑1 (MQ=255)                                                         ccggccgTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGtt                                                                                              <  2:70089/149‑1 (MQ=255)                                                             ccgTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCa                                                                                          <  2:158991/149‑1 (MQ=255)                                                                 tCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAgcgc                                                                                      >  7:18518/1‑149 (MQ=255)                                                                  cACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAgcgc                                                                                      <  8:172582/148‑1 (MQ=255)                                                                     cTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGcc                                                                                   >  6:212637/1‑148 (MQ=255)                                                                     cTGGCGCAAGCGCTGGGCACGCTATTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAgcgcgc                                                                                    <  5:50866/147‑1 (MQ=255)                                                                     cTGGCGCAAGCGCTGGGCACGCTATTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAgcgcgc                                                                                    >  6:50866/1‑147 (MQ=255)                                                                     cTGGAGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGcc                                                                                   >  1:267028/1‑148 (MQ=255)                                                                        gcgcAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCa                                                                                                                       >  1:369909/1‑109 (MQ=255)                                                                        gcgcAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCa                                                                                                                       <  2:369909/109‑1 (MQ=255)                                                                        gcgcAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCaa                                                                                <  6:94672/148‑1 (MQ=255)                                                                           cAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATggc                                                                            <  8:314403/149‑1 (MQ=255)                                                                                     gCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAacac                                                                  <  1:318468/149‑1 (MQ=255)                                                                                      cACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACAcc                                                                 >  3:56654/1‑149 (MQ=255)                                                                                         gCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGc                                                              >  2:59826/1‑149 (MQ=255)                                                                                           tGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCa                                                            <  6:107055/149‑1 (MQ=255)                                                                                             ttGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGcc                                                                                                                               >  8:317773/1‑80 (MQ=255)                                                                                             ttGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGcc                                                                                                                               <  7:317773/80‑1 (MQ=255)                                                                                             ttGGGAGCAGGGGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGACA‑GCGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCATCGCGCCCAATGGCAAGGCAACAGCTGCCAGc                                                          <  4:266844/149‑1 (MQ=255)                                                                                                 gAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTgg                                                      >  3:8774/1‑149 (MQ=255)                                                                                                 gAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTgg                                                      >  2:279644/1‑149 (MQ=255)                                                                                                 gAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGAGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTgg                                                      >  5:239533/1‑149 (MQ=255)                                                                                                  aGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGgc                                                     <  2:267028/149‑1 (MQ=255)                                                                                                     aGGTGCGTTACTGCGCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATACCGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGcgcc                                                  <  5:212637/149‑1 (MQ=255)                                                                                                       gTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATggca                                                                           >  3:219909/1‑122 (MQ=255)                                                                                                       gTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATggca                                                                           <  4:219909/122‑1 (MQ=255)                                                                                                        tGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGAGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCgccgcc                                               <  6:239533/149‑1 (MQ=255)                                                                                                         gCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATggcaaggcaa                                                                     <  5:77496/126‑1 (MQ=255)                                                                                                         gCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATggcaaggcaa                                                                     >  6:77496/1‑126 (MQ=255)                                                                                                           gTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGTGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTAc                                            <  5:77797/149‑1 (MQ=255)                                                                                                            ttACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACa                                           <  3:316305/149‑1 (MQ=255)                                                                                                                  cTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAac                                                                    <  4:144300/118‑1 (MQ=255)                                                                                                                  cTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAac                                                                    >  3:144300/1‑118 (MQ=255)                                                                                                                   tCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGc                                      >  3:136599/1‑147 (MQ=255)                                                                                                                   tCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGc                                      >  5:21323/1‑147 (MQ=255)                                                                                                                   tCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGAGCCGCCTACATAGGc                                      >  3:288607/1‑147 (MQ=255)                                                                                                                     cGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCaga                                  <  8:263389/149‑1 (MQ=255)                                                                                                                     cGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCa                                    >  8:51908/1‑147 (MQ=255)                                                                                                                      gAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCagag                                 <  1:59826/149‑1 (MQ=255)                                                                                                                      gAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCagag                                 <  6:167618/149‑1 (MQ=255)                                                                                                                        gCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCagag                                 >  8:222601/1‑147 (MQ=255)                                                                                                                        gCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCagag                                 <  7:222601/147‑1 (MQ=255)                                                                                                                        gCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAggg                               >  8:193109/1‑149 (MQ=255)                                                                                                                           aTGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGgccg                            >  7:219752/1‑149 (MQ=255)                                                                                                                            tGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGgccgc                           <  8:18518/149‑1 (MQ=255)                                                                                                                                gCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAAt                       >  7:42116/1‑149 (MQ=255)                                                                                                                                 ccGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATc                      <  6:21323/149‑1 (MQ=255)                                                                                                                                 ccGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATc                      <  5:213823/149‑1 (MQ=255)                                                                                                                                    gTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCAGGGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGc                    >  1:69742/1‑149 (MQ=255)                                                                                                                                    gTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGcc                   >  8:205647/1‑149 (MQ=255)                                                                                                                                        ctggcgcgcctcgCGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCca               >  4:47780/5‑149 (MQ=255)                                                                                                                                         tggcgcgcctcgCGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCcag              >  3:86107/4‑149 (MQ=255)                                                                                                                                         tggcgcgcctcgCGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCcag              >  6:322284/4‑149 (MQ=255)                                                                                                                                          ggcgcgcctcgCGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAacac                                                                  <  5:338706/94‑1 (MQ=255)                                                                                                                                          ggcgcgcctcgCGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAacac                                                                  >  6:338706/3‑96 (MQ=255)                                                                                                                                                  gcgAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCCAGCAGTGCGcc     <  7:51908/149‑1 (MQ=255)                                                                                                                                                    gAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCTGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCCAGCAGTGCGCCGa   <  1:279644/149‑1 (MQ=255)                                                                                                                                                     aGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCCAGCAGTGCGCCGAt  <  8:219752/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GGGGCAGGCGTTGTGAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCGAGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCA‑GGGCCAAAGCCCATAACGGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCCAGCAGTGCGCCGAT  >  NC_002947/3116774‑3117070

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