Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,821,794	G→C	9.3%	D29E (GAC→GAG) ‡	ptxD ←	phosphonate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,821,794	0	G	C	9.3%	165.1 / 7.2	54	D29E (GAC→GAG) ‡	ptxD	phosphonate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (19/30);  new base C (3/2);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

GTCCACGCCTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCAGAGGTC  >  NC_002947/3821647‑3821942                                                                                                                                                     |                                                                                                                                                      gTCCACGCCTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGgtggtg                                                                                                                                                       <  1:206482/148‑1 (MQ=255)    cACGCCTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCa                                                                                                                                                   <  6:276835/149‑1 (MQ=255)         cTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAg                                                                                                                                             <  1:262900/149‑1 (MQ=255)         cTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAg                                                                                                                                             <  3:81353/149‑1 (MQ=255)            ccGAAACGCTCGACACCTCTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGTATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTc                                                                                                                                          <  1:54851/149‑1 (MQ=255)            ccGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGt                                                                                                                                           <  5:93317/148‑1 (MQ=255)            ccGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTc                                                                                                                                          <  2:11271/149‑1 (MQ=255)            ccGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCAAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGACACCA‑CCCAgc                                                                                                                                            <  2:164025/147‑2 (MQ=255)             cGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCCACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CACAGGTCa                                                                                                                                         >  6:243519/1‑149 (MQ=255)             cGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCa                                                                                                                                                   >  4:171701/1‑140 (MQ=255)             cGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCa                                                                                                                                                   <  3:171701/140‑1 (MQ=255)             cGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCC‑CCCAGGTCa                                                                                                                                         >  5:218176/1‑149 (MQ=255)               aaaCGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCAc                                                                                                                                        >  4:31199/1‑148 (MQ=255)                       gACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGAcgc                                                                                                                                <  7:34143/148‑1 (MQ=255)                         caccacTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGAt                                                                                                                             <  5:328955/149‑1 (MQ=255)                          accacTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATc                                                                                                                            >  7:129233/1‑149 (MQ=255)                               ttCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGCGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTg                                                                                                                         >  1:231219/1‑148 (MQ=255)                               ttCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGACACCA‑CCCAg                                                                                                                                             <  5:58647/127‑1 (MQ=255)                               ttCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGACACCA‑CCCAg                                                                                                                                             >  6:58647/1‑127 (MQ=255)                                 ccAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACTACCCGTGCGCCCCGGGTAACGCATCGCGCTGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑GCCAGGTCACCTCGACGCGATCTTGCAg                                                                                                                      <  2:238720/148‑1 (MQ=255)                                 ccAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCa                                                                                                                       <  3:58893/147‑1 (MQ=255)                                  cAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGt                                                                                                                                           <  7:233510/126‑1 (MQ=255)                                  cAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGt                                                                                                                                           >  8:233510/1‑126 (MQ=255)                                   aGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCg                                                                                                                   <  1:1670/149‑1 (MQ=255)                                    gTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGt                                                                                                                  <  8:215901/149‑1 (MQ=255)                                    gTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGt                                                                                                                  <  7:193190/149‑1 (MQ=255)                                                   cGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAgtcgtc                                                                                                    <  1:67608/149‑1 (MQ=255)                                                      gCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGa                                                                                                 >  6:127101/1‑148 (MQ=255)                                                      gCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGAc                                                                                                <  2:231219/149‑1 (MQ=255)                                                            tGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGt                                                                                          <  3:31199/149‑1 (MQ=255)                                                                aTCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTtat                                                                                      <  5:23934/149‑1 (MQ=255)                                                                 tCCATGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTtata                                                                                     <  4:92994/149‑1 (MQ=255)                                                                   cAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACa                                                                                   <  4:40143/149‑1 (MQ=255)                                                                      gcgcAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGa                                                                                <  5:126081/149‑1 (MQ=255)                                                                             gCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACCGGACGATTc                                                                          >  6:77737/1‑149 (MQ=255)                                                                                 gCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGc                                                                        <  5:117088/147‑1 (MQ=255)                                                                                     ccAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGcttctt                                                                   >  5:291979/1‑148 (MQ=255)                                                                                         caacCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTg                                                              >  4:256464/1‑149 (MQ=255)                                                                                            cccGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAg                                                           <  3:256464/149‑1 (MQ=255)                                                                                               gTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCa                                                                                                                                                   >  2:262213/1‑59 (MQ=255)                                                                                                      cccGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAAt                                                 >  3:268876/1‑149 (MQ=255)                                                                                                               gCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATgg                                        <  5:127101/149‑1 (MQ=255)                                                                                                               gCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATgg                                        >  3:222114/1‑149 (MQ=255)                                                                                                               gCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCACCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATgg                                        >  7:238860/1‑149 (MQ=255)                                                                                                                    gcgcAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGa                                   >  8:34886/1‑149 (MQ=255)                                                                                                                       cAGGCGGGCCAGGGCATCGGGCTGGCTGGTGTCCC‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGAcgc                                >  3:131271/1‑149 (MQ=255)                                                                                                                        aGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGAcgcg                               <  4:268876/149‑1 (MQ=255)                                                                                                                         ggCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCACCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGa                              <  8:238860/149‑1 (MQ=255)                                                                                                                          gCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAg                             >  6:200739/1‑149 (MQ=255)                                                                                                                             ggCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTa                          <  5:200739/149‑1 (MQ=255)                                                                                                                               ccAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATc                        <  4:222114/149‑1 (MQ=255)                                                                                                                                   ggCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCgcgc                    >  3:292027/1‑149 (MQ=255)                                                                                                                                           gCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCg            <  6:291979/149‑1 (MQ=255)                                                                                                                                            cTGGGTGGACACCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCgg           >  8:4911/1‑149 (MQ=255)                                                                                                                                             tGGGTGGACACCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGc          >  7:89965/1‑149 (MQ=255)                                                                                                                                                   gtgtCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCagag     <  2:323536/148‑1 (MQ=255)                                                                                                                                                     gtCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCAGAGGTc  >  7:124707/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GTCCACGCCTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAG‑GCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCA‑CCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCAGAGGTC  >  NC_002947/3821647‑3821942

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