Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,080,154	C→A	5.6%	A33E (GCA→GAA) ‡	PP_5332 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,080,154	0	C	A	5.6%	295.8 / 5.7	90	A33E (GCA→GAA) ‡	PP_5332	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (41/44);  new base A (2/3);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

GCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAGGC  >  NC_002947/6080009‑6080298                                                                                                                                                  |                                                                                                                                                  gCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGca                                                                                                                                                  >  4:160790/1‑147 (MQ=255)       gCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCt                                                                                                                                         >  1:385143/1‑149 (MQ=255)        cccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCTAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTc                                                                                                                                           <  4:282234/146‑1 (MQ=255)        cccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGTGCGCCCTGCTGGGCGAT‑CAgga                                                                                                                                            <  1:152624/145‑2 (MQ=255)        cccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCTGCTGGGCGAT‑CAgga                                                                                                                                            <  7:279326/145‑2 (MQ=255)          cGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAg                                                                                                                                      >  6:358898/1‑149 (MQ=255)             gTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTa                                                                                                                                   <  6:291475/149‑1 (MQ=255)                 aaGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGca‑c                                                                                                                                                <  4:259233/132‑1 (MQ=255)                 aaGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGca‑c                                                                                                                                                >  3:259233/1‑132 (MQ=255)                     cccACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAgg                                                                                                                           <  4:202075/149‑1 (MQ=255)                      ccACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCC‑CCGGTCCTCAGCTACCAGCAgg                                                                                                                           >  7:341731/1‑148 (MQ=255)                      ccACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAg                                                                                                                                      <  1:180412/137‑1 (MQ=255)                      ccACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAg                                                                                                                                      >  2:180412/1‑137 (MQ=255)                        aCTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAg                                                                                                                        <  8:236565/149‑1 (MQ=255)                            gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCt                                                                                                                                         <  1:118739/128‑1 (MQ=255)                            gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCt                                                                                                                                         >  2:118739/1‑128 (MQ=255)                            gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGc                                                                                                                      <  6:196119/147‑1 (MQ=255)                              aatcaatGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAAc                                                                                                                  >  1:261050/1‑149 (MQ=255)                               atcaatGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCTAACTCAAATGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCACCCAGCAGGGCGCA‑CAGGTCCTCATCTACCAGCAGGTAGGCGaaac                                                                                                                 >  8:184031/1‑147 (MQ=255)                               atcaatGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAAcc                                                                                                                 >  3:313337/1‑149 (MQ=255)                                    tGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACcgccgc                                                                                                            >  1:217305/1‑149 (MQ=255)                                    tGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACcgccgc                                                                                                            >  1:271774/1‑149 (MQ=255)                                      agaagaACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgctg                                                                                                          >  1:183498/1‑149 (MQ=255)                                       gaagaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgctgc                                                                                                         >  7:327370/1‑149 (MQ=255)                                        aagaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgctgct                                                                                                        >  3:279399/1‑149 (MQ=255)                                          gaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCTAACTCCTGATCGCCCATCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCtgt                                                                                                      <  6:184707/149‑1 (MQ=255)                                          gaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCtgt                                                                                                      >  5:189318/1‑149 (MQ=255)                                           aaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCtgtg                                                                                                     <  4:358397/149‑1 (MQ=255)                                              gCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCt                                                                                                  <  5:358898/149‑1 (MQ=255)                                                  tGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGc                                                                                              <  2:217305/149‑1 (MQ=255)                                                      tACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGAc                                                                                           >  1:177786/1‑148 (MQ=255)                                                        ccTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAg                                                                                        <  2:261050/149‑1 (MQ=255)                                                        ccTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAg                                                                                        <  5:318359/149‑1 (MQ=255)                                                          tttCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGCACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGc                                                                                      >  7:171797/1‑149 (MQ=255)                                                            tCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGc                                                                                      <  7:339009/147‑1 (MQ=255)                                                              aGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCAt                                                                                   <  6:189318/148‑1 (MQ=255)                                                                   gTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTcc                                                                              <  5:12312/148‑1 (MQ=255)                                                                   gTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTcc                                                                              <  4:313337/148‑1 (MQ=255)                                                                    tAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGc                                                                            >  7:332858/1‑149 (MQ=255)                                                                    tAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGc                                                                            <  3:160790/149‑1 (MQ=255)                                                                     aaCCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCa                                                                                                                                       >  5:140939/1‑89 (MQ=255)                                                                     aaCCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCa                                                                                                                                       <  6:140939/89‑1 (MQ=255)                                                                     aaCCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGc                                                                            >  6:101284/1‑148 (MQ=255)                                                                         gAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAggg                                                                       <  5:101284/149‑1 (MQ=255)                                                                            cAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCaa                                                                    <  4:279399/149‑1 (MQ=255)                                                                              aCATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGt                                                                  <  7:184031/149‑1 (MQ=255)                                                                              aCATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGt                                                                  <  4:356761/149‑1 (MQ=255)                                                                                 tGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCgg                                                               >  5:242269/1‑149 (MQ=255)                                                                                   gCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTa                                                             <  7:326129/149‑1 (MQ=255)                                                                                      aaCTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTAcgc                                                          <  8:332858/149‑1 (MQ=255)                                                                                        cTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACgcgc                                                        >  1:139725/1‑149 (MQ=255)                                                                                         tCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTa                                                             <  4:220568/143‑1 (MQ=255)                                                                                         tCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTa                                                             >  3:220568/1‑143 (MQ=255)                                                                                            aCTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCggg                                                    <  6:114920/149‑1 (MQ=255)                                                                                              tgctgGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCgggg                                                   >  5:208680/1‑148 (MQ=255)                                                                                                ctgGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTAc                                                <  2:220665/149‑1 (MQ=255)                                                                                                        tCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCGGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGc                                        >  2:63361/1‑149 (MQ=255)                                                                                                        tCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGc                                        <  8:327370/149‑1 (MQ=255)                                                                                                        tCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGc                                        >  1:308313/1‑149 (MQ=255)                                                                                                        tCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGc                                        >  2:117579/1‑149 (MQ=255)                                                                                                           gcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTgg                                     >  5:157431/1‑149 (MQ=255)                                                                                                             agTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGc                                    >  7:260299/1‑148 (MQ=255)                                                                                                              gTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCa                                  >  5:155226/1‑149 (MQ=255)                                                                                                               tGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCAt                                 >  4:22037/1‑149 (MQ=255)                                                                                                                 aCCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATg                                <  4:253747/148‑1 (MQ=255)                                                                                                                    gTGCCGAACTCCTGATCGCCCAGCAGGGCGCAGCAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGtt                             >  1:242935/1‑149 (MQ=255)                                                                                                                    gTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGttt                            <  6:157431/149‑1 (MQ=255)                                                                                                                     tGCCGAAGTCCTGATCGCCCAGCAGGGCGCT‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTg                           <  2:373112/149‑1 (MQ=255)                                                                                                                     tGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCTGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTg                           <  1:63361/149‑1 (MQ=255)                                                                                                                         gAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCaa                       <  6:242269/149‑1 (MQ=255)                                                                                                                          aaCTCCTGATCGCCCCGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAAc                      <  8:341731/149‑1 (MQ=255)                                                                                                                          aaCTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAAc                      <  2:271774/149‑1 (MQ=255)                                                                                                                            cTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCt                    >  7:287289/1‑149 (MQ=255)                                                                                                                               ctgtgcGCCCTGCTGGGCGAT‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTg                   >  1:212198/6‑147 (MQ=255)                                                                                                                               cTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTg                   <  2:183498/147‑1 (MQ=255)                                                                                                                               cTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTg                   <  6:208680/147‑1 (MQ=255)                                                                                                                                  aTCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCg                                              >  4:44831/1‑117 (MQ=255)                                                                                                                                  aTCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCg                                              <  3:44831/117‑1 (MQ=255)                                                                                                                                     gCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCaa           <  2:308313/149‑1 (MQ=255)                                                                                                                                      cccAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAAc                      <  4:50670/137‑1 (MQ=255)                                                                                                                                      cccAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAAc                      >  3:50670/1‑137 (MQ=255)                                                                                                                                      cccAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAc          >  5:72187/1‑149 (MQ=255)                                                                                                                                       ccAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAcc         <  4:85675/149‑1 (MQ=255)                                                                                                                                       ccAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAcc         >  4:325151/1‑149 (MQ=255)                                                                                                                                       ccAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAcc         >  3:85675/1‑149 (MQ=255)                                                                                                                                         tgctgggCGAT‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGttt                            >  1:326625/5‑128 (MQ=255)                                                                                                                                         tgctgggCGAT‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGttt                            <  2:326625/124‑1 (MQ=255)                                                                                                                                         agcagGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCt       <  6:155226/149‑1 (MQ=255)                                                                                                                                            agGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAg    >  7:175614/1‑149 (MQ=255)                                                                                                                                             gggCGCA‑CACGTCCTCAGCTACCAGCAGGTAGGCGAACCGCAGCTGATGTGGCTGCGCTACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGATGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCCgg   >  2:239498/1‑149 (MQ=255)                                                                                                                                              ggCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAgg   <  2:216805/148‑1 (MQ=255)                                                                                                                                              ggCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAGGc  >  7:225163/1‑149 (MQ=255)                                                                                                                                              ggCGCA‑CAGGTCCTCAGCTACCAGCAGGAAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAgg   >  1:216805/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCA‑CAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAGGC  >  NC_002947/6080009‑6080298

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