Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,721,265	C→A	P23Q (CCG→CAG)	insK →	IS150 conserved protein InsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,721,265	0	C	A	100.0%	137.0 / NA	40	P23Q (CCG→CAG)	insK	IS150 conserved protein InsB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (20/20);  total (20/20)

AGTTCGATTTCTTGAAACGCGGCTTATGTATCTAAAAAAGCTGAAAGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACCGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAGGGTAACGCTGTCTCTTC  >  NC_000913/3721129‑3721399                                                                                                                                         |                                                                                                                                       aGTTCGATTTCTTGAAACGCGGCTTATGTATCTAAAAAAGCTGAAAGCCTTAGCTCATCCCACGAAAAGGTGAAAGAACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACagcg                                                                                                                                     <  2:737859/140‑1 (MQ=255) aGTTCGATTTCTTGAAACGCGGCTTATGTATCTAAAAAAGCTGAAAGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACagcg                                                                                                                                     <  2:481537/140‑1 (MQ=255)     cGATTTCTTGAAACGCGGCTTATGTATCTAAAAAAGCTGAAAGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGt                                                                                                                                 <  1:666753/140‑1 (MQ=255)             tGAAACGCGGCTTATGTATCTAAAAAAGCTGAAAGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTtta                                                                                                                         >  1:482779/1‑140 (MQ=255)                   gcgGCTTATGTATCTAAAAAAGCTGAAAGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTatca                                                                                                                   <  2:482779/140‑1 (MQ=255)                                 tAAAAAAGCTGAAAGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCa                                                                                                     >  2:830796/1‑140 (MQ=255)                                   aaaaaGCTGAAAGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGc                                                                                                   <  2:758194/140‑1 (MQ=255)                                      aaGCTGAAAGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAg                                                                                                <  2:247234/140‑1 (MQ=255)                                         cTGATAGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCt                                                                                             >  1:282500/1‑140 (MQ=255)                                           gAAAGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGa                                                                                           >  2:2992/1‑140 (MQ=255)                                            aaaGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGAc                                                                                          <  2:222953/140‑1 (MQ=255)                                                          tCCCACGAAAAAGTGAAGGTACTCAACGAGCGAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGt                                                                            >  2:711608/1‑140 (MQ=255)                                                          tCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGt                                                                            >  1:108682/1‑140 (MQ=255)                                                               cGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACATGTATGCGGACGTTaaaa                                                                       >  1:450861/1‑140 (MQ=255)                                                                   aaaGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCg                                                                   >  1:134090/1‑140 (MQ=255)                                                                   aaaGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCg                                                                   >  1:670593/1‑140 (MQ=255)                                                                          aaGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGt                                                            <  1:906385/140‑1 (MQ=255)                                                                            gTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTga                                                          <  1:282672/140‑1 (MQ=255)                                                                                 cAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGAttt                                                     <  2:109686/140‑1 (MQ=255)                                                                                   aCGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTAt                                                   <  2:450232/140‑1 (MQ=255)                                                                                   aCGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTAt                                                   <  1:562295/140‑1 (MQ=255)                                                                                    cGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATc                                                  >  1:460652/1‑140 (MQ=255)                                                                                     gAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACATGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCa                                                 <  2:450861/140‑1 (MQ=255)                                                                                                 ttttATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAgg                                     >  1:611402/1‑140 (MQ=255)                                                                                                     aTCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGa                                 <  2:282500/140‑1 (MQ=255)                                                                                                        ctctTGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATAc                              >  2:739840/1‑140 (MQ=255)                                                                                                           ttGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGa                           <  1:830796/140‑1 (MQ=255)                                                                                                           ttGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGa                           <  1:2992/140‑1 (MQ=255)                                                                                                           ttGATGAGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGa                           >  1:686347/1‑140 (MQ=255)                                                                                                             gatgaGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATa                         >  1:594116/1‑140 (MQ=255)                                                                                                                gaGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCg                      <  1:739840/140‑1 (MQ=255)                                                                                                                gaGCTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCg                      <  2:686347/140‑1 (MQ=255)                                                                                                                   cTTCTCAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAg                   >  1:677919/1‑140 (MQ=255)                                                                                                                      ctcAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAGGGt                <  2:670593/140‑1 (MQ=255)                                                                                                                       tcAGGGCTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAGGGTa               >  2:477458/1‑140 (MQ=255)                                                                                                                             cTGCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAGGGTAACGCTg         >  1:255875/1‑140 (MQ=255)                                                                                                                               gCGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAGGGTAACGCTGtc       <  1:219517/140‑1 (MQ=255)                                                                                                                                cGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAGGGTAACGCTGtct      >  1:545095/1‑140 (MQ=255)                                                                                                                                cGGAGATACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAGGGTAACGCTGtct      >  2:207859/1‑140 (MQ=255)                                                                                                                                    gaTACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAGGGTAACGCTGTCTCTTc  <  2:493781/140‑1 (MQ=255)                                                                                                                                     aTACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAGGGTAACGCTGTCTCtt   >  1:429439/1‑138 (MQ=255)                                                                                                                                     aTACAGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAGGGTAACGCTGTCTCtt   <  2:429439/138‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       AGTTCGATTTCTTGAAACGCGGCTTATGTATCTAAAAAAGCTGAAAGCCTTAGCTCATCCCACGAAAAAGTGAAAGTACTCAACGAGCTAAGGCAGTTTTATCCTCTTGATGAGCTTCTCAGGGCTGCGGAGATACCGCGCAGTACGTTTTATTATCATCTAAAGGCTCTCAGCAAGCCTGACAAGTATGCGGACGTTAAAAAGCGTATTAGTGAGATTTATCACGAGAATAGAGGCCGATACGGATACCGTAGGGTAACGCTGTCTCTTC  >  NC_000913/3721129‑3721399

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: ‑