Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 392,913 C→T 10.6% intergenic (‑12/‑57) insE1 ← / → yaiT IS3 transposase A/pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913392,9130CT10.6% 120.5 / 5.0 47intergenic (‑12/‑57)insE1/yaiTIS3 transposase A/pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein
Reads supporting (aligned to +/- strand):  ref base C (20/22);  new base T (2/3);  total (22/25)
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.87e-01

CTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCAGTTCGATACTGAAGCTTTCCGTACT  >  NC_000913/392827‑393059
                                                                                      |                                                                                                                                                  
cgtCAGGGCTGCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTGTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAAc                                                                                      <  8:129084/147‑1 (MQ=12)
       gCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATAAGATTACAAACTGGACg                                                                               >  5:367539/1‑149 (MQ=12)
        cTTCACTGCGAAATTCAGGCGAATTCTGTTTACTGGGGTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGc                                                                              <  6:330870/149‑1 (MQ=255)
        cTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGc                                                                              <  5:385300/149‑1 (MQ=38)
            aCTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAAt                                                                          <  7:149538/149‑1 (MQ=255)
             cTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATa                                                                         >  1:21658/1‑149 (MQ=255)
              tGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTg                                                                                   >  1:129005/1‑138 (MQ=34)
              tGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTg                                                                                   <  2:129005/138‑1 (MQ=34)
              tGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATAc                                                                        >  7:226914/1‑149 (MQ=255)
                   aaTTCAGGCGAATGCTGTTTACGGGGGTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGGCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCa                                                                   <  7:112544/149‑1 (MQ=255)
                     ttCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCAcc                                                                 <  1:367720/149‑1 (MQ=255)
                         ggCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaa                                                             <  8:399856/149‑1 (MQ=255)
                             aaTGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaactaa                                                         >  6:216085/1‑149 (MQ=255)
                                 cTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAAcgg                                                      <  3:131235/148‑1 (MQ=255)
                                   gTTTAAGGGGTTTTTTACTGGTTGATACTGTTTATGTCATGTGAGTCACCTCTGCCTGAGAGTTTACTCACTTAGCCGCCTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCAAATACTATCCCCGTACCTAAcggcgg                                                   >  2:371893/1‑149 (MQ=17)
                                    tttACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGgt                                                  >  8:116740/1‑149 (MQ=255)
                                        cGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTAt                                              >  5:102064/1‑149 (MQ=255)
                                         ggggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCAATTAGCCCCGTGTCCACAATTGCTGGGAAAGATAAAATAACAACCTGAACGCCAAAACTAACACCGTAAATAACGGCGGTGAAtt                                             >  4:326865/1‑149 (MQ=17)
                                          gggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCAGCGTGTCCACTATTGCTGGGAAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACTGCGGTGTATTg                                            >  5:320574/1‑149 (MQ=255)
                                           ggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTg                                            >  1:351553/1‑148 (MQ=255)
                                            gtgttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTATCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGa                                          <  5:381215/146‑1 (MQ=255)
                                               ttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGAtg                                        <  7:116740/148‑1 (MQ=255)
                                                   aCTGGTTGATACTGTTTTTGTCATGTGAGTCACATCTGCATGAGAGTTTACTAACTTAGCCGCGTGTCCACTATTGCTGGGTAAGCTCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGaa                                    <  6:367539/148‑1 (MQ=21)
                                                    cTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGATAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAAc                                   <  1:128930/148‑1 (MQ=255)
                                                    cTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAAc                                   <  8:226914/148‑1 (MQ=255)
                                                       gTTGATACTGTTTTTGTCATGTGAGTCACCTTTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGtt                                >  2:255086/1‑148 (MQ=255)
                                                       gTTGATACTGTTTTTGTCATGTGAGTCACCTTTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGtt                                >  2:248765/1‑148 (MQ=255)
                                                       gTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGtt                                >  6:79143/1‑148 (MQ=255)
                                                       gTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGt                                 <  6:223078/147‑1 (MQ=255)
                                                       gTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTAATGGATGTGAAc                                   >  2:128930/1‑145 (MQ=255)
                                                        ttGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGa                              >  3:183447/1‑149 (MQ=255)
                                                         tGATACTGTTTTTGTCATGTGAGTCACCTTTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGa                              <  3:326865/148‑1 (MQ=255)
                                                          gATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGAt                             <  6:102064/148‑1 (MQ=255)
                                                              cTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCAGtt                        >  2:429397/1‑149 (MQ=255)
                                                               tGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaactaa                                                         >  3:381430/1‑115 (MQ=255)
                                                               tGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaactaa                                                         <  4:381430/115‑1 (MQ=255)
                                                               tGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGAtg                                        <  1:206366/132‑1 (MQ=255)
                                                               tGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGAtg                                        >  2:206366/1‑132 (MQ=255)
                                                               tGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCAGTTc                       <  1:423549/149‑1 (MQ=255)
                                                                 tttttGTCATGTGAGTCACCTTTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCAGTTCGa                     <  1:255086/149‑1 (MQ=255)
                                                                 tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTATCCGCGTGTCCACTATTGCTGTGTAAGATGAGATTACAACCTGTACGCAAATACTTTCACCGTAACTAACTGCGGTGTATTGGATTTGAACGTTGATCAGTTCGa                     <  1:371893/149‑1 (MQ=21)
                                                                  ttttGTCATGTGAGTCACCTTTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCAGTTCGa                     <  1:248765/148‑1 (MQ=255)
                                                                         aTGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCAGTTCGATACTGAAg             >  2:14733/1‑149 (MQ=255)
                                                                          tgtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCAGTTCGATACTGAAGc            >  7:70599/1‑149 (MQ=255)
                                                                             gAGTCAGCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCa                                                                   <  2:418917/91‑1 (MQ=255)
                                                                              aGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCAGTTCGATACTGAAGCttt         <  5:79143/148‑1 (MQ=255)
                                                                               gTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGTCGGTGTATTGGATATGAACGTTGATCAGTTCGATACAGAAGATTTcc       >  6:128498/1‑149 (MQ=255)
                                                                                     tctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCAGTTCGATACTGAAGCTTTCCGTACt  >  2:45099/1‑148 (MQ=255)
                                                                                      |                                                                                                                                                  
CTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCAGTTCGATACTGAAGCTTTCCGTACT  >  NC_000913/392827‑393059

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