Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,399,257	G→A	8.9%	R79H (CGC→CAC)	PP_3004 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,399,257	0	G	A	8.9%	158.0 / 4.5	56	R79H (CGC→CAC)	PP_3004	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/28);  new base A (2/3);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

TCGACCTGGAGCAGGCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGGCAGACCGGTCGGCGAGCCACCGGCT  >  NC_002947/3399110‑3399404                                                                                                                                                     |                                                                                                                                                    tCGACCTGGAGCAGGCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGCTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGc                                                                                                                                                    >  3:372185/1‑149 (MQ=255)  cGACCTGGAGCAGGCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGc                                                                                                                                                    <  6:298831/148‑1 (MQ=255)    aCCTGGAGCAGGCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGAGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCACCt                                                                                                                                                  >  4:291907/1‑148 (MQ=255)       tGGAGCAGGCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAg                                                                                                                                              >  1:288620/1‑149 (MQ=255)         gAGCAGGCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGTCTGGTCCGGGTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAg                                                                                                                                              <  4:77697/147‑1 (MQ=255)         gAGCAGGCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAgg                                                                                                                                             <  5:73571/148‑1 (MQ=255)         gAGCAGGCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAgg                                                                                                                                             <  5:198747/148‑1 (MQ=255)         gAGCAGGCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAg                                                                                                                                              <  7:51645/147‑1 (MQ=255)         gAGCAGGCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAg                                                                                                                                              <  4:293658/147‑1 (MQ=255)               gCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGCTGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCcg                                                                                                                                      >  8:78867/1‑149 (MQ=255)                 aaTGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCcgcg                                                                                                                                    >  5:247620/1‑149 (MQ=255)                 aaTGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCcgcg                                                                                                                                    >  6:84907/1‑149 (MQ=255)                  aTGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCGCGGATCACCTCAAGGCGATcccgc                                                                                                                                    <  6:235433/149‑4 (MQ=255)                  aTGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCgcgc                                                                                                                                   <  1:215639/149‑1 (MQ=255)                     aGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCACCTCAAGGCGATcccgccct                                                                                                                                 <  6:362425/148‑7 (MQ=255)                     aGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGAGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCACCTCAAGGCGATcccgccct                                                                                                                                 <  3:291907/148‑7 (MQ=255)                          cGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCa                                                                                                                            >  5:225425/1‑148 (MQ=255)                            aaGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAg                                                                                                                         >  1:76365/1‑149 (MQ=255)                              gCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGt                                                                                                                       <  1:361314/149‑1 (MQ=255)                                     aaGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCt                                                                                                                >  4:327832/1‑149 (MQ=255)                                     aaGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCt                                                                                                                >  4:297324/1‑149 (MQ=255)                                            gTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGaa                                                                                                          <  7:282295/148‑1 (MQ=255)                                            gTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGaa                                                                                                          <  4:253542/148‑1 (MQ=255)                                             tCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAg                                                                                                         <  4:372185/148‑1 (MQ=255)                                               ccAAAGCTGGCTGGTCCTGTTCTCAATTACAAAAATGGTGTCAGGACGTTATCGGATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGcc                                                                                                      <  6:247620/149‑1 (MQ=255)                                               ccAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGcc                                                                                                      >  2:310540/1‑149 (MQ=255)                                                            gTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGtt                                                                                           <  2:86724/147‑1 (MQ=255)                                                             tCCGGTTCTAAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTg                                                                                        <  7:307804/149‑1 (MQ=255)                                                                     tcAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTg                                                                                >  2:324821/1‑149 (MQ=255)                                                                       aaTTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTc                                                                              <  3:295138/149‑1 (MQ=255)                                                                         ttACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCgt                                                                            <  6:225425/149‑1 (MQ=255)                                                                         ttACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCgt                                                                            <  1:205883/149‑1 (MQ=255)                                                                                       tCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACttt                                                               <  3:251308/148‑1 (MQ=255)                                                                                        cAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCa                                                             >  6:280332/1‑149 (MQ=255)                                                                                         aGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCCCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCAc                                                            >  1:309501/1‑149 (MQ=255)                                                                                              gTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACcgcgc                                                       <  7:342311/149‑1 (MQ=255)                                                                                                  tCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGAc                                                    >  7:349694/1‑148 (MQ=255)                                                                                                  tCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGAc                                                    <  3:297324/148‑1 (MQ=255)                                                                                                  tCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGAc                                                    >  7:229019/1‑148 (MQ=255)                                                                                                    gCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCaa                                                 <  5:280332/149‑1 (MQ=255)                                                                                                    gCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCaa                                                 <  5:84907/149‑1 (MQ=255)                                                                                                     cATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAAt                                                <  7:78867/149‑1 (MQ=255)                                                                                                     cATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAAt                                                >  1:65626/1‑149 (MQ=255)                                                                                                     cATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAAt                                                >  3:351407/1‑149 (MQ=255)                                                                                                       tCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACAATAGGTCCAGCCATGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATg                                               >  3:123341/1‑148 (MQ=255)                                                                                                                   gagagaACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGcc                                  <  1:310540/149‑1 (MQ=255)                                                                                                                        aaCGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTc                             >  1:292409/1‑149 (MQ=255)                                                                                                                              aTCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGGCAGa                       <  2:76365/149‑1 (MQ=255)                                                                                                                               tCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGGCAGAc                      >  3:196666/1‑149 (MQ=255)                                                                                                                               tCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGGCAGAc                      <  1:291264/149‑1 (MQ=255)                                                                                                                                 acaGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGGCAGACCg                    >  7:201308/1‑149 (MQ=255)                                                                                                                                       gCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGGCAGACCGGTCGGc              >  3:196656/1‑149 (MQ=255)                                                                                                                                       gCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGGCAGACCGGTCGGc              >  6:297441/1‑149 (MQ=255)                                                                                                                                         aGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGGCAGACCGGTCGGCGa            <  2:288620/149‑1 (MQ=255)                                                                                                                                              gcggATCACCTCAAGGCGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGGCAGACCGGTCGGCGAGCCa        >  5:327378/3‑149 (MQ=255)                                                                                                                                                ggATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGTAAGACCGGTCGGCGAGCCAcc      <  2:292409/149‑1 (MQ=255)                                                                                                                                                     gCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGGCAGACCGGTCGGCGAGCCACCGGCt  <  1:324821/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TCGACCTGGAGCAGGCAATGAGCTTCGAAGCAGGCCAAGCCTTGTCCCAAAGCGGGCTGGTCCGGTTCTCAATTACCAAAATGGTGTCAGGACGTTATCGCATCAGCGATGAAGGAGAGAACGCCATCACAGAGGCAGGGCG‑GGATCGCCTTGAGGTGATCCGCGCCTGACCATAGGTCCAGCCTTGCGAAGGCCGCAGTCGCGTTCTGCGGCCCTGTCGTGTCATATCACTTTCACCGCGCGACCAATGAACTTCACCGGCCCTGTCGGCAGACCGGTCGGCGAGCCACCGGCT  >  NC_002947/3399110‑3399404

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