Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,159,354	C→T	100%	R323R (CGC→CGT)	mdtC →	multidrug efflux system, subunit C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,159,354	0	C	T	98.9%	285.2 / ‑4.4	91	R323R (CGC→CGT)	mdtC	multidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (37/52);  total (37/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

CCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTC  >  NC_000913/2159257‑2159449                                                                                                  |                                                                                                ccAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTc                                                                                               >  2:559346/1‑100 (MQ=255)  cAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTcc                                                                                              <  2:1389136/100‑1 (MQ=255)    atatTATCCAGACGGTTGACAGTATCCGGTTAAAATTACCGGAGTTGTAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTcccc                                                                                            <  1:2646315/100‑1 (MQ=255)       ttatCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCccac                                                                                         <  2:2594688/100‑1 (MQ=255)       ttatCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTACCCAGGATCGTTCCCccac                                                                                         >  1:522576/1‑100 (MQ=255)        tatCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCccacc                                                                                        <  1:1678777/100‑1 (MQ=255)        tatCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCccac                                                                                         <  2:1097761/99‑1 (MQ=255)         atCCAGACGGTTGACAGCATCCAGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGCCGATCGATCTGCAAATTGCCCAATATCGTTCCCccac                                                                                         <  2:698096/98‑1 (MQ=255)          tCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCAt                                                                                      >  2:767892/1‑100 (MQ=255)           ccAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCAtt                                                                                     >  2:1208006/1‑100 (MQ=255)           ccAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCAtt                                                                                     >  2:2454535/1‑100 (MQ=255)            cAGACGGTTGACAGTATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTc                                                                                    <  1:1470509/100‑1 (MQ=255)             aGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTc                                                                                    <  2:383974/99‑1 (MQ=255)             aGACGGTTGACAGCAGCCGGGCAAGATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCccacca                                                                                       <  2:2638958/96‑1 (MQ=255)               aCGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTcgcg                                                                                 >  2:2401292/1‑100 (MQ=255)                cGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCgcgc                                                                                <  2:148482/100‑1 (MQ=255)                   ttGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTc                                                                             <  1:1144463/100‑1 (MQ=255)                    tGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCg                                                                            <  2:2354023/100‑1 (MQ=255)                     gACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGc                                                                           >  1:1105658/1‑100 (MQ=255)                       cAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTg                                                                         <  2:468434/100‑1 (MQ=255)                       cAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTg                                                                         <  2:1683033/100‑1 (MQ=255)                       cAGCATCCGGGCAAAAGTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTg                                                                         <  2:293995/100‑1 (MQ=255)                        aGCATCCGGGCAAAATTACCGGTGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTg                                                                         <  1:1399646/99‑1 (MQ=255)                         gCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGga                                                                       >  2:806282/1‑100 (MQ=255)                           aTCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGgaag                                                                     <  1:1368920/100‑1 (MQ=255)                           aTCCGGGCAAAATTACCGGAGTTGCAGGAAACCATACCGGCCGCGATAGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTAGCGCCTCGCTGgaag                                                                     <  2:1046784/100‑1 (MQ=255)                              cGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGaagaag                                                                  >  1:1450599/1‑100 (MQ=255)                              cGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGaagaag                                                                  >  2:2091748/1‑100 (MQ=255)                              cGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCACAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGaagaag                                                                  <  2:149511/100‑1 (MQ=255)                               gggCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGgaaga                                                                    <  2:738794/97‑1 (MQ=255)                                ggCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTc                                                                >  2:280784/1‑100 (MQ=255)                                 gCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCg                                                               <  1:1619443/100‑1 (MQ=255)                                 gCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCg                                                               >  1:564439/1‑100 (MQ=255)                                  cAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGa                                                              >  2:761454/1‑100 (MQ=255)                                  cAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGa                                                              <  1:559349/100‑1 (MQ=255)                                   aaaaTTACCGTAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAg                                                             <  1:1264661/100‑1 (MQ=255)                                     aaTTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCa                                                           >  1:2235166/1‑100 (MQ=255)                                        tACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAAc                                                        <  1:49713/100‑1 (MQ=255)                                           cGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCt                                                     <  2:1010588/100‑1 (MQ=255)                                            ggAGTTGCAGGAAACCATTCCGGCGGCGATTGCTCTGCAGATCGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTg                                                    <  2:2127633/100‑1 (MQ=255)                                            ggAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTg                                                    >  2:2804536/1‑100 (MQ=255)                                              aGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGAt                                                  <  2:368613/100‑1 (MQ=255)                                              aGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGAt                                                  <  2:1559543/100‑1 (MQ=255)                                               gTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGAtt                                                 >  1:1976056/1‑100 (MQ=255)                                               gTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGAGCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGAtt                                                 <  2:866576/100‑1 (MQ=255)                                                ttGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTa                                                >  1:152739/1‑100 (MQ=255)                                                ttGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCAATCGAGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTa                                                >  1:158515/1‑100 (MQ=255)                                                ttGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTCCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTa                                                <  2:2967287/100‑1 (MQ=255)                                                 tGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAt                                               >  1:1898227/1‑100 (MQ=255)                                                 tGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAt                                               >  1:84071/1‑100 (MQ=255)                                                 tGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAt                                               >  2:280112/1‑100 (MQ=255)                                                  gCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAtc                                              <  2:704792/100‑1 (MQ=255)                                                  gCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAtc                                              >  2:1986904/1‑100 (MQ=255)                                                   cAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAtct                                             >  1:2209166/1‑100 (MQ=255)                                                     ggAAACCATTCCGGCGGCGTTTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCg                                           >  1:2486735/1‑100 (MQ=255)                                                     ggAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCg                                           >  2:1705889/1‑100 (MQ=255)                                                        aaCCATTCCGGCGGCGATTGATCTGCAAATTGCCTAGGATCGTTCCCGCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTg                                        <  2:1887527/100‑1 (MQ=255)                                                        aaCCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTg                                        <  1:1663022/100‑1 (MQ=255)                                                        aaCCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTg                                        <  1:100059/100‑1 (MQ=255)                                                         aCCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAATCGCTGATTATCTCGGTgg                                       <  2:1975327/100‑1 (MQ=255)                                                         aCCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTgg                                       <  1:3048218/100‑1 (MQ=255)                                                         aCCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTgg                                       <  1:1174741/100‑1 (MQ=255)                                                         aCCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTgg                                       <  1:118158/100‑1 (MQ=255)                                                         aCCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCCCCTCGCTGGATGAAGTCGAGCAAACGCTTATTATCTCGGttg                                       <  1:2047676/100‑3 (MQ=255)                                                          ccATTCCGGCGGCGATTGATCTGCAAATTGCTCAGGATCGTTCTCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTgg                                       <  1:280115/99‑1 (MQ=255)                                                            aTTCCGGCGGCGATTGATCTGCACATTGCCCAGGATCGTTCCCCCCCCATTCGCGCCTCGCTGGAAGAAGTAGAACAAACCCTAATTATCTCTGTGgcg                                     >  1:2903937/1‑99 (MQ=255)                                                            aTTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGgcgc                                    >  1:2701714/1‑100 (MQ=255)                                                             ttCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCt                                   <  1:1229593/100‑1 (MQ=255)                                                             ttCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCt                                   <  1:2207715/100‑1 (MQ=255)                                                               ccGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTgg                                 <  1:1165821/100‑1 (MQ=255)                                                                cggcggCGATTGATCTGCAAATTGCCCAGGATCGTTTCCCCATCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGt                                <  1:2019040/100‑1 (MQ=255)                                                                  gcgacgATTGATCTCCAAATCGGCCAGGATTTTATCTCGACCTTTCGCGCCTCGCTGCAAGAAGTCGAGCCAACGCTGATTATCTCGGTGGCGCTgg                                 <  2:851201/97‑1 (MQ=255)                                                                       gATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCTTTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCtg                         <  1:2869026/100‑1 (MQ=255)                                                                          ttgtttGCAGATTGCCCAGGCTCGGTCTCTCTTCTTTCGTGTCTCGCTGGTAGTAGTCTAGCTATCGCTGATTATCTCGGTGGCGCTGGTGATTc                           <  2:1038785/90‑1 (MQ=255)                                                                          tGATCTGCAGAGTGGCCGGGATCCTTCCTCCGCCATTGGCGCCTTGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCAGGTGATTCtggtn                      <  2:2284206/100‑2 (MQ=255)                                                                          tGATCTGCAAATTGCCAAGGATCGTTGCCCCTCCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCtggtg                      <  2:888054/100‑1 (MQ=255)                                                                           gATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGTCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCtggtgg                     >  1:533637/1‑100 (MQ=255)                                                                           gATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCtggtgg                     >  2:2669269/1‑100 (MQ=255)                                                                             tCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGTAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGgtggtg                   <  2:270032/100‑1 (MQ=255)                                                                             tCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCCGGTGATTCTGgtggtg                   >  2:384421/1‑100 (MQ=255)                                                                               tGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGtt                 >  1:1569298/1‑100 (MQ=255)                                                                                gCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATATCGATAGCGCTGGTGATTCTGGTGGTGttt                >  1:523776/1‑100 (MQ=255)                                                                                  aaaTTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGttttt              <  1:2775812/100‑1 (MQ=255)                                                                                     ttGCCTAGGTTCGTTCCCCCACTATTAGCGCCTTGTTGGGAGAAGTCGAGCAAACTCTGATTATCTCGGTTGCGGTGGTGATTCTGGTGGTGTTTTTAtt           <  1:1326804/100‑1 (MQ=255)                                                                                     ttGCCCAGGATCGTTCCCCCACCATTCGTGCCTCTCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTAtt           <  1:2401320/100‑1 (MQ=255)                                                                                       gCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAATAAGTCGAGTAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTcc         <  1:1949913/100‑1 (MQ=255)                                                                                       gCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGAGTTTTTATTcc         >  2:2413436/1‑100 (MQ=255)                                                                                         ccAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATATCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTg       >  1:1149789/1‑100 (MQ=255)                                                                                          cAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGGAGGTGATTCTGGAGGTGATTTTATTCCTgc      >  2:1075041/1‑100 (MQ=255)                                                                                            ggATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTAATTATCTCGGTGGCGCAGGTGATTCTGGTGGTGTTTTTATTCCTgcgc    >  1:223122/1‑100 (MQ=255)                                                                                              aTCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCCGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTc  >  2:2203960/1‑100 (MQ=255)                                                                                                  |                                                                                                CCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTC  >  NC_000913/2159257‑2159449

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