Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,125,097	C→T	69.4%	E572K (GAG→AAG)	priA ←	primosome factor N'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,125,097	0	C	T	69.4%	60.2 / 32.0	57	E572K (GAG→AAG)	priA	primosome factor N'
Reads supporting (aligned to +/- strand):  ref base C (8/9);  new base T (17/21);  total (26/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 = 8.17e-01

CGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGC  >  NC_000913/4124956‑4125235                                                                                                                                                |                                                                                                                                             cGGCAAAGGCGTCGCAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGACGGTGCGTTTGCAGCACCACTTCGCCCCGTTTACCCGCACGCCCGGCACGACCGGCGACCTTGGTGTAAAGCT‑GAGCGAAACGc‑tc‑t                                                                                                                                            <  1:216826/143‑1 (MQ=255)  ggCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGAGTAAAGAT‑GAGCGAACCGC‑TC‑Tg                                                                                                                                           >  2:586920/1‑143 (MQ=255)     aaaGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAAAAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑Cg                                                                                                                                       >  1:154837/1‑143 (MQ=255)       agtcgtcgtAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAg                                                                                                                                     <  1:660133/140‑1 (MQ=255)               tAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATc                                                                                                                             <  1:587824/143‑1 (MQ=255)                aGCCTTTATAGAGCAACGTTTGCAACAGCGCACGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCg                                                                                                                            <  2:429584/143‑1 (MQ=255)                  ccTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCAGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAGAGCT‑GAGCGAAACGC‑TT‑TGC‑CGCGCGAAAATTGGc                                                                                                                          >  1:528743/1‑143 (MQ=255)                   cTTTATAGAGCAACGTTTGCACCAGCGGCTGTTCCGGATGCTGCGTTTGCAGCACCCCTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCAGGAACTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGGACATCGGCa                                                                                                                         >  1:407030/1‑143 (MQ=255)                   cTTCATAGAGCAACGTTTGCAACAGAGGATGTTGCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCa                                                                                                                         <  2:260338/143‑1 (MQ=255)                      tataGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTGCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑TAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGaa                                                                                                                      <  1:138198/143‑1 (MQ=255)                       ataGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGGCCTGGGTGTAAAGCT‑GAGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGaaa                                                                                                                     <  1:345393/143‑1 (MQ=255)                       ataGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGCCCCGGGCGTACAGCT‑GAGCGACACGC‑TC‑TGC‑CGCGCCCAACTCGGCAGgaa                                                                                                                     >  2:530038/1‑143 (MQ=255)                             cAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCCGGGTGTAAAGCT‑GCGCCAGACGC‑TT‑TGC‑CGACCGAAAATAGGCAGAAAACAgcg                                                                                                               >  2:483644/1‑143 (MQ=255)                               aCGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCATCACTTCGCGCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAgcgcg                                                                                                             <  2:325220/143‑1 (MQ=255)                                cGTTTGCAACAGATGCTGTTCTGGATGGTGCGGTTGCCGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAgcgcgc                                                                                                            <  2:197268/143‑1 (MQ=255)                                  tttGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGAGCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGACAATCGGCAGAAAACAGCGCGCcg                                                                                                          >  1:327429/1‑143 (MQ=255)                                    ggCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCcgtc                                                                                                        <  1:701603/142‑1 (MQ=255)                                      cAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCcgtcca                                                                                                      >  1:582394/1‑143 (MQ=255)                                       aaCAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGACACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCcgtccac                                                                                                     >  1:393103/1‑143 (MQ=255)                                            cggctgtgCCGGATGGTGCGTTTGCAGCACCACTGCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCcgtccacgtcca                                                                                                <  2:384368/139‑1 (MQ=255)                                            agGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCcgtccacgtcca                                                                                                <  2:154837/143‑1 (MQ=255)                                                   tCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTCACCCGCACGCCCGGCACGACCG‑‑‑‑CCGGAGCGTTATGTTGGTGGGAAAGACATCTTGTTCTTTCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGc                                                                                         <  1:188105/143‑1 (MQ=255)                                                        aTGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCa                                                                                    <  1:600573/143‑1 (MQ=255)                                                         tGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAg                                                                                   <  1:25832/143‑1 (MQ=255)                                                         tGGTGCGTTTGCAGCACCACTTCGCACTGTTTACCCGCACGTCCGGTCCGACCGGCGACCTGGGTGTATAGCT‑GAGCGACACGC‑TC‑TGC‑CGAGCGTCAATCGGCGGAAAACAGCGCGCCGTCCACGTCCAGTAATGCACCCAg                                                                                   >  2:338692/1‑143 (MQ=255)                                                          tgTGCCTCTGCAGCCTCACTCTGCCCTCTTTACCCGCACGCCCGGCACGACCGTCGACCTGGGTTTAAATCT‑GAGCTCAACGC‑TC‑CGC‑CGAGCGAAAATCGGCAGCAAACAGCCCGCCGTCCACGTCCATTAATGCAACCAGc                                                                                  <  2:53440/142‑1 (MQ=255)                                                          cgTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACACCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGc                                                                                  <  2:436676/142‑1 (MQ=255)                                                            tGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGt                                                                                >  2:175381/1‑143 (MQ=255)                                                            tGCGTTGGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑AAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGt                                                                                <  2:96418/143‑1 (MQ=255)                                                             gCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCGGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTc                                                                               >  1:297142/1‑143 (MQ=255)                                                             gCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTc                                                                               >  2:151526/1‑143 (MQ=255)                                                             gCGTTTGCAGCACCACTTCGCCCGGTTTACGAGCACGCCCGGCACTACCGGCGACCTGGGTGGAAAGCT‑GAGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGGAGTGCAGCCGGCGTc                                                                               >  1:602305/1‑143 (MQ=255)                                                              cGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTca                                                                              >  2:99352/1‑143 (MQ=255)                                                                          cacTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCG‑GAGCGAAACGC‑TC‑GGC‑CGAGCGAACAGCGGCAGAACACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGCCCCACCCGGGcaac                                                                  >  1:786632/1‑142 (MQ=255)                                                                             ttCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGCGCGACAATCGGCAGAAAACAGCGCGCCGTCCGCGCCCAGTACGCCAACCAGCGCCACACCCgggccaccgt                                                               >  2:356675/1‑136 (MQ=255)                                                                             ttCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGCGTAAAGCT‑GAGCGAAACGC‑TA‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACACCCGGGCAACgga                                                               >  2:156746/1‑142 (MQ=255)                                                                               cGCCCCGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGa                                                             <  1:55546/143‑1 (MQ=255)                                                                                gCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGAc                                                            <  2:393103/143‑1 (MQ=255)                                                                                gCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCCCGTCCAGCACTGCACCCAGCGTCCCATCCGGGACAGGGTgcc                                                            >  2:105995/1‑141 (MQ=255)                                                                                gCCCTGTTTACCCGCACGCCCCGCACGACCGGCGACCTGGGTGTAAAGCT‑GCGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGAc                                                            <  2:593708/143‑1 (MQ=255)                                                                                         aCCCGACCGCCCGGCACGACCGGCGACCTGTGTTTAAAGCG‑TAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGCAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCa                                                   <  2:316097/143‑1 (MQ=255)                                                                                           ccGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGc                                                 >  1:756096/1‑143 (MQ=255)                                                                                               aCGGCCGGGGCGACACGCGACTTGGGTGTAAAGCT‑GAGCGGACCGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAGACAGAGCGCCGTCCACGGCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCAttt                                             <  2:407030/143‑1 (MQ=255)                                                                                               aCGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCAttt                                             <  1:99352/143‑1 (MQ=255)                                                                                                cGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAGATGGTGACCCTTCGCCAGCAGTtg                                            >  2:630138/1‑143 (MQ=255)                                                                                                 gCCCGGCACGACCGGCGACCTGGGTGTACAGCT‑GAGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTtgt                                           <  1:716419/143‑1 (MQ=255)                                                                                                   ccGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGCAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACACCCGGGACATGGTGACCTTCAGCCAGCCTTTgtgt                                         >  2:469276/1‑143 (MQ=255)                                                                                                   ccGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCCACCAACGTCACATCCGGGAACTGGTGCCCTCTCGCCAGCAtgcgcgt                                         >  2:553285/1‑137 (MQ=255)                                                                                                    cGGCCCGAGCGGCGGCCTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTa                                        <  2:731816/143‑1 (MQ=255)                                                                                                               gcaaCCTGGGTGTAAGGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAGAATCGGCAGAAAACAGCGCGCCGTCCACGCCCAGTCATGCAACCAGCGTCACATCCGGGAAATGGCGACCTCTAGCCCGCAACCGTGTCCCAATCCAAAt                             >  2:688284/4‑143 (MQ=255)                                                                                                                   ccTGGGTGTAAAGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCcgc                         <  1:271454/143‑1 (MQ=255)                                                                                                                              aGCT‑GAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCAAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGAGCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCcgcg              >  1:250775/1‑143 (MQ=255)                                                                                                                                cT‑GGGCGCAACGC‑GT‑TGC‑CGAGCGAAAAGCGGCAGACATCGGCGCGCCGCCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGGCCTTGCGCCAGCATTGGTGTACCAATCAAAATCCGCGCGCCGCCGCGAt            <  2:200480/143‑1 (MQ=255)                                                                                                                                   gAGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGt          <  1:638509/143‑1 (MQ=255)                                                                                                                                    aGCGGAACGC‑TC‑TGC‑CGGGCGAAAATCGGCAGAAAACAGCGCGCCGGCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTa         <  2:603214/143‑1 (MQ=255)                                                                                                                                    aGCGAAACGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGCCACATCCGGGAAACGGTGACCTTTCGCCAGCATTTGCGTACCGACCAAAATCCCCGCCCCGCCGCGACGTa         >  2:205803/1‑143 (MQ=255)                                                                                                                                        aaaCGC‑TG‑TGC‑GGAGCCAATATCGCCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGTGATGTACTTc     <  2:304025/143‑1 (MQ=255)                                                                                                                                           cGC‑TT‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGc  <  1:175381/143‑1 (MQ=255)                                                                                                                                                |                                                                                                                                             CGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCT‑GAGCGAAACGC‑TC‑TGC‑CGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGC  >  NC_000913/4124956‑4125235

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑