Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,584,131	T→A	L152Q (CTG→CAG)	insB‑6 →	IS1 protein InsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,584,131	0	T	A	100.0%	167.1 / NA	49	L152Q (CTG→CAG)	insB‑6	IS1 protein InsB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (21/28);  total (21/28)

ACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGAT  >  NC_000913/3584062‑3584267                                                                      |                                                                                                                                         aCCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAtt                                                                    <  2:720267/140‑1 (MQ=14)  ccTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTa                                                                   >  2:916443/1‑140 (MQ=17)      aaTCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                                                               >  1:1381817/1‑140 (MQ=25)       aTCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCACTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGt                                                              >  1:238121/1‑140 (MQ=255)         cTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                                                            <  1:657806/140‑1 (MQ=38)           gAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACt                                                          >  1:363014/1‑140 (MQ=255)            aGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtt                                                         >  2:981976/1‑140 (MQ=255)            aGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtt                                                         >  2:1256492/1‑140 (MQ=255)            aGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtt                                                         >  1:1532333/1‑140 (MQ=255)            aGGCAGCACCTGGCACGGCTGGGACGAAAGTCGCTGTCGTTCTCAAAATCGGTGGATCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtt                                                         >  2:1404620/1‑140 (MQ=255)               cagcaCCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                                                      <  1:1104370/140‑1 (MQ=255)                agcaCCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGATTCATTACCGACATTAAAAGTAACtttttt                                                     <  2:509454/140‑1 (MQ=255)                  caCCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa                                                   >  1:853705/1‑140 (MQ=255)                   aCCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaa                                                  >  1:795542/1‑140 (MQ=255)                    ccTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa                                                 <  1:1119930/140‑1 (MQ=255)                     cTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAt                                                <  1:1579498/140‑1 (MQ=255)                     cTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAt                                                <  1:916443/140‑1 (MQ=255)                      tGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                               <  1:685109/140‑1 (MQ=255)                      tGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                               <  2:464071/140‑1 (MQ=255)                      tGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                               <  2:841223/140‑1 (MQ=255)                       ggCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa                                              <  2:1376306/140‑1 (MQ=255)                       ggCACGGCTGGGACGGAAGTCACTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa                                              <  2:238121/140‑1 (MQ=255)                          aCGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag                                           <  2:1062155/140‑1 (MQ=255)                            ggCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat                                         <  2:795542/140‑1 (MQ=255)                            ggCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat                                         <  1:1256492/140‑1 (MQ=255)                            ggCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat                                         >  1:933114/1‑140 (MQ=255)                            ggCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat                                         <  2:363014/140‑1 (MQ=255)                              cTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa                                       >  2:898203/1‑140 (MQ=255)                                gggACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataata                                     <  1:47213/140‑1 (MQ=255)                                   aCGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag                                           >  2:1000782/1‑131 (MQ=255)                                   aCGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag                                           <  1:1000782/131‑1 (MQ=255)                                     ggAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc                                <  2:1532333/140‑1 (MQ=255)                                      gAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                               >  2:1378363/1‑124 (MQ=255)                                      gAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                               <  1:1378363/124‑1 (MQ=255)                                      gAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGca                               <  2:45696/140‑1 (MQ=255)                                       aaGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGcac                              >  1:862374/1‑140 (MQ=255)                                       aaGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGcac                              >  1:103122/1‑140 (MQ=255)                                       aaGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGcac                              <  2:388961/140‑1 (MQ=255)                                         gTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAg                            <  1:1338740/140‑1 (MQ=255)                                         gTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAg                            <  1:560828/140‑1 (MQ=255)                                         gTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAg                            <  2:708365/140‑1 (MQ=255)                                         gTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAg                            <  2:51824/140‑1 (MQ=255)                                            gCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGaaa                         >  2:1372058/1‑140 (MQ=255)                                             cTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGaaaa                        <  2:268056/140‑1 (MQ=255)                                                   ttCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATAtt                  >  2:619121/1‑140 (MQ=255)                                                   ttCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATAtt                  >  2:827529/1‑140 (MQ=255)                                                       cAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGt              <  1:981976/140‑1 (MQ=255)                                                              ggTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTa       >  1:866629/1‑140 (MQ=255)                                                                tGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGt     >  1:1014416/1‑140 (MQ=255)                                                                tGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGt     >  1:638548/1‑140 (MQ=255)                                                                   agcagcATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGAt  <  2:1381817/140‑1 (MQ=255)                                                                      |                                                                                                                                         ACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGAT  >  NC_000913/3584062‑3584267

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