Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,227,340 (C)6→7 pseudogene (1270/1287 nt) mdtQ ← pseudogene; putative channel/filament proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,227,3341.C96.7% 103.1 / ‑2.9 30pseudogene (1276/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (15/14);  total (16/14)
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.44e-01

TCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGACGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGC  >  NC_000913/2227188‑2227478
                                                                                                                                                   |                                                                                                                                                
tCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTAcccc                                                                                                                                              >  2:91055/1‑151 (MQ=255)
    aCCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTAcccccccg                                                                                                                                          >  1:184437/1‑149 (MQ=255)
     ccccTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTAcccccccgc                                                                                                                                         >  2:451625/1‑148 (MQ=255)
          gctgcGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAAc                                                                                                                                     >  1:71531/1‑151 (MQ=255)
            tgcGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGc                                                                                                                                   >  1:403593/1‑151 (MQ=255)
                 cAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCAcc                                                                                                                                >  1:628880/1‑149 (MQ=255)
                   gTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTg                                                                                                                            <  1:625749/151‑1 (MQ=255)
                               aCTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAg                                                                                                                <  2:343192/151‑1 (MQ=255)
                                    ttGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAgga                                                                                                              <  1:257665/148‑2 (MQ=255)
                                           tCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCAc                                                                                                    >  1:433089/1‑151 (MQ=255)
                                            ccccGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCAc                                                                                                    <  1:407277/150‑1 (MQ=255)
                                              ccGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGc                                                                                                 <  1:142362/151‑1 (MQ=255)
                                              ccGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGc                                                                                                 >  1:93777/1‑151 (MQ=255)
                                                           ggTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACa                                                                                    <  2:306078/151‑1 (MQ=255)
                                                                taatCGCAAACTGTGCCATCGGGTTATCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGAGTTTGTACCCCCCACCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCAt                                                                              >  2:293528/1‑151 (MQ=255)
                                                                   tCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTg                                                                            >  1:619176/1‑151 (MQ=255)
                                                                    cGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTg                                                                            <  2:104941/150‑1 (MQ=255)
                                                                         ccTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGc                                                                      >  1:526901/1‑151 (MQ=255)
                                                                                cATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGc                                                                >  2:485876/1‑150 (MQ=255)
                                                                                    gAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAgcgc                                                           >  1:378618/1‑151 (MQ=255)
                                                                                    gAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAgcgc                                                           >  2:53000/1‑151 (MQ=255)
                                                                                           ccACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCgggg                                                       <  1:4952/148‑1 (MQ=255)
                                                                                           ccACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGAtt                                                    <  2:85156/151‑1 (MQ=255)
                                                                                                tAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGc                                               >  1:55348/1‑151 (MQ=255)
                                                                                                   cGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTc                                            <  2:503371/151‑1 (MQ=255)
                                                                                                        aaTAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTg                                        >  1:593798/1‑150 (MQ=255)
                                                                                                                           aGGCGCATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGACGCGGGAACCAGCAatgat                    <  1:295625/151‑1 (MQ=255)
                                                                                                                               gcATATCATCAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGACGCGGGAACCAGCAATGATGccc                <  2:218648/151‑1 (MQ=255)
                                                                                                                                   atcatcAGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAAGCAACCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGATCACGCAGCGCGGGGATTCTGGCTTCGCTGACGCGGGAACCAGCAATGATGCCCGCGt            >  2:359886/1‑151 (MQ=255)
                                                                                                                                     catcaGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGACGCGGGAACCAGCAATGATGCCCGCGtt           <  2:431639/150‑1 (MQ=255)
                                                                                                                                     catcaGCGTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGACGCGGGAACCAGCAATGATGCCCGCGtt           >  1:68066/1‑150 (MQ=255)
                                                                                                                                            gTTTGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGACGCGGGAACCAGCAATGATGCCCGCGTTAAAgcgcg   <  1:464443/151‑1 (MQ=255)
                                                                                                                                             tttGTACCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGACGCGGGAACCAGCAATGATGCCCGCGTTAAAgcgcgc  >  1:470634/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCAAACCCCTGCTGCGCCAGTAATAACGCGCACTCTTTGCCGATCCCCGAATCGGAGGCGGTAATAATCGCAACCTGTGCCATCGAGTTCTCCACTTAACGCTGAATAAACGTTAAGTATAGAAGGCGCATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGACGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGC  >  NC_000913/2227188‑2227478

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