Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,928,343	G→A	A43T (GCA→ACA)	rho →	transcription termination factor Rho

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,928,343	0	G	A	100.0%	187.2 / NA	50	A43T (GCA→ACA)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/23);  total (27/23)

CGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCC  >  REL606/3928270‑3928412                                                                          |                                                                      cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGcacaca                                                                     >  1:979048/1‑76 (MQ=255)  gAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACaa                                                                    >  2:900213/1‑76 (MQ=255)  gAAAATATGGGGCTGGAAAACCTGGCTCGAAAGGGTAAGCAGGAGATTATTTTTTGCATCCTGCACCAGCACACaa                                                                    >  1:837364/1‑76 (MQ=255)    aaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAag                                                                  >  1:419335/1‑76 (MQ=255)      atatGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAagag                                                                >  2:598177/1‑76 (MQ=255)       tatGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGt                                                               <  1:152887/76‑1 (MQ=255)       tatGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGt                                                               >  2:1072041/1‑76 (MQ=255)         tGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTgg                                                             >  1:666051/1‑76 (MQ=255)           gggCTGGAAAACCTGGCTCGTATGCGTAAGCGGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCg                                                           <  2:241646/76‑1 (MQ=255)             gCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGaa                                                         <  1:478275/76‑1 (MQ=255)              cTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAg                                                        <  1:444464/76‑1 (MQ=255)               tGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGa                                                       <  2:808230/76‑1 (MQ=255)               tGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGa                                                       <  2:1033828/76‑1 (MQ=255)                ggAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGat                                                      <  1:223364/76‑1 (MQ=255)                ggAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGat                                                      <  1:223369/76‑1 (MQ=255)                     aCCTGGCTCGTATGCGTAAGCANGACATTATTTGGGCCNTCCTGAAGCAGCACACAACGGGTGGCGAAAANGTCtt                                                 >  1:26411/1‑76 (MQ=255)                      ccTGGCTCGTATGCGTAAGCAGGACATTATTTTTGTCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCttt                                                >  2:466286/1‑76 (MQ=255)                       cTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTg                                               <  1:994116/76‑1 (MQ=255)                       cTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTg                                               <  1:1033058/76‑1 (MQ=255)                          gCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTg                                            >  1:832170/1‑76 (MQ=255)                          gCTCGTATGCGTAAGCAGGACATTATTTGTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTg                                            <  2:847702/76‑1 (MQ=255)                             cGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATg                                         <  2:610861/76‑1 (MQ=255)                             cGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATg                                         >  1:582151/1‑76 (MQ=255)                              gTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATgg                                        >  2:662450/1‑76 (MQ=255)                                aTGCGTAAGCAGGACATTATTTTTGCCATCCTNAAACAGCANACAAACAGTGGCGAAGATANCTTTGGTGATGGCg                                      >  2:1189446/1‑76 (MQ=255)                                aTGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGNTGATGGCg                                      >  1:1267832/1‑76 (MQ=255)                                aTGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCg                                      <  1:776700/76‑1 (MQ=255)                                aTGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCg                                      >  1:775836/1‑76 (MQ=255)                                    gTAAGCAGGACATTATTTTTGCCNTCCTGNAGCAGCACACAAAGNGTAGCGTAGATATCTTAGGTGATGGCGTACt                                  >  2:192110/1‑76 (MQ=255)                                    gTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACt                                  <  1:989387/76‑1 (MQ=255)                                    gTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACt                                  <  1:989383/76‑1 (MQ=255)                                     tAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTg                                 >  1:567416/1‑76 (MQ=255)                                     tAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTg                                 >  2:1434880/1‑76 (MQ=255)                                      aaGCAGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTgg                                <  2:1245321/76‑1 (MQ=255)                                          aGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGat                            >  2:816506/1‑76 (MQ=255)                                          aGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGat                            >  2:232242/1‑76 (MQ=255)                                          aGGACATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGNACTGGAGat                            <  1:35691/76‑1 (MQ=255)                                             aCATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATAtt                         <  1:1145177/76‑1 (MQ=255)                                             aCATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATAtt                         >  2:638521/1‑76 (MQ=255)                                              cATTATTTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTg                        >  2:110751/1‑76 (MQ=255)                                                ttattTTTGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCa                      >  1:465372/1‑76 (MQ=255)                                                     tttGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatg                 <  2:893513/76‑1 (MQ=255)                                                       tGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgga               >  2:1098030/1‑76 (MQ=255)                                                       tGCCATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgga               >  2:1041274/1‑76 (MQ=255)                                                       tGCCATCCNGTAGCAGCACNCAAAGAGTGGCGAAGATATCTTNGGTGATGGCGTACTGGAGATATTGCAggatgga               <  2:1211318/76‑1 (MQ=255)                                                        gCCATCCTGAAGCAGCACACAAAGAGTCGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatggat              >  2:99029/1‑76 (MQ=255)                                                          cATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGAttt            <  2:271294/76‑1 (MQ=255)                                                          cATCCTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGAttt            >  1:748215/1‑76 (MQ=255)                                                               tGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGttt       >  2:1260661/1‑76 (MQ=255)                                                               tGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGttt       <  2:1110310/76‑1 (MQ=255)                                                                gAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtc      >  2:378292/1‑76 (MQ=255)                                                                gAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtc      <  1:1252058/76‑1 (MQ=255)                                                                gAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtc      <  1:1141181/76‑1 (MQ=255)                                                                    cagcaCACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcctcc  >  1:306615/1‑76 (MQ=255)                                                                          |                                                                      CGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCC  >  REL606/3928270‑3928412

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