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,3660CT98.2% 165.5 / ‑4.2 54*75* (TAG→TAAcspDinhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/25);  total (29/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.45e-01

CGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTGGAAT  >  NC_000913/922246‑922499
                                                                                                                        |                                                                                                                                     
cGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCt                                                                                                                     >  2:1307976/1‑139 (MQ=255)
    aTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                            <  2:823015/128‑1 (MQ=255)
    aTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                            >  1:823015/1‑128 (MQ=255)
         gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                            <  2:554277/139‑1 (MQ=255)
                       aGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACt                                                                                              <  1:1307976/139‑1 (MQ=255)
                         tAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg                                                                                            >  1:417270/1‑139 (MQ=255)
                              gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGa                                                                                       >  2:898410/1‑139 (MQ=255)
                              gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTAGACg                                                                                                            >  1:1092878/1‑118 (MQ=255)
                              gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTGTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                            <  2:1092878/118‑1 (MQ=255)
                                         agagTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTg                                                                                                                             <  1:1282312/90‑1 (MQ=255)
                                         agagTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTg                                                                                                                             >  2:1282312/1‑90 (MQ=255)
                                          gagTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                           <  1:898410/139‑1 (MQ=255)
                                                    tGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGAc                                                                                                             <  2:596313/95‑1 (MQ=255)
                                                    tGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGAc                                                                                                             >  1:596313/1‑95 (MQ=255)
                                                    tGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGTGCACAATAACACTgg                                                                                            <  2:1129168/112‑1 (MQ=255)
                                                    tGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg                                                                                            >  1:1129168/1‑112 (MQ=255)
                                                      agaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACa                                                               <  2:417270/139‑1 (MQ=255)
                                                           aaaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCaaa                                                          >  2:672574/1‑139 (MQ=255)
                                                                     ccGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                            <  1:1284748/79‑1 (MQ=255)
                                                                     ccGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                            >  2:1284748/1‑79 (MQ=255)
                                                                     ccGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGAt                                                >  2:1516093/1‑139 (MQ=255)
                                                                              cTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCg                                                                                                               >  2:309293/1‑67 (MQ=255)
                                                                              cTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCg                                                                                                               <  1:309293/67‑1 (MQ=255)
                                                                               tGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                            >  2:1531324/1‑53 (MQ=255)
                                                                               tGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                            <  1:1531324/53‑1 (MQ=255)
                                                                               tGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCtt                                      >  2:833236/1‑139 (MQ=255)
                                                                                  cATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                            >  2:845653/1‑50 (MQ=255)
                                                                                  cATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                            <  1:845653/50‑1 (MQ=255)
                                                                                                   gTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             <  1:477947/80‑1 (MQ=255)
                                                                                                   gTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             >  2:477947/1‑80 (MQ=255)
                                                                                                   gTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAAccatcc                  >  2:713521/1‑139 (MQ=255)
                                                                                                    tACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             >  1:179801/1‑79 (MQ=255)
                                                                                                    tACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             <  2:179801/79‑1 (MQ=255)
                                                                                                      cacaATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACcatccatc               >  1:981672/1‑139 (MQ=255)
                                                                                                       acaATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa                                                                                                  <  1:927785/55‑1 (MQ=255)
                                                                                                       acaATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa                                                                                                  >  2:927785/1‑55 (MQ=255)
                                                                                                       acaATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCt              <  1:713521/139‑1 (MQ=255)
                                                                                                         aaTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggt                                                                    <  2:597254/83‑1 (MQ=255)
                                                                                                         aaTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggt                                                                    >  1:597254/1‑83 (MQ=255)
                                                                                                         aaTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGa            >  2:1365059/1‑139 (MQ=255)
                                                                                                          aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           <  1:672574/139‑1 (MQ=255)
                                                                                                          aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCCTTTGGCCCCTGGTGGACAACAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           <  2:1011710/139‑1 (MQ=255)
                                                                                                              gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCg                                                                                                               <  2:1388160/35‑1 (MQ=255)
                                                                                                              gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCg                                                                                                               >  1:1388160/1‑35 (MQ=255)
                                                                                                              gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                                 <  2:84490/97‑1 (MQ=255)
                                                                                                              gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                                 >  1:84490/1‑97 (MQ=255)
                                                                                                              gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           >  1:612093/1‑135 (MQ=255)
                                                                                                              gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           <  2:612093/135‑1 (MQ=255)
                                                                                                              gaCAGAAGAGCTATGCGACTACCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtggt       >  2:855420/1‑139 (MQ=255)
                                                                                                               aCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGa                                                                                                              >  2:99760/1‑35 (MQ=255)
                                                                                                               aCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGa                                                                                                              <  1:99760/35‑1 (MQ=38)
                                                                                                                 agaagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGa                                                      >  1:217639/1‑89 (MQ=255)
                                                                                                                 agaagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGa                                                      <  2:217639/89‑1 (MQ=255)
                                                                                                                   aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTGGaat  >  1:365918/1‑139 (MQ=255)
                                                                                                                        |                                                                                                                                     
CGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTGGAAT  >  NC_000913/922246‑922499

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