Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,068,952	C→T	*302* (TAG→TAA)	insD1 ←	IS2 transposase TnpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,068,952	0	C	T	100.0%	175.3 / NA	55	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/29);  total (26/29)

TGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCT  >  NC_000913/2068820‑2069075                                                                                                                                     |                                                                                                                            tGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTAGAtt                                                                                                                       >  2:518362/1‑139 (MQ=255) tGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATACCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATAtt                                                                                                                       >  2:489034/1‑139 (MQ=255)      aacaaacaTTGTAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAg                                                                                                                  <  1:320945/139‑1 (MQ=255)      aacaaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAg                                                                                                                  >  1:414297/1‑139 (MQ=255)        caaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGAc                                                                                                                <  1:203726/139‑1 (MQ=255)        caaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGAc                                                                                                                >  1:281005/1‑139 (MQ=255)         aaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACa                                                                                                               <  2:250743/139‑1 (MQ=255)          aacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                              <  2:54001/139‑1 (MQ=255)                gCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                        <  1:384535/139‑1 (MQ=255)                      cTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGa                                                                                                                 >  1:50656/1‑124 (MQ=255)                      cTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGa                                                                                                                 <  2:50656/124‑1 (MQ=255)                         caccaATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAc                                                                                               >  2:530584/1‑139 (MQ=255)                           ccaATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAccc                                                                                             >  1:273066/1‑139 (MQ=255)                                      aTAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGccc                                                                                  >  2:740165/1‑139 (MQ=255)                                       tAGCCGATGATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                                 <  1:489034/139‑1 (MQ=255)                                       tAGCCGATAATTTGCAGGACATCCCCGGATTAGAGATAACAGACATCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGTCCg                                                                                 <  1:518362/139‑1 (MQ=255)                                       tAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                                 <  2:783806/139‑1 (MQ=255)                                       tAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                                 <  2:384027/139‑1 (MQ=255)                                           cGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGc                                                                             >  2:138855/1‑139 (MQ=255)                                             aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATc                                                                                                             <  1:306173/105‑1 (MQ=255)                                             aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATc                                                                                                             >  2:306173/1‑105 (MQ=255)                                             aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCg                                                                           <  2:273066/139‑1 (MQ=255)                                             aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCg                                                                           <  1:305840/139‑1 (MQ=255)                                              tAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGc                                                                          <  1:138855/139‑1 (MQ=255)                                                    gCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGatat                                                                    >  2:184456/1‑139 (MQ=255)                                                           aTCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTg                                                             >  1:95138/1‑139 (MQ=255)                                                             ccccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTta                                                                                                                           <  2:719253/75‑1 (MQ=255)                                                             ccccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTta                                                                                                                           >  1:719253/1‑75 (MQ=255)                                                             ccccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgc                                                           >  1:792634/1‑139 (MQ=255)                                                                     tAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATa                                                   <  2:281005/139‑1 (MQ=255)                                                                       gcgaTAGCAGACAAGCCCTTGGCAGGCCGTAAACAAACCCGCAGGCGTGGCTGAATTTGCCCTTATCTTCCCACACCTCTGTGATCGCTTAAAACATTACAAGCCCGCCGCCTCCGACATTCCTGTTGCTAGCCATCAc                                                 >  1:615467/3‑139 (MQ=255)                                                                       gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                         <  2:324393/83‑1 (MQ=255)                                                                       gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                         >  1:324393/1‑83 (MQ=255)                                                                       gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAAc                                                 >  2:546703/1‑139 (MQ=255)                                                                       gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAAc                                                 >  1:575570/1‑139 (MQ=255)                                                                        agaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGt                                                                                                          <  1:598568/81‑1 (MQ=255)                                                                        agaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGt                                                                                                          >  2:598568/1‑81 (MQ=255)                                                                           tAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAg                                             <  2:414297/139‑1 (MQ=255)                                                                              cAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcgc                                          <  1:530584/139‑1 (MQ=255)                                                                               aGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                         >  1:391544/1‑139 (MQ=255)                                                                                gACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTCGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTACCGTGGCGAGCGATAACCCAGCGCAc                                        <  1:792335/139‑1 (MQ=255)                                                                                 aCAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCGGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGTCGAGCGATAACCCAGCGCACt                                       >  1:467550/1‑139 (MQ=255)                                                                                    aCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATg                                    >  1:127052/1‑139 (MQ=255)                                                                                          tGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                              <  2:467550/139‑1 (MQ=255)                                                                                            gCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                            <  1:546703/139‑1 (MQ=255)                                                                                             cAGGCCGTATACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                              >  1:127077/1‑136 (MQ=255)                                                                                             cAGGCCGTATACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                              <  2:127077/136‑1 (MQ=255)                                                                                                gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTc                        >  1:259639/1‑139 (MQ=255)                                                                                                       aCAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATaa                 <  2:742815/139‑1 (MQ=255)                                                                                                               gCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGaa         <  2:127052/139‑1 (MQ=255)                                                                                                               gCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGaa         <  1:497351/139‑1 (MQ=38)                                                                                                                 aTTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                                 <  2:588617/65‑1 (MQ=255)                                                                                                                 aTTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                                 >  1:588617/1‑65 (MQ=255)                                                                                                                  ttCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGc      <  2:792634/139‑1 (MQ=25)                                                                                                                      tGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCctct  >  2:409016/1‑139 (MQ=17)                                                                                                                                     |                                                                                                                            TGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCT  >  NC_000913/2068820‑2069075

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