Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 922,366 C→T *75* (TAG→TAA cspD ← inhibitor of DNA replication, cold shock protein homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913922,3660CT97.8% 138.5 / ‑6.9 46*75* (TAG→TAAcspDinhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/20);  minor base A (1/0);  total (26/20)
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.99e-01

GCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTGGAAT  >  NC_000913/922243‑922499
                                                                                                                           |                                                                                                                                     
gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCt                                                                                                                        >  1:1167330/1‑139 (MQ=255)
              aCGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                          <  2:10871/139‑1 (MQ=255)
                            tAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg                                                                                            >  2:991791/1‑139 (MQ=255)
                                 gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGAc                                                                                                             <  2:598292/117‑1 (MQ=255)
                                 gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGAc                                                                                                             >  1:598292/1‑117 (MQ=255)
                                      ggCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGcc                                                                                  >  1:315148/1‑139 (MQ=255)
                                              agTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcc                                                                          >  2:333056/1‑139 (MQ=255)
                                               gTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGATATGCTACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCtttttaccc                                                                         >  1:583390/1‑134 (MQ=255)
                                                         agaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggt                                                                    >  1:1192547/1‑134 (MQ=255)
                                                         agaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggt                                                                    <  2:1192547/134‑1 (MQ=255)
                                                         agaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACa                                                               >  1:1079460/1‑139 (MQ=255)
                                                                      aGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGc                                                                                                                         >  2:575382/1‑68 (MQ=255)
                                                                      aGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGc                                                                                                                         <  1:575382/68‑1 (MQ=255)
                                                                                  tGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCtt                                      >  1:1281078/1‑139 (MQ=255)
                                                                                   ggCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGCGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGa                                                      <  1:689123/122‑1 (MQ=255)
                                                                                   ggCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGa                                                      >  2:689123/1‑122 (MQ=255)
                                                                                              tttAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTAGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGt                          <  2:980646/139‑1 (MQ=255)
                                                                                                             aggagaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGAAATCAAACTGAAGGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           >  2:777038/3‑139 (MQ=255)
                                                                                                             aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGc                                                                                                         <  1:522449/45‑1 (MQ=255)
                                                                                                             aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGc                                                                                                         >  2:522449/1‑45 (MQ=255)
                                                                                                             aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           <  1:784955/139‑1 (MQ=255)
                                                                                                             aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTAGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           >  2:784955/1‑139 (MQ=255)
                                                                                                              tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTc                             <  1:556345/120‑1 (MQ=255)
                                                                                                              tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTc                             >  2:556345/1‑120 (MQ=255)
                                                                                                               gagaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtg         <  1:903941/139‑1 (MQ=255)
                                                                                                               gagaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtg         <  2:1281078/139‑1 (MQ=255)
                                                                                                                 gaCAGCAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACaa                                                                                                     >  1:890118/1‑45 (MQ=255)
                                                                                                                 gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                          >  2:289730/1‑40 (MQ=255)
                                                                                                                 gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                          >  2:724017/1‑40 (MQ=255)
                                                                                                                 gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                          <  1:289730/40‑1 (MQ=39)
                                                                                                                 gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                          <  1:724017/40‑1 (MQ=39)
                                                                                                                 gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACaa                                                                                                     <  2:890118/45‑1 (MQ=255)
                                                                                                                 gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtgg                                                                  <  2:465393/80‑1 (MQ=255)
                                                                                                                 gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtgg                                                                  >  1:465393/1‑80 (MQ=255)
                                                                                                                 gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGt                               <  2:340421/115‑1 (MQ=255)
                                                                                                                 gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGt                               >  1:340421/1‑115 (MQ=255)
                                                                                                                 gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTaa                        <  1:527159/122‑1 (MQ=255)
                                                                                                                 gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTaa                        >  2:527159/1‑122 (MQ=255)
                                                                                                                      aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                                 <  2:1149142/92‑1 (MQ=255)
                                                                                                                      aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                                 >  1:1149142/1‑92 (MQ=255)
                                                                                                                      aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           <  2:1133500/130‑1 (MQ=255)
                                                                                                                      aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           >  1:1133500/1‑130 (MQ=255)
                                                                                                                      aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           >  2:909572/1‑130 (MQ=255)
                                                                                                                      aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           <  1:909572/130‑1 (MQ=255)
                                                                                                                      aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTGGaat  >  2:616946/1‑139 (MQ=255)
                                                                                                                        gagTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAAcca                     >  1:251409/1‑118 (MQ=255)
                                                                                                                           |                                                                                                                                     
GCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTGGAAT  >  NC_000913/922243‑922499

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