Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,661,834	T→G	L474R (CTG→CGG)	mdtF →	multidrug efflux pump RND permease MdtF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,661,834	0	T	G	100.0%	232.3 / NA	64	L474R (CTG→CGG)	mdtF	multidrug efflux pump RND permease MdtF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (30/34);  total (30/34)

GGAAGCGACGCATAAATCGATGGGGCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCTGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCACGCTTCAACACGCTG  >  NC_000913/3661697‑3661964                                                                                                                                          |                                                                                                                                   ggaagCGACGCATAAATCGATGGGGCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTGCTTTCCTCATGGTTTAGGCCGTTGGCCTTTATGAGCGGTGCACCCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGa                                                                                                                                  <  2:134103/140‑1 (MQ=255)  gaagCGAGGCATAAATCGATGGGGCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGAGGGCTTTTATGAGGGGTGCAACCGGGGAGATTTACCGCCAGTTCTCCATCACGCGGAt                                                                                                                                 <  2:1208886/140‑1 (MQ=255)   aagCGACGCATAAATCGATGGGGCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGAtc                                                                                                                                >  1:144366/1‑140 (MQ=255)    agCGACGCATAAATCGATGGGGCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGAtct                                                                                                                               >  2:331847/1‑140 (MQ=255)     gCGACGCATAAATCGATGGGGCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGAtctc                                                                                                                              >  1:300833/1‑140 (MQ=255)            aTAAATCGATGGGGCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCAtg                                                                                                                       <  2:77485/140‑1 (MQ=255)                 tCGATGGGGCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATgctgct                                                                                                                  <  2:854303/140‑1 (MQ=255)                 tCGATGGGGCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATgctgct                                                                                                                  <  1:941840/140‑1 (MQ=255)                 tCGATGGGGCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATgctgct                                                                                                                  <  1:45693/140‑1 (MQ=255)                       gggCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGt                                                                                                            <  2:567971/140‑1 (MQ=255)                       gggCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGt                                                                                                            <  1:1185969/140‑1 (MQ=255)                       gggCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGt                                                                                                            <  1:280174/140‑1 (MQ=255)                          cAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTAtt                                                                                                         >  2:692711/1‑140 (MQ=255)                           aGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTAttt                                                                                                        >  1:987240/1‑140 (MQ=255)                                  aCGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCaa                                                                                                 <  1:571981/140‑1 (MQ=255)                                  aCGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCaa                                                                                                 >  2:174574/1‑140 (MQ=255)                                         cTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCt                                                                                          <  1:454753/140‑1 (MQ=255)                                         cTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCt                                                                                          <  1:1201065/140‑1 (MQ=255)                                         cTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCt                                                                                          <  2:631393/140‑1 (MQ=255)                                         cTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCt                                                                                          <  1:837219/140‑1 (MQ=255)                                           ggTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGa                                                                                        <  2:760443/140‑1 (MQ=255)                                            gTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCTATGAGCCTGAc                                                                                       >  2:617335/1‑140 (MQ=255)                                                 tATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACccctg                                                                                  >  1:294398/1‑140 (MQ=255)                                                 tATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACccctg                                                                                  <  2:477558/140‑1 (MQ=255)                                                    tGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCt                                                                                          >  1:1189273/1‑129 (MQ=255)                                                    tGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCt                                                                                          <  2:1189273/129‑1 (MQ=255)                                                      ccGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACccctgccct                                                                              >  1:1213296/1‑139 (MQ=255)                                                      ccGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACccctgccct                                                                              <  2:1213296/139‑1 (MQ=255)                                                       cGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCtgt                                                                            >  2:1110192/1‑140 (MQ=255)                                                                      aGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGa                                                             <  2:1038772/140‑1 (MQ=255)                                                                      aGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGAGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGa                                                             >  2:1120066/1‑140 (MQ=255)                                                                         gtTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAg                                                          <  1:858747/140‑1 (MQ=255)                                                                             aTGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCcgc                                                      >  2:7442/1‑140 (MQ=255)                                                                              tGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCcgcg                                                     <  1:695988/140‑1 (MQ=255)                                                                               gCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCgcgc                                                    <  1:597646/140‑1 (MQ=255)                                                                                  gATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCgg                                                 >  1:142589/1‑140 (MQ=255)                                                                                  gATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCgg                                                 >  1:1142765/1‑140 (MQ=255)                                                                                    tGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGaa                                               >  1:952506/1‑140 (MQ=255)                                                                                    tGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGaa                                               >  1:1079965/1‑140 (MQ=255)                                                                                    tGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGaa                                               >  1:194011/1‑140 (MQ=255)                                                                                       ccTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGc                                            >  1:201261/1‑140 (MQ=255)                                                                                       ccTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGc                                            <  2:266672/140‑1 (MQ=255)                                                                                               aGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACaa                                    <  1:855383/140‑1 (MQ=255)                                                                                               aGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACaa                                    <  1:1017864/140‑1 (MQ=255)                                                                                               aGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACaa                                    <  2:144366/140‑1 (MQ=255)                                                                                               aGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACaa                                    <  2:294398/140‑1 (MQ=255)                                                                                                gCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACaaa                                   >  1:280612/1‑140 (MQ=255)                                                                                                   gTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCt                                >  1:1008904/1‑140 (MQ=255)                                                                                                    tGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTa                               >  2:465566/1‑140 (MQ=255)                                                                                                    tGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTa                               >  2:581182/1‑140 (MQ=255)                                                                                                      cAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAAc                             >  1:801029/1‑140 (MQ=255)                                                                                                         ccGGGGAGATCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGcc                          >  1:1140509/1‑140 (MQ=255)                                                                                                                gaTCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCg                   <  2:987240/140‑1 (MQ=255)                                                                                                                  tCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCa                 <  1:331847/140‑1 (MQ=255)                                                                                                                  tCTACCGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCa                 <  1:174574/140‑1 (MQ=255)                                                                                                                       cGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCACg               >  1:896234/1‑137 (MQ=255)                                                                                                                       cGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCACg               <  2:896234/137‑1 (MQ=255)                                                                                                                       cGCCAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCACGCtt            >  2:29506/1‑140 (MQ=255)                                                                                                                         ccAGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCACGCTTCa          >  2:1060689/1‑140 (MQ=255)                                                                                                                           aGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCACGCTTCAac        >  2:33071/1‑140 (MQ=255)                                                                                                                           aGTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCACGCTTCAac        >  2:351478/1‑140 (MQ=255)                                                                                                                            gTTCTCCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCACGCTTCAaca       <  1:1001344/140‑1 (MQ=255)                                                                                                                              tctcCATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCACGCTTCAACACg     <  2:1079965/140‑1 (MQ=255)                                                                                                                                 ccATCACGCGGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCACGCTTCAACACGCTg  <  2:952506/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                   GGAAGCGACGCATAAATCGATGGGGCAGATCCAACGTGCGCTGGTCGGTATTGCCGTTGTTCTTTCCGCAGTGTTTATGCCGATGGCCTTTATGAGCGGTGCAACCGGGGAGATCTACCGCCAGTTCTCCATCACGCTGATCTCCTCCATGCTGCTTTCAGTATTTGTGGCAATGAGCCTGACCCCTGCCCTGTGCGCCACCATTCTGAAAGCCGCGCCGGAAGGCGGTCACAAACCTAACGCCCTGTTCGCACGCTTCAACACGCTG  >  NC_000913/3661697‑3661964

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑