Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 367,644 G→T intergenic (‑96/‑20) nanA ← / → nagC N‑acetylneuraminate lyase/possible transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730367,6440GT97.4% 125.1 / ‑2.7 39intergenic (‑96/‑20)nanA/nagCN‑acetylneuraminate lyase/possible transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (16/22);  total (16/23)
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.17e-01

AGTAACGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGT  >  CP000730/367511‑367778
                                                                                                                                     |                                                                                                                                      
aGTAACGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATAtt                                                                                                                                 >  2:113084/1‑141 (MQ=255)
aGTAACGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATAtt                                                                                                                                 >  2:4828/1‑141 (MQ=255)
aGTAACGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATAtt                                                                                                                                 >  1:15217/1‑141 (MQ=255)
 gTAACGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATAttt                                                                                                                                >  2:288441/1‑141 (MQ=255)
             ataAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTAtaattaa                                                                                                                    >  2:55964/1‑141 (MQ=255)
              taAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTAtaattaat                                                                                                                   <  1:191937/141‑1 (MQ=255)
                aaacTTTTAAATCTTTATTGATTATAAAGTCACTGCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTATACGCTATCCTAAAATGTAACAATATCTTTTTGAAAATTGTATTGATATTTATGTATAATTAATGa                                                                                                                 <  2:77430/138‑1 (MQ=255)
                 aCCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCACTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGaa                                                                                                                <  2:283981/141‑1 (MQ=255)
                     tttAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATg                                                                                                            <  2:62279/141‑1 (MQ=255)
                           tCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTc                                                                                                      <  1:109244/141‑1 (MQ=255)
                             tttGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCAT‑AAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCaaa                                                                                                   >  2:592/1‑141 (MQ=255)
                                   cATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCaa                                                                                              <  2:140928/141‑1 (MQ=255)
                                   cATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCaa                                                                                              <  1:6411/141‑1 (MQ=255)
                                        tGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGaa                                                                                         <  1:15044/141‑1 (MQ=255)
                                           aGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATTTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAAATAATGAAAATGTTTTTCAAAATCAATAGAAATg                                                                                      >  1:8856/1‑141 (MQ=255)
                                            gTCCCTCCAAAAAGCATTTGATAATATCGTTTACAGCTATAATTGTAAACGCTTACATTAATTGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATgg                                                                                     >  2:85669/1‑141 (MQ=255)
                                             tCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGga                                                                                    >  1:279547/1‑141 (MQ=255)
                                                 tCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGttttt                                                                                                       >  1:216164/1‑118 (MQ=255)
                                                 tCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGttttt                                                                                                       <  2:216164/118‑1 (MQ=255)
                                                                ataataTAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACa                                                                 <  1:5693/141‑1 (MQ=255)
                                                                 taataTAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACAt                                                                <  2:103389/141‑1 (MQ=255)
                                                                      taGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCaa                                                           <  2:184401/141‑1 (MQ=255)
                                                                        gTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATc                                                         <  1:4828/141‑1 (MQ=255)
                                                                         tttACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCg                                                        <  1:224242/141‑1 (MQ=255)
                                                                                  ataATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAg                                               <  2:264022/141‑1 (MQ=255)
                                                                                         tAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTTTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCa                                        >  2:19768/1‑141 (MQ=255)
                                                                                               cTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTaa                                  >  2:117946/1‑141 (MQ=255)
                                                                                                  tCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATc                               <  2:61760/141‑1 (MQ=255)
                                                                                                    aTAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCgg                             <  1:55964/141‑1 (MQ=255)
                                                                                                    aTAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCgg                             >  2:227785/1‑141 (MQ=255)
                                                                                                     tAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGaaa                                                                                        >  2:19957/1‑81 (MQ=255)
                                                                                                     tAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGaaa                                                                                        <  1:19957/81‑1 (MQ=255)
                                                                                                                 atatCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATaa                <  1:227785/141‑1 (MQ=255)
                                                                                                                      tttttGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAAt           <  2:170413/141‑1 (MQ=255)
                                                                                                                      tttttAAAAAGTGTATTGATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGGGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTAAAATTAAATCGGCAGTTTTTGCTAAGCAAt           <  2:188377/141‑1 (MQ=255)
                                                                                                                          tGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGaa       <  2:192229/141‑1 (MQ=255)
                                                                                                                           gAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAat      <  1:78144/141‑1 (MQ=255)
                                                                                                                               aTTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTGAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCAATCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGt  >  2:181152/1‑141 (MQ=255)
                                                                                                                               aTTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGt  >  1:109509/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                      
AGTAACGCTGCATATAAACCTTTTAAATCTTTGTTCATTATGAAGTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGGCACTCAAATTAAATCGGCAGTTATTGATAAGCAATTGAATATGT  >  CP000730/367511‑367778

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