Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,793,377 Δ1 bp coding (182/750 nt) btuD ← vitamin B12 ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,793,3770T.92.9% 89.8 / 1.0 28coding (182/750 nt)btuDvitamin B12 ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (10/16);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

ACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATC  >  NC_000913/1793169‑1793617
                                                                                                                                                                                                                |                                                                                                                                                                                                                                                
aCCGCCGGAAAGTTGATTGGTGCTACGGCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCTTCTGCTGGTGTGAAAGATAAGCGCGGCGCAGCGCGAGGTCTGGAGCGGACCATGCTTCCAGTGGTGCCCCGGGAACTGAATGCTTCCCTTACCGCTGGTCATTCCgcc                                                                                                                                                                                                       <  1:415047/251‑3 (MQ=255)
              gATTGGTGCTACGTCCGAGATTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTAAGTAGCTCGGTACGCGTTTTATAGTGCTTATGCAGTGTCCGGTAGTGCCAGACCGGCGTTGCACACGGCGGCGTCTGCTGTTGTGAAAGATATGCGCGATGCCGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCCTTCGCGCCAGt                                                                                                                                                                                          >  1:431145/1‑250 (MQ=255)
                         cGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACTCGTTTTATCGTGCTGATGCCGTGTCAGGTCGTGCCAGACCTGCGTTGCAAACGGCGTCGTCTGCTGTTGTGCAAGATCGGCGCGATGCCGCGCGCGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGACCTGAATGCTTCCCTTACCGCTGTTCATTCCGGCCATTCGCGCCAGtcagttcctctt                                                                                                                                                                               >  1:75143/1‑239 (MQ=255)
                           nccGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGGGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTAATCTTAcc                                                                                                                                                                            >  1:45862/2‑251 (MQ=255)
                                         caAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTTGTGCCCCGCTAGCTGAATGCTTCCCTTACCGGTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCggn                                                                                                                                                               <  1:93419/250‑2 (MQ=255)
                                            gagCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGTCTTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCAGGCTTCCAGTGGTGCCCCGCGAACTTAATGCTTCCCTTACCGCTTGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCcca                                                                                                                                                           <  1:397061/251‑1 (MQ=255)
                                                   cgcCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCATCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGTCAGTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGc                                                                                                                                                               <  1:437939/240‑1 (MQ=255)
                                                     ccccTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCccac                                                                                                                                                          <  1:363027/243‑1 (MQ=255)
                                                      cccTGCGACATCATTCAGTAGTGCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATTCTTACAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGTTCATGCCGGCCAGTCGCGCCAGTAAGCTACTCGTACCCGCGCCATTCGGCCCCACCAGGTGCAg                                                                                                                                                 <  1:286269/251‑1 (MQ=255)
                                                               atcatTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGGAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGGTTTGTTGCGGACCTTGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCAGTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTcccn                                                                                                                                        <  1:52864/251‑2 (MQ=255)
                                                                                          aTCGTGCTGATGCAGGTTCAGGTGGTGCCAGGCCGGCGTTGCAAGCGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGGGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCTATTCGGCCACACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCggac                                                                                                              <  1:425416/250‑3 (MQ=255)
                                                                                                                tAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAgn                                                                                                                                                      <  1:279808/188‑2 (MQ=255)
                                                                                                                tAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATNCCGNNNNNTcnn                                                                                                                                                                                                >  2:279808/1‑144 (MQ=255)
                                                                                                                  gTGCCAGTCCGGCGTAGCCAACGGCGGCGTCGTCTGTTGTGAAAGAGAGGCTCGTTGCAGCGCGAGTTTTGTTGCGGACCAGGTTTCCAGTGGTTCACGGCGGACTGAATGCGGCCCTTACCGTTGGTCATTCCGGCCATGCGCGCCAGTATGTTACTCTTAcc                                                                                                                                                                            <  1:433445/164‑1 (MQ=255)
                                                                                                                                  gCAAACGGCGGCGTCTGCTGTTGTGGAAGATAGGCGCGATCCAGCGCGAGTTTTGTTGCGGACCCTGCTTCCCGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGGAATTCCCGCCCTTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGGGCAGGATCTCCCCAGCCCGAACCTCGCCAGACAACCGCCCCAGCCGGGTCGAATTCGCCACATCTTGTAACTTCATcacc                                                                     >  1:407421/1‑250 (MQ=255)
                                                                                                                                   cAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTAcc                                                                                                                                                                                                                      <  1:299203/105‑1 (MQ=255)
                                                                                                                                                     gttgtGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAAGAGACATTATTTTGcc                                                    <  1:357080/249‑1 (MQ=255)
                                                                                                                                                          gAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCCTGCTTCCAGTGGTGACCCGCGAACTGAATGCTTACCTTACCGATGGCCATTCCGGCCATTCGAGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTg                                                                                                                                                    >  1:317495/1‑148 (MQ=255)
                                                                                                                                                          gAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCNCTANCTGANTGNTTCCCTTACCGCTGGNCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATNNGNNNCNACCAGGTg                                                                                                                                                    <  2:317495/148‑1 (MQ=255)
                                                                                                                                                           aaaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAgg                                                                                                                                                      <  1:178221/145‑1 (MQ=255)
                                                                                                                                                                  ggCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGCTCTGAATGCCTCCCTTACCGCTGGTAATTCCGGCCATTCGCGCCAGTAACGACCTCTTAcc                                                                                                                                                                            >  1:70827/1‑116 (MQ=255)
                                                                                                                                                                      cgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTNc                                                     <  1:347978/231‑1 (MQ=255)
                                                                                                                                                                      cgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAAAGCCAGTTTAATACttt                                 >  1:169825/1‑251 (MQ=255)
                                                                                                                                                                      cgATGCAGCGCGAGCTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGATACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGa                                                                                                 >  1:375707/1‑187 (MQ=255)
                                                                                                                                                                                            cGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCGCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCgn           <  1:90825/251‑2 (MQ=255)
                                                                                                                                                                                             ggACCNNGCNNNCANNGGTGCCCCGAGAACTGAATGCTTNCCTTCCCGCTGGTCATTCCGGCCCTTCGCNCNANTAAGNNACTNTTACCTGCGCNATTCGGCCCCACCAgg                                                                                                                                                      >  2:322573/1‑111 (MQ=255)
                                                                                                                                                                                             ggACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGTTCATTCCGGCCATTCGCGCCAGTAAGGTACTCGTACCCGCGCCATTCGGCCCCACCAgg                                                                                                                                                      <  1:322573/111‑1 (MQ=255)
                                                                                                                                                                                                       tCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATc  <  1:356995/249‑1 (MQ=255)
                                                                                                                                                                                                              ggtGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCNNCCAGGtnn                                                                                                                                                   >  2:12734/3‑96 (MQ=255)
                                                                                                                                                                                                                tGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGGAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCg            <  1:66441/231‑1 (MQ=255)
                                                                                                                                                                                                                |                                                                                                                                                                                                                                                
ACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATC  >  NC_000913/1793169‑1793617

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