Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,531,361	T→C	H188R (CAT→CGT)	b1462 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,531,361	0	T	C	100.0%	220.6 / NA	70	H188R (CAT→CGT)	b1462	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (44/26);  total (44/26)

TTTGATAAAGAAACCTGGTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGATGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAAC  >  NC_000913/1531273‑1531450                                                                                         |                                                                                          tttGATAAAGAAACCTGGTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCaa                                                                                       >  1:609015/1‑93 (MQ=255) tttGATAAAGAAACCTGGTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCaa                                                                                       >  1:554248/1‑93 (MQ=255)   tGATAAAGAAACCTGGTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACg                                                                                     >  2:2091399/1‑93 (MQ=255)    gATAAAGAAACCTGGTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGt                                                                                    >  2:3523220/1‑93 (MQ=255)      tAAAGAAACCTGGTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTAt                                                                                  >  2:1681132/1‑93 (MQ=255)      tAAAGAAACCTGGTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTAt                                                                                  >  1:4334810/1‑93 (MQ=255)            aaCCTGGTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCt                                                                            >  1:1969309/1‑93 (MQ=255)               cTGGTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTat                                                                         >  1:2498583/1‑93 (MQ=255)               cTGGTTATAACATGAATTCTACCCCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTat                                                                         <  2:868551/93‑1 (MQ=255)                tGGTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTatc                                                                        >  2:1859198/1‑93 (MQ=255)                  gTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTatcat                                                                      >  1:1452436/1‑93 (MQ=255)                      taGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTatcatcatc                                                                  >  2:2823331/1‑93 (MQ=255)                      taGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTatcatcatc                                                                  >  1:2319076/1‑93 (MQ=255)                      taGCATGAATGCTTAACCCGCCGAAACACGCTTGCCGCTGGTAACAAACGTCCCAGTACCTAAATGACGCAACGTATTTAGCTTatcatcatc                                                                  >  2:3505858/1‑93 (MQ=255)                        gCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCaa                                                                >  2:4281594/1‑93 (MQ=255)                           tGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACt                                                             >  2:1889324/1‑93 (MQ=255)                             aaTGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGTCGCAACGTATTTAGCTTATCATCATCAAACTGc                                                           >  1:427654/1‑93 (MQ=255)                               tGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCa                                                         >  1:1781766/1‑93 (MQ=255)                                  ttAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTACCAGCg                                                      >  1:2721692/1‑93 (MQ=255)                                    aaCCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGc                                                    >  2:1194463/1‑93 (MQ=255)                                      cccGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCt                                                  <  2:347582/93‑1 (MQ=255)                                      cccGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCt                                                  <  2:4147531/93‑1 (MQ=255)                                      cccGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCt                                                  <  1:1712623/93‑1 (MQ=255)                                      cccGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCt                                                  <  1:1302647/93‑1 (MQ=255)                                      cccGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCt                                                  <  2:1343361/93‑1 (MQ=255)                                       ccgccgTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCtt                                                 >  2:3357922/1‑93 (MQ=255)                                          ccgTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCATCGGCCTTCGa                                              >  1:2825486/1‑93 (MQ=255)                                           cgTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGAc                                             <  1:795792/93‑1 (MQ=255)                                            gTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACa                                            >  2:882321/1‑93 (MQ=255)                                            gTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACa                                            >  1:4474711/1‑93 (MQ=255)                                             tAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACaa                                           >  1:2764320/1‑93 (MQ=255)                                             tAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACaa                                           >  1:986187/1‑93 (MQ=255)                                              aaCACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACaaa                                          <  2:2978605/93‑1 (MQ=255)                                              aaCACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACaaa                                          <  1:791838/93‑1 (MQ=255)                                               acacGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAAt                                         >  2:3649473/1‑93 (MQ=255)                                                cacGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATa                                        <  1:134076/93‑1 (MQ=255)                                                 acGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATAc                                       <  1:3925355/93‑1 (MQ=255)                                                     ttGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGt                                   <  2:3342834/93‑1 (MQ=255)                                                     ttGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGt                                   >  2:1848662/1‑93 (MQ=255)                                                     ttGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGt                                   <  2:802240/93‑1 (MQ=255)                                                     ttGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGt                                   >  1:1234584/1‑93 (MQ=255)                                                     ttGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGt                                   >  1:401473/1‑93 (MQ=255)                                                     ttGCCGCTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGt                                   >  2:1479211/1‑93 (MQ=255)                                                           cTGGTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCt                             <  1:2001845/93‑1 (MQ=255)                                                             ggTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATGCCCGTGCGTCTGc                           <  2:1096051/93‑1 (MQ=255)                                                              gTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGcc                          >  2:148077/1‑93 (MQ=255)                                                              gTAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGcc                          <  2:4635179/93‑1 (MQ=255)                                                               tAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCg                         >  1:977218/1‑93 (MQ=255)                                                               tAACAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCg                         >  2:3902228/1‑93 (MQ=255)                                                                 aCAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCgct                       >  1:1647798/1‑93 (MQ=255)                                                                  cAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCgctg                      <  2:2708544/93‑1 (MQ=255)                                                                  cAAACGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCgctg                      >  1:3729158/1‑93 (MQ=255)                                                                   aaaCGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCgctgc                     <  1:4468910/93‑1 (MQ=255)                                                                   aaaCGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCgctgc                     <  1:3171578/93‑1 (MQ=255)                                                                   aaaCGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCgctgc                     >  2:151991/1‑93 (MQ=255)                                                                   aaaCGTCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCgctgc                     <  1:4311106/93‑1 (MQ=255)                                                                       gTCCCAACACCTAAATGACCCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTCCGACAAATCCCCGTGCGTCTGCCGCTGCTGAt                 <  1:471350/93‑1 (MQ=255)                                                                        tCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGAta                <  1:1709537/93‑1 (MQ=255)                                                                        tCCCAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGAta                <  1:1769937/93‑1 (MQ=255)                                                                          ccAGCACCTAAATGACGCAACGTATTGAGCTTATCCTCATCAAACTGCCAGCGGCCTTCGACAAATACCCCTGCGTCTGCCGCTGCTGAtact              >  1:291316/1‑93 (MQ=255)                                                                          ccAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGAtact              >  2:3132566/1‑93 (MQ=255)                                                                          ccAGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGAtact              <  2:2980185/93‑1 (MQ=255)                                                                            aGCACCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGAtactac            <  2:71950/93‑1 (MQ=255)                                                                               aCCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTc         >  1:4411902/1‑93 (MQ=255)                                                                               aCCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTc         >  1:80301/1‑93 (MQ=255)                                                                               aCCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTc         >  2:584834/1‑93 (MQ=255)                                                                               aCCTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTc         <  2:1772520/93‑1 (MQ=255)                                                                                 cTAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGc       <  1:96413/93‑1 (MQ=255)                                                                                  tAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGcc      >  1:4339588/1‑93 (MQ=255)                                                                                  tAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGcc      >  2:2909087/1‑93 (MQ=255)                                                                                  tAAATGACGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGcc      >  2:2455746/1‑93 (MQ=255)                                                                                      tgacgCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAAc  <  2:2720344/93‑1 (MQ=255)                                                                                      tgacgCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAAc  >  1:1069368/1‑93 (MQ=255)                                                                                         |                                                                                          TTTGATAAAGAAACCTGGTTATAGCATGAATGCTTAACCCGCCGTAACACGCTTGCCGCTGGTAACAAACGTCCCAGCACCTAAATGATGCAACGTATTGAGCTTATCATCATCAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAAC  >  NC_000913/1531273‑1531450

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