Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,390,969	A→T	6.3%	S354T (TCC→ACC)	tldD ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,390,969	0	A	T	6.3%	199.5 / 3.8	63	S354T (TCC→ACC)	tldD	putative peptidase
Reads supporting (aligned to +/- strand):  ref base A (26/33);  new base T (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

GCCGGAGGTGATATCCACCTGACCGCCACCAAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGCCACCGAACCTCGGC  >  NC_000913/3390826‑3391115                                                                                                                                                |                                                                                                                                                   gCCGGAGGTGATATCCACCTGACCGCCACCAAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCa                                                                                                                                                <  1:41492/148‑1 (MQ=255)   cGGAGGTGATATCCACCTGACCGCCACCAAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACgg                                                                                                                                             <  7:300043/149‑1 (MQ=255)   cGGAGGTGATATCCACCTGACCGCCACCAAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACgg                                                                                                                                             <  8:55343/149‑1 (MQ=255)           atatCCACCTGACCGCCACCAAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGtt                                                                                                                                     <  7:294410/149‑1 (MQ=255)            tatCCACCTGACCGCCACCAAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTg                                                                                                                                    >  8:5889/1‑149 (MQ=255)                 aCCTGACCGCCACCAAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGt                                                                                                                               <  7:2036/149‑1 (MQ=255)                    tGACCGCCACCAAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCgg                                                                                                                            >  6:330091/1‑149 (MQ=255)                       ccGCCACCAAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGt                                                                                                                         >  8:44574/1‑149 (MQ=255)                        cGCCACCAAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACg                                                                                                                                              >  4:39983/1‑127 (MQ=255)                           caccaAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATc                                                                                                                     >  3:67452/1‑149 (MQ=255)                              caAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATcccc                                                                                                                  <  5:297685/149‑1 (MQ=255)                               aaaGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCa                                                                                                                 <  5:330091/149‑1 (MQ=255)                                aaGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCAt                                                                                                                <  2:210579/149‑1 (MQ=255)                                      ggcgCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAAc                                                                                                           <  7:168460/145‑1 (MQ=255)                                         gCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAAcgcgc                                                                                                       <  6:267585/149‑1 (MQ=255)                                           aTAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGCCCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGt                                                                                                     >  6:116383/1‑149 (MQ=255)                                           aTAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGt                                                                                                     >  2:179197/1‑149 (MQ=255)                                           aTAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGt                                                                                                     <  8:269827/149‑1 (MQ=255)                                                  ccGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGttt                                                                                              <  2:322040/149‑1 (MQ=255)                                                   cGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTg                                                                                                                                    <  2:46310/110‑1 (MQ=255)                                                   cGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTg                                                                                                                                    >  1:46310/1‑110 (MQ=255)                                                     tACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGAGGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATc                                                                                           <  4:341432/149‑1 (MQ=255)                                                       cTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATcc                                                                                          <  5:21485/148‑1 (MQ=255)                                                        tCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTg                                                                                        >  8:328782/1‑149 (MQ=255)                                                          aaCGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGCCCGTTGCCAGTCGGCGGCATCCCCATCAAACGCCCGTTGAGGTTATCCTGCa                                                                                      <  4:270929/149‑1 (MQ=255)                                                          aaCGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGc                                                                                       <  4:67452/148‑1 (MQ=255)                                                            cGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATAAGCGGCATGTGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATg                                                                                    >  7:166148/1‑149 (MQ=255)                                                             ggATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATg                                                                                    >  2:134950/1‑148 (MQ=255)                                                               aTTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAg                                                                                 <  1:13628/149‑1 (MQ=255)                                                                   aataatTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGa                                                                                                  <  2:267373/128‑1 (MQ=255)                                                                   aataatTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGa                                                                                                  >  1:267373/1‑128 (MQ=255)                                                                   aataatTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCtt                                                                             >  3:207759/1‑149 (MQ=255)                                                                   aataatTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCtt                                                                             >  2:172852/1‑149 (MQ=255)                                                                     taatTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTc                                                                           <  7:5889/149‑1 (MQ=255)                                                                            cTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTACCAGTCGGCGTCATCCCCATCAAAcgcgcg                                                                                                      >  7:213709/1‑115 (MQ=255)                                                                            cTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTACCAGTCGGCGTCATCCCCATCAAAcgcgcg                                                                                                      <  8:213709/115‑1 (MQ=255)                                                                                      gATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAAt                                                          <  1:172852/149‑1 (MQ=255)                                                                                       aTTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCa                                                                                      >  4:17224/1‑120 (MQ=255)                                                                                               ggCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGtt                                                 <  4:207759/149‑1 (MQ=255)                                                                                                gcagcaTATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTg                                                <  1:179197/149‑1 (MQ=255)                                                                                                    caTATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTAct                                            <  4:77817/149‑1 (MQ=255)                                                                                                          ggTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTActggcctg                                      <  1:134950/149‑1 (MQ=255)                                                                                                            tgtTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTActggcctggc                                    <  7:328782/149‑1 (MQ=255)                                                                                                                 gTCATACGCGGGATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGGCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTAc                                <  8:166148/148‑1 (MQ=255)                                                                                                                    aTACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTc                            <  5:116383/149‑1 (MQ=255)                                                                                                                    aTACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTc                            <  7:44574/149‑1 (MQ=255)                                                                                                                       cgcgGCATGGGCAGATGGGCGTAGGTTTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCa                                                                                      >  4:253072/1‑88 (MQ=255)                                                                                                                       cgcgGCATGGGCAGATGGGCGTAGGTTTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCa                                                                                      >  4:248464/1‑88 (MQ=255)                                                                                                                       cgcgGCATGGGCAGATGGGCGTAGGTTTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCa                                                                                      <  3:248464/88‑1 (MQ=255)                                                                                                                       cgcgGCATGGGCAGATGGGCGTAGGTTTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCa                                                                                      <  3:253072/88‑1 (MQ=255)                                                                                                                        gcgGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCTGCGTCATCCCCATCAACCGCGCGGTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGGTCTCAATCAGCCCGTTGTACTGGCCGGGCGTACCTTCGtca                        >  2:96695/1‑149 (MQ=255)                                                                                                                        gcgGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGACTGCCGTGCTCAATCAGCACGTTGTACTGGCCTGGGGTCCCGTCGtca                        >  2:91670/1‑149 (MQ=255)                                                                                                                         cgGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGtcat                       <  5:355260/149‑1 (MQ=255)                                                                                                                                ggCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGc                >  8:69193/1‑149 (MQ=255)                                                                                                                                ggCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGc                >  8:62804/1‑149 (MQ=255)                                                                                                                                ggCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGc                <  8:226372/149‑1 (MQ=255)                                                                                                                                 gCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCACACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGc                >  1:171245/1‑148 (MQ=255)                                                                                                                                 gCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGc                >  5:108828/1‑148 (MQ=255)                                                                                                                                            tAGGATTCACGGCGCCCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGCCACCGAACCTCg    >  7:335505/1‑149 (MQ=255)                                                                                                                                             aGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGCCACCGAACCTCgg   >  3:73494/1‑149 (MQ=255)                                                                                                                                              ggATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGCCACCGAACCTCgg   <  7:69193/148‑1 (MQ=255)                                                                                                                                               gATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGCCACCGAACCTCGGc  >  1:70116/1‑148 (MQ=255)                                                                                                                                               gATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGCCACCGAACCTCGGc  >  6:357085/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                   GCCGGAGGTGATATCCACCTGACCGCCACCAAAGTTCGGTGCATAGATACCGTACTCAACGGATTCAATAATTTCCTGCGGGGTCGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGCCACCGAACCTCGGC  >  NC_000913/3390826‑3391115

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