Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	291,343	C→T	A2T (GCT→ACT)	insA‑3 ←	IS1 protein InsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	291,343	0	C	T	100.0%	142.1 / NA	45	A2T (GCT→ACT)	insA‑3	IS1 protein InsA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/25);  total (20/25)

GATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCA  >  NC_000913/291265‑291478                                                                               |                                                                                                                                        gATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGaa                                                                            >  1:556469/1‑140 (MQ=11)  aTGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAAt                                                                           >  1:24835/1‑140 (MQ=14)            gTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCg                                                                 >  1:429285/1‑140 (MQ=255)            gTGCTTTTGCCATTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCg                                                                 >  1:603360/1‑140 (MQ=255)                ttttGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATg                                                             <  2:125889/140‑1 (MQ=255)                 tttGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGt                                                            <  2:66384/140‑1 (MQ=255)                  ttGCCGTTACGCACCACCCCGTGAGTAGCTGCACAGGAGGGACAGCTGTTAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGtt                                                           <  2:247830/140‑1 (MQ=255)                  ttGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGtt                                                           >  2:675986/1‑140 (MQ=255)                    gCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTc                                                         <  1:578037/140‑1 (MQ=255)                     ctgTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACTCCATCATACACTAAATCAGTAAGTTGGCAGCATCGCCGAATTTTCAGCCCGGATGTTTCa                                                        <  1:142970/138‑1 (MQ=255)                     ccGTTACGCACCACCTCGTCAGCAGCTGATCAGGAGGTACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCa                                                        <  1:242160/140‑1 (MQ=14)                       gTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACg                                                      >  2:514755/1‑140 (MQ=255)                         tACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCACGGATGTTTCACGTa                                                    >  1:397721/1‑140 (MQ=255)                          aCGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTAc                                                   <  2:707534/140‑1 (MQ=255)                            gCACCACCGCGTCAGTAGCTGAACAGGAGGGTCCGTTGAGAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCt                                                 <  2:397721/140‑1 (MQ=17)                             caccacCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTa                                                <  1:469140/140‑1 (MQ=255)                             caccacCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTa                                                <  2:24835/140‑1 (MQ=255)                                caccCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGt                                             >  2:99661/1‑140 (MQ=255)                                     cGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAat                                        >  2:648040/1‑140 (MQ=255)                                          gTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTc                                            <  2:190834/131‑1 (MQ=255)                                          gTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTc                                            >  1:190834/1‑131 (MQ=255)                                          gTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTg                                   <  2:429285/140‑1 (MQ=255)                                             gCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGcc                                >  1:110780/1‑140 (MQ=255)                                             gCTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGcc                                >  2:790650/1‑140 (MQ=255)                                              cTGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACAGAATTTTCGGCCGGGATGTTTCACGAACCTATGTCGAATATCTGGCCg                               >  2:375316/1‑140 (MQ=18)                                               tGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGa                              <  2:110780/140‑1 (MQ=255)                                               tGAACAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGa                              <  1:648040/140‑1 (MQ=255)                                                 acagggaggGACAGCGGATAGAGACAGAAGTCACTGGTGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATa                            <  1:421134/136‑1 (MQ=21)                                                 aaCAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATa                            >  2:496299/1‑140 (MQ=255)                                                 aaCAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATa                            >  2:827739/1‑140 (MQ=255)                                                 aaCAGGAGGGACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATa                            <  1:514755/140‑1 (MQ=255)                                                        gggACAGCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGc                     <  2:603360/140‑1 (MQ=255)                                                              gCTGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGc                     <  2:853869/134‑1 (MQ=255)                                                              gCTGATAGAAACAGAAGTCACTGAAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGGTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGc                     >  1:853869/1‑134 (MQ=255)                                                                tGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGa                              <  2:513917/123‑1 (MQ=255)                                                                tGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGa                              >  1:513917/1‑123 (MQ=255)                                                                tGATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATATTTCACGTACCTGTGTGTAATATCTCGCCGCCAAGGTTGAGGATctgt             >  2:421134/1‑137 (MQ=18)                                                                 gATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATa                            <  2:745601/124‑1 (MQ=255)                                                                 gATAGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATa                            >  1:745601/1‑124 (MQ=255)                                                                    aGAAACAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATATCGAATACCTGGCCGATAAGGGTGCGGATCGagaaga         >  1:374230/1‑140 (MQ=255)                                                                       aaCAGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGa      <  1:790650/140‑1 (MQ=255)                                                                          aGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGc   <  1:491769/140‑1 (MQ=255)                                                                          aGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGc   <  2:556469/140‑1 (MQ=255)                                                                          aGAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGc   <  1:186900/140‑1 (MQ=255)                                                                           gAAGTCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCa  <  1:675986/140‑1 (MQ=255)                                                                               |                                                                                                                                        GATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGAGAAGAGGATGCA  >  NC_000913/291265‑291478

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