Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,750,766 A→T 9.1% K42M (AAG→ATG) ‡ PP_1551 → replication protein O

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,750,7660AT 9.1% 118.8 / 3.6 45K42M (AAG→ATG) ‡PP_1551replication protein O
Reads supporting (aligned to +/- strand):  ref base A (22/18);  new base T (2/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

CAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGT  >  NC_002947/1750620‑1750914
                                                                                                                                                  |                                                                                                                                                    
cAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCaa                                                                                                                                                      <  6:135480/147‑1 (MQ=255)
  gCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGAGGAAGATTGTGCTGTACGTGGCGAAGGcc                                                                                                                                                  >  7:149370/1‑149 (MQ=255)
   cATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGcca                                                                                                                                                 <  6:337166/149‑1 (MQ=255)
        gAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCAt                                                                                                                                             >  4:162683/1‑148 (MQ=255)
          aCGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCTTGGccacgtacag                                                                                                                                          <  2:126846/149‑7 (MQ=255)
          aCGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCTTGGccacgtacag                                                                                                                                          <  2:209812/149‑7 (MQ=255)
                tAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATATCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCAACAACTTCgg                                                                                                                                    >  3:232675/1‑149 (MQ=255)
                   aGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGc                                                                                                                                 <  8:105140/149‑1 (MQ=255)
                       tGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGgc                                                                                                                             >  3:147222/1‑149 (MQ=255)
                          cAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGgcgc                                                                                                                           >  2:60301/1‑148 (MQ=255)
                             aCATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCaa                                                                                                                       <  3:125422/149‑1 (MQ=255)
                                   tCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATc                                                                                                                 <  1:60301/149‑1 (MQ=255)
                                        cTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGc                                                                                                            <  1:67593/149‑1 (MQ=255)
                                            aCAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTAc                                                                                                         <  6:236540/148‑1 (MQ=255)
                                              aaGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGa                                                                                                      >  6:298398/1‑149 (MQ=255)
                                                gTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACa                                                                                                    >  6:20829/1‑149 (MQ=255)
                                                   cAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATcg                                                                                                 >  8:172628/1‑149 (MQ=255)
                                                        ggTTCACCCGGATGGACAAACAGCTCATGGATGGCCTGCTGGCTATCCATCTCCCGGCCCGGGAGATCAAGATTGTCCTGTACGTGGCCAAGGCCCCCATCAACTTCGGTGCGCGCGCTCAACGCATCCCGCCTCACGCCCTCCCGaaa                                                                                            >  8:148879/1‑149 (MQ=255)
                                                         gTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTAcc                                                                                                        <  2:193651/136‑1 (MQ=255)
                                                         gTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTAcc                                                                                                        >  1:193651/1‑136 (MQ=255)
                                                               ccGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGCCATCGCGAAAGCCATc                                                                                      >  5:52248/1‑148 (MQ=255)
                                                                     ggACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCt                                                                                >  4:308693/1‑148 (MQ=255)
                                                                      gACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTc                                                                               <  1:201148/148‑1 (MQ=255)
                                                                         aaCCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCaa                                                                                                                       <  4:81110/105‑1 (MQ=255)
                                                                          aCCAGCTCATGGATGGCCTTCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCTAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCt                                                                          <  8:149370/149‑1 (MQ=255)
                                                                            cAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCGACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGa                                                                        <  5:20829/149‑1 (MQ=255)
                                                                                  atggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgt                                                                  >  6:171633/1‑149 (MQ=255)
                                                                                  atggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACACATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgt                                                                  <  4:147222/149‑1 (MQ=255)
                                                                                     gatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCACGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgtgt                                                                <  7:172628/148‑1 (MQ=255)
                                                                                          ccTGCTTGCTATCTATCTCCCGGCCCTGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAgg                                                          <  6:52248/149‑1 (MQ=255)
                                                                                           ctgctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCt                                                                          >  2:144366/1‑132 (MQ=255)
                                                                                           ctgctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCt                                                                          <  1:144366/132‑1 (MQ=255)
                                                                                           ctgctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGc                                                         >  5:131191/1‑149 (MQ=255)
                                                                                             gctgGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCgg                                                       >  7:234235/1‑149 (MQ=255)
                                                                                                 gCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTc                                                   >  2:292213/1‑149 (MQ=255)
                                                                                                  cTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTcc                                                  <  5:171633/149‑1 (MQ=255)
                                                                                                  cTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTc                                                   >  5:190375/1‑148 (MQ=255)
                                                                                                       gATCTCCCGGCCCGGGAGATGAAGATTGTGATGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAGGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCct                                             >  6:297597/1‑149 (MQ=255)
                                                                                                             ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGc                                                         >  2:315032/1‑131 (MQ=255)
                                                                                                             ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGc                                                         <  1:315032/131‑1 (MQ=255)
                                                                                                             ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgct                                          <  5:298398/146‑1 (MQ=255)
                                                                                                                  ccggGAGATGATGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGACAGCCATCCACGCTCACCCTGCCACCGTGTCCAAGGCGGGTTCCAGCCTGCTGCGTCGcc                                  <  7:148879/149‑1 (MQ=255)
                                                                                                                                     ttgatggTGGCCTTGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAgg               >  4:156710/7‑149 (MQ=255)
                                                                                                                                     ttgatggTGGCCTTGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACAACGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAgg               >  8:237317/7‑149 (MQ=255)
                                                                                                                                              gCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCg      >  5:78708/1‑149 (MQ=255)
                                                                                                                                                cAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCggg    >  7:250274/1‑149 (MQ=255)
                                                                                                                                                  aGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGt  <  3:162683/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGT  >  NC_002947/1750620‑1750914

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: ‑