Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,785,180	Δ1 bp	coding (413/495 nt)	b2655 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,785,180	0	T	.	100.0%	420.8 / NA	92	coding (413/495 nt)	b2655	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (44/48);  total (44/48)

AGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCGTCCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTG  >  NC_000913/2785109‑2785250                                                                        |                                                                       agaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGAc                                                                   <  2:12434/76‑1 (MQ=255) agaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGAc                                                                   <  1:1549336/76‑1 (MQ=255) agaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGAc                                                                   <  2:5840695/76‑1 (MQ=255) agaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGAc                                                                   <  2:1130723/76‑1 (MQ=255)  gaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACg                                                                  >  2:3124766/1‑76 (MQ=255)  gaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACg                                                                  <  2:4259254/76‑1 (MQ=255)  gaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACg                                                                  <  1:2715613/76‑1 (MQ=255)  gaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACg                                                                  >  2:1610338/1‑76 (MQ=255)      atatGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGAcc                                                              <  1:2385093/76‑1 (MQ=255)       tatGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCa                                                             <  2:740771/76‑1 (MQ=255)       tatGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCa                                                             >  1:669299/1‑76 (MQ=255)        atGTCTGAATTTGCGCGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAg                                                            <  2:2239895/76‑1 (MQ=255)        atGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAg                                                            >  1:5695122/1‑76 (MQ=255)        atGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAg                                                            >  1:1390821/1‑76 (MQ=255)        atGTCTGAATGTGCCAGTTGCCATTAGGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAg                                                            <  2:5691232/76‑1 (MQ=255)         tGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGa                                                           >  2:4272845/1‑76 (MQ=255)           tCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGaaa                                                         >  2:658637/1‑76 (MQ=255)            cTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAAc                                                        <  2:2532267/76‑1 (MQ=255)             tGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACt                                                       >  1:4966824/1‑76 (MQ=255)             tGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACt                                                       >  2:2703189/1‑76 (MQ=255)              gAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTc                                                      <  1:1487569/76‑1 (MQ=255)              gAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTc                                                      >  1:4874782/1‑76 (MQ=255)              gAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTc                                                      <  1:4068352/76‑1 (MQ=255)              gAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTc                                                      >  1:4138470/1‑76 (MQ=255)                aTTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTccc                                                    <  2:1341111/76‑1 (MQ=255)                 tttGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCg                                                   <  2:5298372/76‑1 (MQ=255)                 nttGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCg                                                   >  1:2242870/2‑76 (MQ=255)                  ttGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCga                                                  <  2:3906754/76‑1 (MQ=255)                  ttGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCga                                                  >  2:5280779/1‑76 (MQ=255)                   tGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaa                                                 >  2:4168100/1‑76 (MQ=255)                   tGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaa                                                 <  2:4493362/76‑1 (MQ=255)                    gcgcGTTGCCATTATGTGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaag                                                <  2:1491186/73‑1 (MQ=255)                    gCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaag                                                <  2:2219593/76‑1 (MQ=255)                    gCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaag                                                <  2:2128216/76‑1 (MQ=255)                     cACGTTGCCATGAGGGGGATATACAAATCCCGCCAGTGACAACTTGCCGCG‑CCGCCGGACGAGAAACTCCcgacga                                               >  2:3365894/1‑76 (MQ=255)                      gcGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGCCGGACCAGAAACTCCCgaataa                                              <  2:4334721/75‑1 (MQ=255)                      aCGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaa                                              >  1:3552863/1‑76 (MQ=255)                      aCGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaa                                              <  1:1527096/76‑1 (MQ=255)                       ctttGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaag                                             <  2:1912703/74‑1 (MQ=255)                       cGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaag                                             <  2:3589375/76‑1 (MQ=255)                       cGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaag                                             >  2:1368903/1‑76 (MQ=255)                        gatGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTGCCGAG‑CCGACGGACCAGAAACGCCCgaagaaga                                            >  1:5403191/3‑76 (MQ=255)                        gTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaaga                                            >  2:1541658/1‑76 (MQ=255)                         ttGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaagaa                                           >  1:4904384/1‑76 (MQ=255)                            ccATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAg                                        <  2:4874781/76‑1 (MQ=255)                             ctgtATGGGTGTATACAAAGCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAgg                                       <  2:543262/73‑1 (MQ=255)                             cATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAgg                                       <  2:5651739/76‑1 (MQ=255)                             cATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAgg                                       <  1:3124766/76‑1 (MQ=255)                              aTTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAggg                                      >  2:1792937/1‑76 (MQ=255)                              aTTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAggg                                      <  2:1804955/76‑1 (MQ=255)                              aTTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAggg                                      <  2:4904383/76‑1 (MQ=255)                               ttATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGt                                     >  1:2350925/1‑76 (MQ=255)                                tATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGtt                                    <  1:1151387/76‑1 (MQ=255)                                  tGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTAc                                  >  2:3859647/1‑76 (MQ=255)                                   gggTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACt                                 <  2:5603870/76‑1 (MQ=255)                                      tatataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGt                              <  1:4964486/76‑1 (MQ=255)                                         ataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTc                           <  1:1849757/76‑1 (MQ=255)                                         ataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTc                           >  1:4607694/1‑76 (MQ=255)                                          taCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGATACTGGTTTCa                          >  1:3407412/1‑76 (MQ=255)                                           aCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCAt                         >  1:5239727/1‑76 (MQ=255)                                            cAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAGGAAAGGGGAGACTGGTTTCATg                        >  2:5325466/1‑76 (MQ=255)                                             aaaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGc                       >  2:2219819/1‑76 (MQ=255)                                               accccTCCAATGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCtt                     <  2:133418/74‑1 (MQ=255)                                                tCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTa                    <  2:3177559/76‑1 (MQ=255)                                                  ccTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATc                  >  2:3798123/1‑76 (MQ=255)                                                  ccTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATc                  <  1:1610338/76‑1 (MQ=255)                                                   cTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCa                 >  1:5670458/1‑76 (MQ=255)                                                     ccAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAc               <  2:3352937/76‑1 (MQ=255)                                                      cAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAcc              <  1:5280780/76‑1 (MQ=255)                                                       agggACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAccc             <  2:5530073/73‑1 (MQ=255)                                                       aGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAccc             <  2:1026886/76‑1 (MQ=255)                                                       aGTGACAACTTTCCGCG‑CCGACGGACCAGAAACGCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAccc             <  2:4811785/76‑1 (MQ=255)                                                        gTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCg            <  1:1541658/76‑1 (MQ=255)                                                        gTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCg            >  1:2168409/1‑76 (MQ=255)                                                        gTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCg            >  2:4556541/1‑76 (MQ=255)                                                        gTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCg            >  1:2258895/1‑76 (MQ=255)                                                         tGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGa           >  1:5714996/1‑76 (MQ=255)                                                          gACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGCAAATGTGTACTGGTTTCATGCTTATCAACCCGaa          >  1:355890/1‑76 (MQ=255)                                                          gACAACTTTCCGCG‑CCGACGGACCAGAAACTACCGAAGAAGAAAAGGGTTACTGGGTTCATGCTTATCAACCCGaa          >  1:313431/1‑76 (MQ=255)                                                           aCAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAg         >  1:1625345/1‑76 (MQ=255)                                                            cAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCa                 >  1:1887920/1‑67 (MQ=255)                                                            cAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCa                 <  2:1887920/67‑1 (MQ=255)                                                            cAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAgt        <  1:606139/76‑2 (MQ=255)                                                             aaCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAccc             <  2:4370589/70‑1 (MQ=255)                                                             aaCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAccc             >  1:4370590/1‑70 (MQ=255)                                                             aaCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGAt       >  1:3315871/1‑76 (MQ=255)                                                               cTTTCCGCG‑CCGACGGACCAGAACCTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCt     <  1:4714018/76‑1 (MQ=255)                                                                tttCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCtt    >  2:891613/1‑76 (MQ=255)                                                                tttCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCtt    >  2:903129/1‑76 (MQ=255)                                                                  tCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTg  <  2:5668395/76‑1 (MQ=255)                                                                  tCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTg  >  2:5136906/1‑76 (MQ=255)                                                                  tCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTg  >  2:2934282/1‑76 (MQ=255)                                                                        |                                                                       AGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCGTCCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTG  >  NC_000913/2785109‑2785250

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