Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,583,506	T→C	C8C (TGT→TGC)	insA‑6 →	IS1 protein InsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,583,506	0	T	C	100.0%	77.6 / NA	24	C8C (TGT→TGC)	insA‑6	IS1 protein InsA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/7);  total (17/7)

CAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATC  >  NC_000913/3583371‑3583564                                                                                                                                        |                                                           cAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCcct                                                        >  2:59499/1‑140 (MQ=255)        ggTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTc                                                 >  1:663115/1‑140 (MQ=255)           tAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGcccc                                                         >  1:105449/1‑129 (MQ=255)           tAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGcccc                                                         <  2:105449/129‑1 (MQ=255)           tAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGc                                              >  2:619193/1‑140 (MQ=255)           tAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGc                                              >  2:294003/1‑140 (MQ=255)                   ggggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACg                                      <  1:619193/140‑1 (MQ=255)                   ggggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACg                                      >  1:407760/1‑140 (MQ=255)                   ggggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACg                                      >  1:14918/1‑140 (MQ=255)                      gATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACgggg                                   <  1:189254/140‑1 (MQ=255)                          ataGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTAc                                           <  1:111678/128‑1 (MQ=255)                          ataGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTAc                                           >  2:111678/1‑128 (MQ=255)                                 ttCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTATCTATCAGCTGCCCCTCCTGTTCAGCTACTGACgggg                                   >  2:687484/1‑129 (MQ=255)                                 ttCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTATCTATCAGCTGCCCCTCCTGTTCAGCTACTGACgggg                                   <  1:687484/129‑1 (MQ=255)                                   cAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACgg                      >  1:423443/1‑140 (MQ=255)                                   cAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGAGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACgg                      >  1:432003/1‑140 (MQ=255)                                     atatTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCa                    >  1:104686/1‑140 (MQ=255)                                           aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAAc                        >  2:479110/1‑130 (MQ=255)                                           aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAAc                        <  1:479110/130‑1 (MQ=255)                                            cTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAc             <  1:294003/140‑1 (MQ=255)                                            cTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAc             <  1:319007/140‑1 (MQ=255)                                                tGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAccgcc         >  1:671349/1‑140 (MQ=37)                                                tGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAccgcc         >  1:721159/1‑140 (MQ=37)                                                tGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAccgcc         >  1:590231/1‑140 (MQ=37)                                                      tAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACAt   >  1:272581/1‑140 (MQ=17)                                                       aaaGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGCCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATc  >  2:62194/1‑140 (MQ=14)                                                                                                                                        |                                                           CAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATC  >  NC_000913/3583371‑3583564

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