Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,300,860	A→G	6.0%	T9A (ACC→GCC)	PP_3774 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,300,860	0	A	G	6.0%	196.9 / 2.3	66	T9A (ACC→GCC)	PP_3774	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base G (2/2);  total (34/32)
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.56e-01

AGCTATTGACCTAGATTTAGCTTGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGA  >  NC_002947/4300740‑4301007                                                                                                                         |                                                                                                                                                     aGCTATTGACCTAGATTTAGCTTGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGtctc                                                                                                                             >  1:42688/1‑146 (MQ=255) aGCTATTGACCTAGATTTAGCTTGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGtctc                                                                                                                             <  2:42688/146‑1 (MQ=255)      ttGACCTAGATTTAGCTTGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTcaac                                                                                                                       <  3:234144/147‑1 (MQ=255)          ccTAGATTTAGCTTGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAAccc                                                                                                                 >  6:261959/1‑149 (MQ=255)          ccTAGATTTAGCTTGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAAccc                                                                                                                 >  8:183967/1‑149 (MQ=255)           cTAGATTTAGCTTGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTcaaccaac                                                                                                                   <  1:275658/146‑1 (MQ=255)           cTAGATTTAGCTTGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTcaaccaac                                                                                                                   >  2:275658/1‑146 (MQ=255)              gATTTAGCTTGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGt                                                                                                              <  8:138687/148‑1 (MQ=255)                      ttGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTg                                                                                                       >  1:142925/1‑147 (MQ=255)                      ttGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTg                                                                                                       >  1:292009/1‑147 (MQ=255)                        gAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGa                                                                                                     <  5:256790/147‑1 (MQ=255)                         aGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGc                                                                                                   <  2:292009/148‑1 (MQ=255)                         aGGCTAAGATTGCGTGCACCCGTCGCTTGTACGGCGGCTATGCCACGAAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCt                                                                                                  <  8:250357/149‑1 (MQ=255)                               aGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACg                                                                                             <  4:129543/148‑1 (MQ=255)                                   tGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGc                                                                                        <  7:183967/148‑1 (MQ=255)                                      gTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTg                                                                                      >  3:41926/1‑147 (MQ=255)                                       ttCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCGCCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAg                                                                                   >  7:183192/1‑149 (MQ=255)                                       ttCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCGCCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAg                                                                                   >  5:255026/1‑149 (MQ=255)                                        tCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACTATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGc                                                                                  >  6:272217/1‑149 (MQ=255)                                         cAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGc                                                                                  <  4:259855/148‑1 (MQ=255)                                            ccGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGAtt                                                                              >  3:98141/1‑149 (MQ=255)                                              gTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGc                                                                            >  1:29547/1‑149 (MQ=255)                                                   ttGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGcc                                                                                                                                  >  7:258537/1‑91 (MQ=255)                                                   ttGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGcc                                                                                                                                  <  8:258537/91‑1 (MQ=255)                                                   ttGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCt                                                                       >  6:167828/1‑149 (MQ=255)                                                       aCGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGcc                                                                           >  7:82700/1‑141 (MQ=255)                                                       aCGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGcc                                                                           <  8:82700/141‑1 (MQ=255)                                                       aCGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTcc                                                                   >  7:50986/1‑149 (MQ=255)                                                       aCGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTc                                                                    >  4:2800/1‑148 (MQ=255)                                                        cggcggCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTcc                                                                   <  3:2800/148‑1 (MQ=255)                                                        cgccggcTATGCCACTCAATCCTTTTCGCCAGGACTTCCCCATGGCCGTAGAAATCAGCCACCTCACCATCGCCACGCCGGTAGCCGCCTCCTCAACCTACCCTGTTGCTCTGGAGCTGAACCG‑TGCTGAAGCGATTGCCCGCTTTcc                                                                   <  2:41506/145‑1 (MQ=255)                                                                 tGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCGCCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGt                                                          <  6:255026/148‑1 (MQ=255)                                                                 tGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCGCCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGt                                                          <  8:183192/148‑1 (MQ=255)                                                                 tGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTg                                                         <  4:231739/149‑1 (MQ=255)                                                                 tGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACACTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCATCTTTCCCAGCAGCGt                                                          >  8:260496/1‑148 (MQ=255)                                                                  gCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGa                                                        >  7:284719/1‑149 (MQ=255)                                                                    cACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACg                                                      >  2:153959/1‑149 (MQ=255)                                                                         aaaCCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTcc                                                                   <  4:116243/131‑1 (MQ=255)                                                                         aaaCCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCt                                                 >  5:169979/1‑149 (MQ=255)                                                                          aaCCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCt                                                 >  2:98930/1‑148 (MQ=255)                                                                          aaCCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTc                                                <  4:98141/149‑1 (MQ=255)                                                                           aCCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGGCGGTGCTCa                                               >  6:203424/1‑149 (MQ=255)                                                                              ttttCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑t                                                                                          <  4:7449/103‑1 (MQ=255)                                                                              ttttCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑t                                                                                          >  3:7449/1‑103 (MQ=255)                                                                              ttttCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGc                                             >  2:15519/1‑148 (MQ=255)                                                                                ttCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCg                                            <  2:142925/147‑1 (MQ=255)                                                                                     cAGGAGTTCCCCATGGCCGTAGATATCCGCCAGCTCACCTTCGCCACGCCGGTAGCCGGCTCCTCACCCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCGCAGCGATGGCTc                                     <  7:260496/149‑1 (MQ=255)                                                                                      acgAGTTCCCCATGGCCGTAGATATCACCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGGCTGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTc                                     <  1:98930/147‑1 (MQ=255)                                                                                      aGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCg                                                                                 <  2:216949/104‑1 (MQ=255)                                                                                      aGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCg                                                                                 >  1:216949/1‑104 (MQ=255)                                                                                             ccccATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGt                             >  8:60054/1‑149 (MQ=255)                                                                                                cATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTtct                          >  8:78153/1‑149 (MQ=255)                                                                                                cATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTtct                          <  7:78153/149‑1 (MQ=255)                                                                                                    gCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGc                      <  3:46281/149‑1 (MQ=255)                                                                                                     ccGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGcc                     <  8:284719/149‑1 (MQ=255)                                                                                                          gATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCccaac                <  1:153959/149‑1 (MQ=255)                                                                                                             atCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCccaacca              >  2:260009/1‑148 (MQ=255)                                                                                                             atCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCccaacca              >  3:27628/1‑148 (MQ=255)                                                                                                               cAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAg           >  7:112296/1‑149 (MQ=255)                                                                                                                aGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCt                                      >  6:19834/1‑121 (MQ=255)                                                                                                                aGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCt                                      <  5:19834/121‑1 (MQ=255)                                                                                                                 gCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGAt                                          >  4:89948/1‑116 (MQ=255)                                                                                                                 gCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGAt                                          <  3:89948/116‑1 (MQ=255)                                                                                                                 gCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAggg         >  1:147546/1‑149 (MQ=255)                                                                                                                     ccTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCtt     <  4:41926/149‑1 (MQ=255)                                                                                                                      cTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGGCGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTc    <  5:203424/149‑1 (MQ=255)                                                                                                                         aCCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGa  <  1:15519/148‑1 (MQ=255)                                                                                                                         |                                                                                                                                                     AGCTATTGACCTAGATTTAGCTTGAGGCTAAGTTTGCGTTCAGCCGTCGCTTGTACGGCGGCTATGCCACTCAAACCTTTTCGCCAGGAGTTCCCCATGGCCGTAGATATCAGCCACCTCACCATCGCCACGCCGGTAGCCGTCTCCTCAACCAACCCTGTTGCTCTGGAGCTGAACGG‑TGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGA  >  NC_002947/4300740‑4301007

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