Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,650,355	T→C	F274L (TTC→CTC)	b1579 →	putative transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,650,355	0	T	C	100.0%	223.7 / NA	84	F274L (TTC→CTC)	b1579	putative transposase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (42/42);  total (42/42)

CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTA  >  NC_000913/1650283‑1650426                                                                         |                                                                        ccAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCa                                                                      <  1:903666/76‑1 (MQ=255) ccAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCa                                                                      >  2:1796937/1‑76 (MQ=255)  cAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCag                                                                     >  1:2991057/1‑76 (MQ=255)  cAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCag                                                                     <  1:4072129/76‑1 (MQ=255)  cAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCag                                                                     <  2:2549688/76‑1 (MQ=255)   aGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCaga                                                                    >  1:2485489/1‑76 (MQ=255)    gCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagag                                                                   >  1:1509539/1‑76 (MQ=255)    gCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagag                                                                   <  2:3227222/76‑1 (MQ=255)     cAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagaga                                                                  <  2:2151632/76‑1 (MQ=255)      aaTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagagag                                                                 >  1:2574420/1‑76 (MQ=255)       aTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGt                                                                >  2:672121/1‑76 (MQ=255)       aTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGt                                                                <  1:1846595/76‑1 (MQ=255)       aTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGt                                                                >  1:1907722/1‑76 (MQ=255)        tCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTa                                                               >  1:1238168/1‑76 (MQ=255)         tgACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTAt                                                              >  2:1920701/2‑76 (MQ=255)          gACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATg                                                             <  1:1406524/76‑1 (MQ=255)            cGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGc                                                           <  2:1900928/76‑1 (MQ=255)             gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGca                                                          >  1:2225678/1‑76 (MQ=255)             gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGca                                                          <  2:2360762/76‑1 (MQ=255)              ctctCCGGAGTCAGTGGACATTAACGAGACTGAGCAAATAGCATATCATTGATGTTCACCTCAGAGAGGATGGcag                                                         <  2:629225/76‑1 (MQ=255)               tctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcaga                                                        >  2:2485251/1‑76 (MQ=255)               tctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcaga                                                        >  1:1915222/1‑76 (MQ=255)               gctcCGGAGTCAGATGACATTAAAGAGACTGAGCAAAGAGCATAGCATTGATGCTCACCACAGACAGTATGGcaca                                                        >  1:864216/2‑74 (MQ=255)                ctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaa                                                       >  1:1014381/1‑76 (MQ=255)                ctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaa                                                       >  1:768898/1‑76 (MQ=255)                cgccGGGGTCAGATGACATTAACGAGTCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaa                                                       <  2:1445764/74‑1 (MQ=255)                 tcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaac                                                      >  1:3817374/1‑76 (MQ=255)                  gcGGAGGCAGATGACATTAACGAGACGGAGCAAAGAGCATATCATTGATGGTCACCTCAGAGAGTATGGcagaaga                                                     <  2:2945/75‑3 (MQ=255)                   cGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacag                                                    >  1:1634883/1‑76 (MQ=255)                   cGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacag                                                    <  2:4497670/76‑1 (MQ=255)                    ggAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacaga                                                   >  1:882876/1‑76 (MQ=255)                     gAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacagaa                                                  <  2:3251224/76‑1 (MQ=255)                       gTCAGCTGATATTAACGAGACTGAGCAATGAGCATATCATTGATTTTCCCCTCAGATGGTGTGGCCGCACAGaaac                                                <  1:3971821/76‑2 (MQ=255)                       gTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaa                                                <  1:3036208/76‑1 (MQ=255)                        tCAGATGACATTAACGAGACTGAGCAACGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCATAACAGaaaaa                                               <  2:2975971/76‑1 (MQ=255)                        tCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaaa                                               >  1:274233/1‑76 (MQ=255)                          aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACa                                             <  2:989175/76‑1 (MQ=255)                          aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACa                                             >  2:1476874/1‑76 (MQ=255)                           gATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaa                                            >  2:2309466/1‑76 (MQ=255)                            aTGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAACACaaa                                           >  2:1122943/1‑76 (MQ=255)                              gACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaaaa                                         >  2:4855138/1‑76 (MQ=255)                               aCATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAt                                        <  1:251420/76‑1 (MQ=255)                                cATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATg                                       <  1:4519291/76‑1 (MQ=255)                                  ttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACGAAAATGCa                                     >  1:3268155/1‑76 (MQ=255)                                  ttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCa                                     <  2:330763/76‑1 (MQ=255)                                  ttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCa                                     >  2:3744442/1‑76 (MQ=255)                                  ttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCa                                     >  1:3132235/1‑76 (MQ=255)                                   tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGACCAGAAAAACAAAAATGCAc                                    <  2:531131/76‑1 (MQ=255)                                   tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCAc                                    <  1:1945323/76‑1 (MQ=255)                                   tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCAc                                    >  1:1879934/1‑76 (MQ=255)                                   tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCAc                                    >  1:2209469/1‑76 (MQ=255)                                   tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCAc                                    <  1:2966094/76‑1 (MQ=255)                                      cGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCtt                                 <  1:459423/76‑1 (MQ=255)                                        tgaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTg                               >  2:3869132/2‑76 (MQ=255)                                          aCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtt                             >  1:2037442/1‑76 (MQ=255)                                          aCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtt                             >  1:4552303/1‑76 (MQ=255)                                           cTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttt                            >  1:2891825/1‑76 (MQ=255)                                              aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTc                         >  2:965489/1‑76 (MQ=255)                                              aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTc                         >  2:842698/1‑76 (MQ=255)                                              aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTc                         <  2:463837/76‑1 (MQ=255)                                               gCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCa                        >  1:2858087/1‑76 (MQ=255)                                               gCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCa                        >  2:4304481/1‑76 (MQ=255)                                                 aaaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAc                      >  2:4363522/1‑76 (MQ=255)                                                   agagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                    >  1:208679/1‑76 (MQ=255)                                                      gCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGaa                 >  1:250481/1‑76 (MQ=255)                                                      gCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGaa                 <  2:1238167/76‑1 (MQ=255)                                                          atCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgtgt             >  2:2252330/1‑76 (MQ=255)                                                           tCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTc            <  2:4563438/76‑1 (MQ=255)                                                           tCATTGATGTTCACCTCAGAGAGTAAGGCAGAACAGAAAAACAAAAAAGCAACTTTGTTTTTTCAACTGAAGTGTc            >  1:2305977/1‑76 (MQ=255)                                                            cATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCa           <  1:4006087/76‑1 (MQ=255)                                                            cATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCa           <  1:2355824/76‑1 (MQ=255)                                                              ttGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGc         <  2:3394058/76‑1 (MQ=255)                                                              ttGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGc         <  1:3639031/76‑1 (MQ=255)                                                              ttGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGc         <  1:3007373/76‑1 (MQ=255)                                                              atGATGTTCACCGCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGc         <  2:3386019/75‑1 (MQ=255)                                                                 aTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgag      <  1:2752780/76‑1 (MQ=255)                                                                 aTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgag      <  1:3776375/76‑1 (MQ=255)                                                                 aTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgag      <  1:328301/76‑1 (MQ=255)                                                                  tGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgaga     <  2:2991057/76‑1 (MQ=255)                                                                   gTTCCCCTCAGAGAGTATGGCAGAACAGAAAAGCAAAAATGCACCTTTGTTTTTCAAGCTGAAGTGTCAGCGagag    <  2:242411/76‑1 (MQ=255)                                                                   gTTCACCTCAGAGAGTATGGCAGAACAGAGACGTGGATGTGGACCTTTGTTTTCCAGGCCGAATTGTCAGCGagag    <  2:239983/76‑1 (MQ=255)                                                                   gTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGagag    <  2:3532308/76‑1 (MQ=255)                                                                    ttCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGt   >  2:3127256/1‑76 (MQ=255)                                                                     tCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTa  <  1:4800319/76‑1 (MQ=255)                                                                         |                                                                        CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTA  >  NC_000913/1650283‑1650426

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