Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,021,896 (T)9→8 coding (789/1140 nt) ECB_02821 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,021,8880T.97.1% 161.9 / ‑1.0 35coding (797/1140 nt)ECB_02821hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (12/22);  total (12/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

ATTCTTCGAGAGTGAAACTGTACCTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACCTTTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACA  >  REL606/3021813‑3021959
                                                                           |                                                                       
aTTCTTCGAGAGTGAAACTGTACCTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACCt                                                                         >  1:311152/1‑76 (MQ=255)
  tcttcGAGAGTGAAACTGTACCTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACCttt                                                                       >  1:1050455/1‑76 (MQ=255)
   cttcGAGAGTGAAACTGTACCTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACCtttt                                                                      >  1:635035/1‑76 (MQ=255)
        agagTGAAACTGTACCTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACCttttttttc                                                                 >  2:22378/1‑75 (MQ=255)
        agagTGAAACTGTACCTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACCttttttttc                                                                 >  2:1220799/1‑75 (MQ=255)
              aaaCTGTACCTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCtt                                                          >  2:1638273/1‑76 (MQ=255)
                   gTACCTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCaa                                                     >  2:1311433/1‑76 (MQ=255)
                    tACCTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCaaa                                                    <  2:1106007/76‑1 (MQ=255)
                     aCCTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCaaaa                                                   <  2:245511/76‑1 (MQ=255)
                       cTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCaaaaaa                                                 <  1:844886/76‑1 (MQ=255)
                         ccATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACa                                               <  2:1412179/76‑1 (MQ=255)
                         ccATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACa                                               <  2:1412176/76‑1 (MQ=255)
                           aTAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACaaa                                             >  1:1003376/1‑76 (MQ=255)
                            tAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACaaaa                                            >  2:1542410/1‑76 (MQ=255)
                              gCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAaata                                          <  1:763588/76‑1 (MQ=255)
                               cAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAaataa                                         <  1:1401150/76‑1 (MQ=255)
                                aGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACCTTTTTTTTTCTGCCTTTCCAAAAAACAAaataa                                         <  2:1486880/76‑1 (MQ=255)
                                aGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAaataat                                        <  2:191597/76‑1 (MQ=255)
                                aGAAAGATCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAaataat                                        <  1:1228171/76‑1 (MQ=255)
                                    atatCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAataataata                                    <  1:517747/76‑1 (MQ=255)
                                    atatCAGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAataataata                                    <  1:1553614/76‑1 (MQ=255)
                                         aGTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTc                               >  2:118216/1‑76 (MQ=255)
                                          gTAATTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCg                              <  1:606512/76‑1 (MQ=255)
                                            aaTTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCATa                            <  1:587896/76‑1 (MQ=255)
                                             aTTGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTAc                           <  1:1823975/76‑1 (MQ=255)
                                               tGCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACaa                         >  1:1654913/1‑76 (MQ=255)
                                                gCAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACaaa                        >  1:806056/1‑76 (MQ=255)
                                                 cAAATGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAat                       >  1:1533625/1‑76 (MQ=255)
                                                  aaaTGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAatc                      >  2:81963/1‑76 (MQ=255)
                                                    aTGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAatcca                    <  2:721308/76‑1 (MQ=255)
                                                     tGAAACGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAACAATATTGTCGTACAAatccat                   <  2:206572/76‑1 (MQ=255)
                                                     tGAAACGCCGATATTTTCTACC‑TTTTGTTTCTGCCTTTCCAAAAAACAAAATAACAATATTGTCGTACAAatccat                   <  2:206584/76‑1 (MQ=255)
                                                       aaaCGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAatccatcc                 <  1:181027/76‑1 (MQ=255)
                                                         aCGCCGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCctc               <  1:1526152/76‑1 (MQ=255)
                                                            ccGATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTAc            <  1:1668661/76‑1 (MQ=255)
                                                              gATATTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCa          >  1:1561156/1‑76 (MQ=255)
                                                                 atTTTCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGa       >  1:424628/1‑76 (MQ=255)
                                                                  ttttCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGaa      <  2:1682409/76‑1 (MQ=255)
                                                                     tCTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAAc   <  1:1460195/76‑1 (MQ=255)
                                                                      cTACC‑TTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACa  >  2:1378609/1‑76 (MQ=255)
                                                                           |                                                                       
ATTCTTCGAGAGTGAAACTGTACCTCCATAGCAGAAATATCAGTAATTGCAAATGAAACGCCGATATTTTCTACCTTTTTTTTTCTGCCTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACA  >  REL606/3021813‑3021959

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