Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,720,772	T→C	R31R (CGT→CGC)	insJ →	IS150 transposase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,720,772	0	T	C	100.0%	263.5 / NA	72	R31R (CGT→CGC)	insJ	IS150 transposase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (39/33);  total (39/33)

CACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGA  >  NC_000913/3720640‑3720905                                                                                                                                     |                                                                                                                                      cACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGgt                                                                                                                                 >  2:1002319/1‑139 (MQ=255)   cccATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGgtgt                                                                                                                               <  2:1289971/139‑1 (MQ=255)                             aTAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATggg                                                                                                     <  1:1002319/139‑1 (MQ=255)                                 ccccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGc                                                                                                 >  1:313732/1‑139 (MQ=255)                                    cTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCct                                                                                              <  2:313732/139‑1 (MQ=255)                                           gTCAAGGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAAcc                                                                                                                     <  1:864497/109‑1 (MQ=255)                                           gTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAAcc                                                                                                                     >  2:864497/1‑109 (MQ=255)                                            tCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaaa                                                                                      >  1:1031982/1‑139 (MQ=255)                                                aGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                                >  1:811836/1‑109 (MQ=255)                                                aGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                                <  2:811836/109‑1 (MQ=255)                                                aGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACa                                                                                  >  1:1095897/1‑139 (MQ=255)                                                  ccAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGAc                                                                                                          >  1:242422/1‑113 (MQ=255)                                                  ccAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGAc                                                                                                          <  2:242422/113‑1 (MQ=255)                                                                     aGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTc                                                                                                               <  2:28884/89‑1 (MQ=255)                                                                     aGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTc                                                                                                               >  1:28884/1‑89 (MQ=255)                                                                      gcgcCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTccc                                                            <  2:448116/139‑1 (MQ=255)                                                                       cgcCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                                <  2:598927/86‑1 (MQ=255)                                                                       cgcCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                                >  1:598927/1‑86 (MQ=255)                                                                         ccTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGAc                                                                                                          <  2:1252257/90‑1 (MQ=255)                                                                         ccTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGAc                                                                                                          >  1:1252257/1‑90 (MQ=255)                                                                         ccTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCaaa                                                         >  1:1349663/1‑139 (MQ=255)                                                                         ccTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCaaa                                                         >  1:1249343/1‑139 (MQ=255)                                                                          cTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAc                                                        <  1:816348/139‑1 (MQ=255)                                                                            tGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCt                                                      >  2:1403298/1‑139 (MQ=255)                                                                                        tCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATggg                                                                                                     <  2:664502/80‑1 (MQ=255)                                                                                        tCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATggg                                                                                                     >  1:664502/1‑80 (MQ=255)                                                                                          actactTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGc                                        <  2:1349663/139‑1 (MQ=255)                                                                                             actTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCct                                                                                              >  2:828609/1‑82 (MQ=255)                                                                                             actTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCct                                                                                              <  1:828609/82‑1 (MQ=255)                                                                                                   caACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCaga                               >  2:484035/1‑139 (MQ=255)                                                                                                    aaCTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCagag                              >  1:727820/1‑139 (MQ=255)                                                                                                       tgatgaTGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGTGTTAGTGCTGATCCAGAGTTg                           <  2:1031982/139‑1 (MQ=255)                                                                                                       tgatgaTGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTg                           <  2:1095897/139‑1 (MQ=255)                                                                                                        gatgatGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaa                                                                                       >  2:555711/1‑78 (MQ=255)                                                                                                        gatgatGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaa                                                                                       <  1:555711/78‑1 (MQ=255)                                                                                                         atgatgGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGcc                                                                                                <  1:900190/68‑1 (MQ=255)                                                                                                         atgatgGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGcc                                                                                                >  2:900190/1‑68 (MQ=255)                                                                                                         atgatgGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                         >  2:721186/1‑139 (MQ=255)                                                                                                         atgatgGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                         >  1:414398/1‑139 (MQ=255)                                                                                                           gatgGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                                <  2:318266/50‑1 (MQ=255)                                                                                                           gatgGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                                >  1:318266/1‑50 (MQ=255)                                                                                                              ggTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTa                    >  1:1405236/1‑139 (MQ=255)                                                                                                                ttACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCct                                                                                              >  2:314209/1‑63 (MQ=255)                                                                                                                ttACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCct                                                                                              <  1:314209/63‑1 (MQ=255)                                                                                                                      gATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTg            <  2:1295202/139‑1 (MQ=255)                                                                                                                        tcatcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                                <  2:20533/37‑1 (MQ=255)                                                                                                                        tcatcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                                >  1:20533/1‑37 (MQ=255)                                                                                                                        tcatcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGt                                                                                                    >  1:893001/1‑49 (MQ=255)                                                                                                                        tcatcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGt                                                                                                    <  2:893001/49‑1 (MQ=255)                                                                                                                        tcatcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGaa          <  2:618626/139‑1 (MQ=255)                                                                                                                        tcatcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGaa          >  1:826466/1‑139 (MQ=255)                                                                                                                        tcatcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGaa          >  1:618626/1‑139 (MQ=255)                                                                                                                          atcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTccc                                                            <  2:1150486/87‑1 (MQ=255)                                                                                                                          atcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTccc                                                            >  1:1150486/1‑87 (MQ=255)                                                                                                                          atcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTaa                                                    <  2:1234142/95‑1 (MQ=255)                                                                                                                          atcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTaa                                                    >  1:1234142/1‑95 (MQ=255)                                                                                                                           tctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTccc                                                            >  2:965839/1‑86 (MQ=255)                                                                                                                           tctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTccc                                                            >  1:1000543/1‑86 (MQ=255)                                                                                                                           tctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTccc                                                            <  2:1000543/86‑1 (MQ=255)                                                                                                                           tctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTaaa                                                   <  2:285899/95‑1 (MQ=255)                                                                                                                           tctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTaaa                                                   >  1:285899/1‑95 (MQ=255)                                                                                                                           tctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGc       >  2:215475/1‑139 (MQ=255)                                                                                                                           tctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGc       >  2:172626/1‑139 (MQ=255)                                                                                                                           tctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGc       >  1:587766/1‑139 (MQ=255)                                                                                                                           tctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGc       >  1:1026995/1‑139 (MQ=255)                                                                                                                           tctcGGCACGCTTTGGTGTCCCCCGAACCCAGGCCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTccc                                                            <  1:965839/86‑1 (MQ=255)                                                                                                                                gCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTg                                         >  2:1100151/1‑100 (MQ=255)                                                                                                                                gCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTg                                         <  1:1100151/100‑1 (MQ=255)                                                                                                                                gCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTg            >  2:113114/1‑129 (MQ=255)                                                                                                                                gCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTg            <  1:113114/129‑1 (MQ=255)                                                                                                                                gCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGa  <  1:172626/139‑1 (MQ=255)                                                                                                                                gCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGa  <  1:1403298/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      CACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGA  >  NC_000913/3720640‑3720905

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