Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,686,443	G→T	8.4%	G284V (GGC→GTC)	PP_3251 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,686,443	0	G	T	8.4%	217.8 / 8.2	71	G284V (GGC→GTC)	PP_3251	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/34);  new base T (4/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

CACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGTCTCATGGC  >  NC_002947/3686298‑3686591                                                                                                                                                  |                                                                                                                                                     cACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCAGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCg                                                                                                                                                    >  6:274203/1‑148 (MQ=255)  aCGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                                  <  7:322526/149‑1 (MQ=255)    gcccgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGCCTTGCAGGTCGAAc                                                                                                                                                 <  8:310102/148‑1 (MQ=255)    gcccgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGCCTTGCAGGTCGAAc                                                                                                                                                 >  7:310102/1‑148 (MQ=255)    gcccgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACt                                                                                                                                                >  7:18281/1‑149 (MQ=255)          cTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                                  >  6:281588/1‑141 (MQ=255)          cTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                                  <  5:281588/141‑1 (MQ=255)                 ttCATGAGCTGAACCAGCTAAGCCCCTCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCtt                                                                                                                                   <  8:202979/149‑1 (MQ=255)                 ttCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCtt                                                                                                                                   <  8:197641/149‑1 (MQ=255)                 ttCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCtt                                                                                                                                   <  7:258769/149‑1 (MQ=255)                     ggCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAAc                                                                                                                               <  6:311602/149‑1 (MQ=255)                     ggCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAAc                                                                                                                               <  8:149066/149‑1 (MQ=255)                      gcgcTGAACCAGCTAAGCCCCGCGCTTGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGTCCGGGGTATCGCCACGCGAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTg                                                                                                                                           <  5:30090/136‑1 (MQ=255)                      gcgcTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACGCTTTg                                                                                                                                           >  6:30090/1‑136 (MQ=255)                           gAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGc                                                                                                                         >  3:157455/1‑149 (MQ=255)                            aaCCAGCTAAGCCCCGCGCATGATCCGCTGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCg                                                                                                                        >  1:222275/1‑149 (MQ=255)                            aaCCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAg                                                                                                                          <  8:18281/147‑1 (MQ=255)                                      gCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGc                                                                                                              >  6:364803/1‑149 (MQ=255)                                      gCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGc                                                                                                              <  7:100663/149‑1 (MQ=255)                                       ccccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGACCGCAAGCgg                                                                                                                       <  3:227569/139‑1 (MQ=255)                                       ccccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCgg                                                                                                                       <  3:232358/139‑1 (MQ=255)                                       ccccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCgg                                                                                                                       >  4:232358/1‑139 (MQ=255)                                       ccccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCgg                                                                                                                       >  4:227569/1‑139 (MQ=255)                                              cATGATCTGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCTCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGat                                                                                                      <  5:274203/149‑1 (MQ=255)                                              cATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGat                                                                                                      <  8:32021/149‑1 (MQ=255)                                              cATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGat                                                                                                      <  2:371623/149‑1 (MQ=255)                                               aTGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGatg                                                                                                     <  4:157455/149‑1 (MQ=255)                                               aTGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGat                                                                                                      >  3:238583/1‑148 (MQ=255)                                                      gCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGcc                                                                                               >  4:138644/1‑148 (MQ=255)                                                       cAGGCGGTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGa                                                                                             <  2:140644/149‑1 (MQ=255)                                                          gCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCa                                                                                          <  3:138644/149‑1 (MQ=255)                                                            gTTCTACGCCAGGGCTCGCTTAGGCTACCGCCTGAATGAGCGGGGTATCGCCACGCCAGACCTGATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCAtt                                                                                        <  6:202932/149‑1 (MQ=255)                                                            gTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCAtt                                                                                        >  6:156950/1‑149 (MQ=255)                                                                 aCGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTgccgc                                                                                   >  4:126360/1‑149 (MQ=255)                                                                  cGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTgccgcc                                                                                  <  1:262582/149‑1 (MQ=255)                                                                     cAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTg                                                                                >  3:101139/1‑148 (MQ=255)                                                                           gcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCagcagc                                                                         >  4:194153/1‑149 (MQ=255)                                                                              cTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCa                                                                      <  5:364803/149‑1 (MQ=255)                                                                                 aGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGa                                                                   >  2:102534/1‑149 (MQ=255)                                                                                 aGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGa                                                                   >  8:25797/1‑149 (MQ=255)                                                                                          ccTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATcaatcaat                                                          >  2:368870/1‑149 (MQ=255)                                                                                                tGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtggctgg                                                    <  4:238583/149‑1 (MQ=255)                                                                                                   ccGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCTGCCTGCAGCAGCGCACGATCAATCAATGGCTGGcc                                                  >  1:19672/1‑148 (MQ=255)                                                                                                   ccGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCg                                                 >  6:60993/1‑149 (MQ=255)                                                                                                       ggTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAAcc                                             <  5:60993/149‑1 (MQ=255)                                                                                                          aTCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGCTGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCgg                                          >  5:199484/1‑149 (MQ=255)                                                                                                           tCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCgg                                          <  5:156950/148‑1 (MQ=255)                                                                                                           tCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGgc                                         <  7:25797/149‑1 (MQ=255)                                                                                                            cgccacgccaGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGgcg                                        >  5:153650/1‑149 (MQ=255)                                                                                                              ccacgccaGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGACCCCGgcgc                                       >  6:222768/1‑148 (MQ=255)                                                                                                               cacgccaGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCa                                     >  7:292660/1‑149 (MQ=255)                                                                                                                 cgccaGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATg                                   <  7:102816/149‑1 (MQ=255)                                                                                                                  gccaGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGc                                  <  8:292660/149‑1 (MQ=255)                                                                                                                    caGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGGGCCATGc                                  >  2:153738/1‑147 (MQ=255)                                                                                                                            tgcaaacAGTTCGCCTTGCAGGTCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGt                        >  1:103697/7‑149 (MQ=255)                                                                                                                            aTCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTAGATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGt                        >  1:105993/1‑149 (MQ=255)                                                                                                                              ccGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCgg                                                                                                                       >  4:341944/1‑52 (MQ=255)                                                                                                                              ccGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCgg                                                                                                                       <  3:341944/52‑1 (MQ=255)                                                                                                                                 ccAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCg                   >  1:199974/1‑149 (MQ=255)                                                                                                                                  cAGTTCGCCTTGCAGGTCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCg                   >  2:315203/1‑148 (MQ=255)                                                                                                                                     ttCGCCTTGCAGGTCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCAtt                           >  6:377813/1‑137 (MQ=255)                                                                                                                                     ttCGCCTTGCAGGTCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCAtt                           <  5:377813/137‑1 (MQ=255)                                                                                                                                       cGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTa              <  6:153650/148‑1 (MQ=255)                                                                                                                                        gACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTa              >  3:271833/1‑147 (MQ=255)                                                                                                                                        gACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACa            <  2:222275/149‑1 (MQ=255)                                                                                                                                         aCCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAg           >  1:318894/1‑149 (MQ=255)                                                                                                                                           cTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGtc         <  3:104920/149‑1 (MQ=255)                                                                                                                                             gCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCTATGGCTGGCCGAACCCGGCGCCATGccc                                <  1:219249/124‑1 (MQ=255)                                                                                                                                             gCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGccc                                >  2:219249/1‑124 (MQ=255)                                                                                                                                              cAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCTATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGtctc       <  1:102534/148‑1 (MQ=255)                                                                                                                                                aGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGTCTCATg    <  2:337280/149‑1 (MQ=255)                                                                                                                                                  gCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGTCTCATGGc  >  1:13325/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGTCTCATGGC  >  NC_002947/3686298‑3686591

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