Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,733,378 A→G 10.2% N1431S (AAT→AGT)  lhr → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,733,3780AG10.2% 106.6 / 4.3 39N1431S (AAT→AGT) lhrputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (16/19);  new base G (2/2);  total (18/21)
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.49e-01

AATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGT  >  NC_000913/1733232‑1733513
                                                                                                                                                  |                                                                                                                                       
aaTTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAAt                                                                                                                                        <  4:202748/148‑1 (MQ=255)
 aTTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAAtg                                                                                                                                       >  6:131340/1‑148 (MQ=255)
  ttttGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAAtgt                                                                                                                                      <  5:395431/145‑1 (MQ=255)
  ttATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAAtgt                                                                                                                                      >  8:71331/1‑148 (MQ=255)
  ttATGCGAGGTCGTTTTGGAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAAtgt                                                                                                                                      >  4:244025/1‑148 (MQ=255)
              gTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAAtct                                                                                                                          <  1:408125/148‑1 (MQ=255)
              gTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAAtct                                                                                                                          <  3:401798/148‑1 (MQ=255)
                       aaGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAGTGTGTGGGGAAATCTTCTGCCCTg                                                                                                                 >  8:9111/1‑148 (MQ=255)
                       aaGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAGTGTGTGGGGAAATCTTCTGCCCTg                                                                                                                 >  8:2902/1‑148 (MQ=255)
                        aGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTg                                                                                                                      <  1:267850/142‑1 (MQ=255)
                        aGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTg                                                                                                                      >  2:267850/1‑142 (MQ=255)
                                     gcgcAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCgg                                                                                                   <  7:271881/148‑1 (MQ=255)
                                     gcgcAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCg                                                                                                    >  4:256406/1‑147 (MQ=255)
                                              gCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGCCACCAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGCTCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCCACGCTgg                                                                                          >  3:319868/1‑148 (MQ=255)
                                                  aaCGTCTGACTATTGCCCGATTGCGCTATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGtt                                                                                        <  7:210166/146‑1 (MQ=255)
                                                          aCTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAGTGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAAcgcg                                                                                <  7:9111/146‑1 (MQ=255)
                                                          aCTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAGTGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAAcgcg                                                                                <  7:2902/146‑1 (MQ=255)
                                                              ttGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTcgggc                                                                          <  2:426971/148‑1 (MQ=255)
                                                                 aCCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCGCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTcgggcggg                                                                       <  5:74554/148‑1 (MQ=255)
                                                                  ccGATTGCGCGATCTCGCGACACACCCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCCACGATCCGCCTAATGTCTGGGGAAATCTTCTGCCCTGGCCTGCACAACCGGCAACGCTGCTTCCAACGCGTcgggcgc                                                                        >  4:352479/1‑145 (MQ=255)
                                                                     aTTGCGCGATCTGCCGACACAAGCCACGCAAACGCGCCACTATACACCAGGGGCGCTCTCTGCCACCGATCCGGCTAATGTGTGGGGACCTCTTCTGCCCTGGCCGGCCCATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGGgcg                                                                   <  6:80428/148‑1 (MQ=255)
                                                                           cgATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGGACATCCGGCCACGCTGGTTCCAACGCGTCGGGCGGGTGCGCtggtg                                                             <  5:131340/148‑1 (MQ=255)
                                                                           cgATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGGTCCAACGCGTCGGGCGGGTGCGCtggtg                                                             <  1:227310/148‑1 (MQ=255)
                                                                               cTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAAtgtg                                                                                                                                     <  5:192982/72‑1 (MQ=255)
                                                                               cTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAAtgtg                                                                                                                                     >  6:192982/1‑72 (MQ=255)
                                                                                 ggCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGtt                                                       >  6:358156/1‑148 (MQ=255)
                                                                                  gCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGttt                                                      <  7:71331/148‑1 (MQ=255)
                                                                                     acacAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTg                                                   <  5:358156/148‑1 (MQ=255)
                                                                                           gCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTggcggc                                              >  1:440983/1‑147 (MQ=255)
                                                                                             cACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGGGGTCGTTTCTGGCGGCaaa                                           <  3:256406/148‑1 (MQ=255)
                                                                                               cGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCaaact                                         <  4:42270/148‑3 (MQ=255)
                                                                                                gCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAAtgtg                                                                                                                                     <  6:120862/55‑1 (MQ=255)
                                                                                                gCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAAtgtg                                                                                                                                     >  5:120862/1‑55 (MQ=255)
                                                                                                 cAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCtg                                                                <  6:147361/123‑1 (MQ=255)
                                                                                                 cAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCtg                                                                >  5:147361/1‑123 (MQ=255)
                                                                                                            cTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTg                            >  3:124788/1‑148 (MQ=255)
                                                                                                                        ggCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTgg                 >  3:314638/1‑147 (MQ=255)
                                                                                                                                  gCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGc      >  3:402253/1‑148 (MQ=255)
                                                                                                                                      aCGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGt  >  3:315315/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                       
AATTATGCGAGGTCGTTTTGTAGAAGGTCTGGGTGGCGCGCAATTCGCTGAACGTCTGACTATTGACCGATTGCGCGATCTGGCGACACAAGCCACGCAAACGCGCCACTATACACCAGTGGCGCTCTCTGCCAACGATCCGGCTAATGTGTGGGGAAATCTTCTGCCCTGGCCTGCACATCCGGCAACGCTGGTTCCAACGCGTCGGGCGGGTGCGCTGGTGGTCGTTTCTGGCGGCAAATTGTTACTCTATCTGGCGCAAGGTGGCAAAAAAATGCTGGT  >  NC_000913/1733232‑1733513

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