Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,733,506 (A)7→8 coding (4420/4617 nt) lhr → putative ATP‑dependent helicase Lhr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,733,4991.A92.9% 92.7 / ‑0.5 28coding (4413/4617 nt)lhrputative ATP‑dependent helicase Lhr
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (10/16);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCA  >  NC_000913/1733420‑1733571
                                                                                |                                                                        
cGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCa                                                                         <  1:311480/81‑1 (MQ=255)
  cTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCaaa                                                                       <  2:311481/81‑1 (MQ=255)
  cTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCaa                                                                        >  1:311515/1‑80 (MQ=255)
  cTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCaa                                                                        >  2:311529/1‑80 (MQ=255)
    ggTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTCTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGCGGCaaaa                                                                      >  2:311528/1‑80 (MQ=255)
        ccAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCaaaaaaaat                                                                 >  2:311513/1‑79 (MQ=255)
            cgcgTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTg                                                              <  2:311482/81‑1 (MQ=255)
             gcgTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTgg                                                             >  1:311522/1‑81 (MQ=255)
                   ggcgggTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTgg                                                        >  2:311520/1‑80 (MQ=255)
                       ggTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGaa                                                  <  2:311483/81‑1 (MQ=255)
                                gtggtCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAgg                                           <  2:311484/80‑1 (MQ=255)
                                     cGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTAc                                     >  2:311521/1‑81 (MQ=255)
                                       tttCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTc                                   <  1:311485/81‑1 (MQ=255)
                                       tttCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTc                                   <  1:311486/81‑1 (MQ=255)
                                         tCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCg                                  <  1:311487/80‑1 (MQ=255)
                                            ggcggcAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGccc                               >  1:311519/1‑80 (MQ=255)
                                               ggcAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCaga                           <  1:311488/81‑1 (MQ=255)
                                                gcAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAATGAGGAATTACTCGCCCCagag                          >  1:311523/1‑81 (MQ=255)
                                                 cAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAgg                         <  1:311489/81‑1 (MQ=255)
                                                  aaaTTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGt                        <  2:311490/81‑1 (MQ=255)
                                                   aaTTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGtt                       <  2:311491/81‑1 (MQ=255)
                                                     ttGTTACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGtttt                     <  1:311492/81‑1 (MQ=255)
                                                        ttACTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCa                  <  1:311493/81‑1 (MQ=255)
                                                          aCTCTATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCAc                 <  2:311494/80‑1 (MQ=255)
                                                             ctATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACgcgc             >  1:311539/1‑81 (MQ=255)
                                                             ctATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACgcgc             >  1:311527/1‑81 (MQ=255)
                                                              tATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCAcgcg              <  1:311495/79‑1 (MQ=255)
                                                              tATCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCt            >  2:311534/1‑81 (MQ=255)
                                                                tCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCt            <  1:311497/79‑1 (MQ=255)
                                                                tCTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCt            <  1:311496/79‑1 (MQ=255)
                                                                 cTGGCGCAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGAc         >  2:311526/1‑81 (MQ=255)
                                                                    gcgtAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAGTTACTCGCCCCAGAGGTTTTCCACGCGCTGACt        <  2:311498/78‑1 (MQ=255)
                                                                     cgcAAGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTAcc     <  1:311499/81‑1 (MQ=255)
                                                                        aaGGTGGCAAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCa  >  1:311532/1‑81 (MQ=255)
                                                                                |                                                                        
CGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGTCTGGCAGGAAAAAGAGGAATTACTCGCCCCAGAGGTTTTCCACGCGCTGACTACCGCA  >  NC_000913/1733420‑1733571

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