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%	302.4 / NA	66	coding (413/495 nt)	b2655	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (36/30);  total (36/30)

AGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCGTCCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCT  >  NC_000913/2785109‑2785247                                                                        |                                                                    agaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGAc                                                                <  1:732624/76‑1 (MQ=255)   aaaaTATGTCTGAATTTGCACGTTGCCATTATGGGTATATACGAATCCCTCCAGTGACAACTTTCCGCG‑CCGACgg                                                              <  2:1327932/76‑1 (MQ=255)   aaaaTATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACgg                                                              <  1:484406/76‑1 (MQ=255)   aaaaTATGTCTGAATTTGCACGTTGCCATTATGGGGATATACAAATCCCTCCAGTGACAACTTTACGCG‑CCGACgg                                                              >  1:4007534/1‑76 (MQ=255)     aaTATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGAc                                                            >  1:4792163/1‑76 (MQ=255)       tatGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCa                                                          >  2:2478285/1‑76 (MQ=255)       tatGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCa                                                          >  2:3746461/1‑76 (MQ=255)        atGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAg                                                         >  1:4305053/1‑76 (MQ=255)         tGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGa                                                        <  2:3570656/76‑1 (MQ=255)         tGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGa                                                        >  1:2987592/1‑76 (MQ=255)           tCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGaaa                                                      <  1:1869214/76‑1 (MQ=255)            cTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAAc                                                     >  2:1964204/1‑76 (MQ=255)             tGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACt                                                    >  2:1861414/1‑76 (MQ=255)              gtaTTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTc                                                   >  2:1029364/3‑76 (MQ=255)              aaaTTTGCACGTTGCCATTATGTTTATATACAATTCCCTCCCGTGAGCACTTTCCGCG‑CCGGCGGACCAGAAACTc                                                   <  1:3690182/75‑1 (MQ=255)                aTTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTccc                                                 >  1:3841075/1‑76 (MQ=255)                 tttGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCg                                                >  1:947835/1‑76 (MQ=255)                 tttGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCg                                                >  1:3464024/1‑76 (MQ=255)                 tttGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCg                                                >  2:1720708/1‑76 (MQ=255)                   tGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaa                                              <  2:646919/76‑1 (MQ=255)                    gCACTTTGCCATTATGGGTATATACAACTCCCTCCAGTGACACCTTTCCGCG‑CCGACGGACCAGAAACTCCCgaag                                             <  1:2846479/76‑1 (MQ=255)                    gCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaag                                             >  1:3839616/1‑76 (MQ=255)                    gCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaag                                             >  1:1605521/1‑76 (MQ=255)                     cACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaaga                                            >  2:2680722/1‑76 (MQ=255)                     cACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaaga                                            >  2:1295162/1‑76 (MQ=255)                      aCGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaa                                           >  1:299888/1‑76 (MQ=255)                      aCGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaa                                           <  1:4750349/76‑1 (MQ=255)                      aCGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaa                                           <  1:1586978/76‑1 (MQ=255)                         ttGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaagaa                                        <  1:3723459/76‑1 (MQ=255)                         ttGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaagaa                                        <  2:3556528/76‑1 (MQ=255)                          tGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGaaa                                       <  2:3222114/76‑1 (MQ=255)                          tGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGaaa                                       >  1:1580966/1‑76 (MQ=255)                            ccATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAg                                     >  1:1879103/1‑76 (MQ=255)                            ccATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAg                                     >  2:438809/1‑76 (MQ=255)                             cATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAgg                                    <  2:892613/76‑1 (MQ=255)                              aTTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAggg                                   >  2:1084117/1‑76 (MQ=255)                              aTTATGGGTATATAAAAATTACTCCAGTGACAACTTTCCGCG‑CCGACGGGCCAGACACTCCCGAAGAAGACAAggg                                   >  1:2249908/1‑76 (MQ=255)                                tATGGGTATATACAAATCCCTCCAGTGACAGCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGtt                                 <  2:2102888/76‑1 (MQ=255)                                tATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGtt                                 <  2:2587809/76‑1 (MQ=255)                                tATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGtt                                 >  2:2141183/1‑76 (MQ=255)                                  tGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTAc                               >  1:2448141/1‑76 (MQ=255)                                  tGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTAc                               >  1:3226657/1‑76 (MQ=255)                                   gggTATATACAAGTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAGGAAAAGGGTTACt                              <  1:215608/76‑1 (MQ=255)                                   gggTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACt                              <  1:697867/76‑1 (MQ=255)                                    ggTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTg                             <  2:3839616/76‑1 (MQ=255)                                     gTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTgg                            >  2:469429/1‑76 (MQ=255)                                      tatataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGt                           >  1:1041609/1‑76 (MQ=255)                                         ataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGCTACTGGttta                        >  2:3466571/1‑75 (MQ=255)                                              aaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGAGCAGAAGCTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCt                   <  2:3602983/76‑1 (MQ=255)                                              aaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGCAAAGGGTTACTGGTTTCATGCt                   >  1:986068/1‑76 (MQ=255)                                                  ccTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATc               <  2:2916663/76‑1 (MQ=255)                                                  ccTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATc               <  2:1257800/76‑1 (MQ=255)                                                  ccTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATc               >  2:3257136/1‑76 (MQ=255)                                                  ccTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTGATc               <  1:878611/76‑1 (MQ=255)                                                    tCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCaa             <  2:2064603/76‑1 (MQ=255)                                                    tCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCaa             >  1:4441599/1‑76 (MQ=255)                                                    tCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGCTTCATGCTTATCaa             >  1:1292480/1‑76 (MQ=255)                                                      cAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAcc           <  1:3061298/76‑1 (MQ=255)                                                       aGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAccc          >  1:3370656/1‑76 (MQ=255)                                                       aGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAccc          <  1:246147/76‑1 (MQ=255)                                                       aGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAccc          >  2:795247/1‑76 (MQ=255)                                                         tGACAACTTTCCGCG‑CCGACGGACCAGAAACTCTCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGa        <  1:2159568/76‑1 (MQ=255)                                                           aCAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAg      >  1:3107190/1‑76 (MQ=255)                                                            cAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGa     <  1:371190/76‑1 (MQ=255)                                                            cAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGa     <  1:111182/76‑1 (MQ=255)                                                               cTTTCCGCG‑CCGACGGACCAGACACTCCCGAAGAATAAAAGGGTTACTGTTTTCATGCTTATCAGCCCGAATATCt  <  2:2550684/76‑1 (MQ=255)                                                                        |                                                                    AGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCGTCCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCT  >  NC_000913/2785109‑2785247

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