Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,661 T→C N285N (AAT→AAC yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,6610TC88.5% 52.6 / ‑2.5 26N285N (AAT→AACyaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (12/11);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTA  >  NC_000913/46418‑46884
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
gCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAAGCGTTGGTTTGGTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGa                                                                                                                                                                                                                           <  1:128598/251‑1 (MQ=255)
              gCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATTAAGAGCAGCCGCAGCAATCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCCGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCttt                                                                                                                                                                                                             <  2:121423/251‑1 (MQ=255)
               cGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCttt                                                                                                                                                                                                             >  1:97677/1‑250 (MQ=255)
                                                                                 ccGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCATTGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGAGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTgacg                                                                                                                                          >  2:143940/1‑248 (MQ=255)
                                                                                        cgcgATCTTTTTAATCTCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTTATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGTTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGc                                                                                                                                   <  1:173232/251‑1 (MQ=255)
                                                                                              cTGTTTAATCGCCGCCATTTTCCTTTTGTTCTATTTGGTGCCGCCATCTGGACCTGCCAGGTGCTCCCATTGTTAGCCATTTACACCTTTGGCCCGCATATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGGACGTGGTGATTAGCCTGTTCTCTATGCTCGGCTGTCTTCCGCCGGTTCTGTGGTGAAGCACTACAGGACGGAGGCAATTTGTGATTGGCTGGTTTTCCatgat                                                                                                                             >  2:259084/1‑251 (MQ=255)
                                                                                                gTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCatga                                                                                                                              <  1:234108/248‑1 (MQ=255)
                                                                                                    aaTCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgct                                                                                                                       >  1:272165/1‑251 (MQ=255)
                                                                                                     aTCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGGTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGCCGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgctg                                                                                                                      >  2:235666/1‑251 (MQ=255)
                                                                                                        gccgccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgctggcg                                                                                                                   >  1:104774/1‑251 (MQ=255)
                                                                                                                  tCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGAGCCCAATGTTCGCTATTTACACCTTTTGCCCGCAAAGCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTGCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGtt                                                                                                         <  1:381689/250‑1 (MQ=255)
                                                                                                                         gttttgtgtgatgccgcactCTGGACCTGCCGGGTGATCCCACTGGTCGCCATTATCCCCTTGGGCGCGCAAATCGGTGGTTTGTTGGTATTGGGTGTGTGCAAAAACGCGGCGCTGGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACGGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGc                                                                                                  <  2:104774/231‑1 (MQ=255)
                                                                                                                         gTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGcg                                                                                                 >  1:321761/1‑250 (MQ=255)
                                                                                                                              gtttgttGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATc                                                                                             <  1:56748/250‑1 (MQ=255)
                                                                                                                                ttgttgCCGCCATCTGGACCTGCCAGGTGATTCCAATGTTCGCCATTTACACCGTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTGGGGGCTGATccc                                                                                           <  1:261656/250‑1 (MQ=255)
                                                                                                                                ttgttgCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGAGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATccc                                                                                           <  1:143940/250‑1 (MQ=255)
                                                                                                                                          cATCTGGACCTACCAGGTGATCCAAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGATGGTTTGTTGGGATTGGGTGTTGGCAAACACGCGGCGCTGGGGAACGTGGTGCTTAGCCTGGTCTTTATGCTAGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGGCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTtgggggggtcccggatgagggg                                                                                 >  2:295174/1‑232 (MQ=255)
                                                                                                                                           aTCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGGTTTTggggtggtcccgggaaggggggg                                                                                >  2:163277/1‑232 (MQ=255)
                                                                                                                                                                         gCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGgtgtt                                                 >  1:274083/1‑250 (MQ=255)
                                                                                                                                                                                      ttGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGGTCTGGCTGGGAGTGATGGCCTTTGCGGTGTATGCCttttttttt                                    >  1:205337/1‑247 (MQ=255)
                                                                                                                                                                                       tGGCCCGCAAAGCTTGGGGTTTTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGGGGTGATTATCCTTTTCTTTATGCTCGGCTGTAGTCCGCCGATGCTGTGGTTAAACACGGCCTGACGTCGTCCATTGTTGATTGGCAGCTTTGCTATGATGACGGTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTctct                                    <  2:178428/250‑1 (MQ=255)
                                                                                                                                                                                          cccGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCAtgatga                                                                                                                           >  1:348282/1‑161 (MQ=255)
                                                                                                                                                                                          cccGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATGCCGCCGATGCTGTGGTGAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCGTTGCCAggatga                                                                                                                           <  2:348282/161‑1 (MQ=255)
                                                                                                                                                                                            cGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGGTCCCGGATATTGGGGTCTGGCTGGTAGTGGTGGCCTTTGCGGTGTATGCCTTTTTCTTTTGCggg                              >  1:18703/1‑251 (MQ=255)
                                                                                                                                                                                                                        tgTTGGCAAATACGCGGCTCTGGGGCACGTGGTGATTAGCCTGTTCGTTATGCGCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTGCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCt   <  2:272165/249‑1 (MQ=255)
                                                                                                                                                                                                                        ggTTTGCAATAACGCGGGGCTGGGGAACGGGGTGATTTGCCTGTTCTTTATGCTCGGCTGTATTCCTCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTTGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTa  <  2:274083/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                
GCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTA  >  NC_000913/46418‑46884

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: