Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,188,260	G→T	8.2%	A52A (GCG→GCT)	PP_3684 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,188,260	0	G	T	8.2%	161.7 / 2.6	56	A52A (GCG→GCT)	PP_3684	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/28);  new base T (3/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

TTTTTCGCGTGGTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCT  >  NC_002947/4188113‑4188402                                                                                                                                                    |                                                                                                                                               tttttCGCGTGGTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCgg                                                                                                                                               >  4:94481/1‑149 (MQ=255)   tttCGCGTGGTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGt                                                                                                                                              <  6:37701/148‑1 (MQ=255)      cgcgTGGTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGaa                                                                                                                                           <  2:265116/148‑1 (MQ=255)         gTGGTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGc                                                                                                                                       <  8:249849/149‑1 (MQ=255)            gTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCTGTTCACCGCCTc                                                                                                                                    <  6:191877/149‑1 (MQ=255)             tttAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATACGGCATGAGAGCCGATTCATCGAGGCGGTCAACAGCCTc                                                                                                                                    >  4:191651/1‑148 (MQ=255)             tttAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTc                                                                                                                                    >  2:161186/1‑148 (MQ=255)                       accaTGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCAATATACAAACTTTGGGATTACGTCTTCCCATACAGAAGACCAGAAGCAGAGATAATGCATGAGCGCCAATTCAACGAGGCTGTGAACAGCAACAAAACCTTGAg                                                                                                                         >  3:92375/1‑149 (MQ=255)                        ccaTGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAgg                                                                                                                        <  8:43904/149‑1 (MQ=255)                         caTGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAggg                                                                                                                       >  6:324982/1‑149 (MQ=255)                               cAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACt                                                                                                                  >  8:112728/1‑148 (MQ=255)                               cAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGAAAGAGCGCCAATTAATAGAGGCGGTAAAAAGCCTCAAAACCTTGAGGGTCaat                                                                                                                  >  7:284387/1‑146 (MQ=255)                                  aCCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCaaaa                                                                                                                                <  5:18228/131‑1 (MQ=255)                                  aCCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCaaaa                                                                                                                                >  6:18228/1‑131 (MQ=255)                                     atccatACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAg                                                                                                           >  3:46599/1‑149 (MQ=255)                                      tccatACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAgg                                                                                                          <  3:94481/149‑1 (MQ=255)                                        catACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGc                                                                                                         <  7:112728/148‑1 (MQ=255)                                        catACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCg                                                                                                        >  7:41699/1‑149 (MQ=255)                                         atACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCgg                                                                                                       <  7:116034/149‑1 (MQ=255)                                         atACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCgg                                                                                                       >  2:122539/1‑149 (MQ=255)                                           aCTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGCCCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGa                                                                                                     >  7:237820/1‑149 (MQ=255)                                             tCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCg                                                                                                                                                <  5:262204/104‑1 (MQ=255)                                             tCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCg                                                                                                                                                >  6:262204/1‑104 (MQ=255)                                             tCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCgg                                                                                                       <  2:59124/145‑1 (MQ=255)                                             tCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCgg                                                                                                       >  1:59124/1‑145 (MQ=255)                                               ggCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATg                                                                                                   <  8:41699/147‑1 (MQ=255)                                               ggCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGtc                                                                                                 <  8:122138/149‑1 (MQ=255)                                               ggCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTAACCGAGGCGGGATGtc                                                                                                 >  4:125483/1‑149 (MQ=255)                                               ggCAGGTGGAAGATGAAATCACTATAACAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAATGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATg                                                                                                   <  7:213880/147‑1 (MQ=255)                                                        aaGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGAccc                                                                                        <  1:122539/149‑1 (MQ=255)                                                            tGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGaa                                                                                    <  2:144605/149‑1 (MQ=255)                                                                  cACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGtcg                                                                                                >  7:225360/1‑131 (MQ=255)                                                                  cACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGtcg                                                                                                <  8:225360/131‑1 (MQ=255)                                                                         ccAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCGATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGGCGATCGACCCTGAAGAGATACgtgat                                                                        <  3:191651/148‑1 (MQ=255)                                                                         ccAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACgtgatg                                                                       <  8:38265/149‑1 (MQ=255)                                                                         ccAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACgtgatg                                                                       <  8:43235/149‑1 (MQ=255)                                                                                   gATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGaa                                                                                    >  6:261994/1‑126 (MQ=255)                                                                                   gATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGaa                                                                                    <  5:261994/126‑1 (MQ=255)                                                                                        cggcTTCCCAGACCGCAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGATGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGcc                                                        <  8:237820/146‑1 (MQ=255)                                                                                              cccAGACAGAAGAACATAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAg                                                   <  8:284387/148‑1 (MQ=255)                                                                                                cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTTAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTgaga                                                           <  7:134530/138‑1 (MQ=255)                                                                                                cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTgaga                                                           >  8:134530/1‑138 (MQ=255)                                                                                                     cGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCaa                                           <  5:324982/149‑1 (MQ=255)                                                                                                            cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGATCAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTg                                    >  1:223444/1‑149 (MQ=255)                                                                                                            cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATc                                                                  >  6:11868/1‑119 (MQ=255)                                                                                                            cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATc                                                                  <  5:11868/119‑1 (MQ=255)                                                                                                             aGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTAACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTgg                                   <  3:125483/149‑1 (MQ=255)                                                                                                             aGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCAATGACCGAGGCCGCATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACttcg                                   >  1:215291/1‑147 (MQ=255)                                                                                                               aaGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTc                                 >  7:289873/1‑149 (MQ=255)                                                                                                                aGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTc                                 <  8:289873/148‑1 (MQ=255)                                                                                                                            gCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATc                    >  3:307638/1‑149 (MQ=255)                                                                                                                            gCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATc                    >  1:35328/1‑149 (MQ=255)                                                                                                                                   cgcTAATTCATCTAGGCGGTGTACAGCCTCAAATCCTTGAGGGTCAGTGACCGAGGCGGGATGTCGATCGACCCTGATGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGGTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGc             <  2:215291/149‑1 (MQ=255)                                                                                                                                          tCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACTTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACaa       >  6:311779/1‑148 (MQ=255)                                                                                                                                            tttgAGGCTGTTCACCGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAAc    >  8:287902/4‑149 (MQ=255)                                                                                                                                            aTCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGCCCGAGGCGGGATGTCGCTCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAACGCGAACaaacc    >  1:246810/1‑147 (MQ=255)                                                                                                                                               gAGGCTGTTCACCGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCt  >  6:201859/1‑148 (MQ=255)                                                                                                                                               gAGGCTGTGAACAGCCTCAAAACCTTGAGGTTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCt  <  7:104232/148‑1 (MQ=255)                                                                                                                                               gAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCt  <  7:108685/148‑1 (MQ=255)                                                                                                                                                    ggTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCaa                                           >  5:75676/1‑102 (MQ=255)                                                                                                                                                    ggTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCaa                                           <  6:75676/102‑1 (MQ=255)                                                                                                                                                    ggTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGt                                  >  8:231157/1‑111 (MQ=255)                                                                                                                                                    ggTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGt                                  <  7:231157/111‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               TTTTTCGCGTGGTTTAACATCCACCATGACCAAACCATCCATACTCGGCAGGTGGAAGATGAAATCACTATACCAACTTTGGGATTACGTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCT  >  NC_002947/4188113‑4188402

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑