Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,372,254 G→A G107G (GGC→GGT nanT ← sialic acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,372,2540GA93.9% 95.4 / 1.9 33G107G (GGC→GGTnanTsialic acid transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/14);  minor base C (1/1);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

TGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGA  >  NC_000913/3372122‑3372374
                                                                                                                                     |                                                                                                                        
tGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGgc                                                                                                                    >  1:1358277/1‑139 (MQ=255)
               aCTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGcc                                                                                                     <  2:567686/139‑1 (MQ=255)
                   gCTTTGTTACGCAGATGTGTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCg                                                                                                 >  1:895675/1‑139 (MQ=255)
                        gTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCg                                                                                            >  2:113037/1‑139 (MQ=255)
                           aCGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGcc                                                                                         <  1:907357/139‑1 (MQ=255)
                                      ttGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTg                                                                                                                      >  1:1308432/1‑99 (MQ=255)
                                      ttGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTg                                                                                                                      <  2:1308432/99‑1 (MQ=255)
                                           cAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCg                                                                                                           <  2:137806/105‑1 (MQ=255)
                                           cAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCg                                                                                                           >  1:137806/1‑105 (MQ=255)
                                           cAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATg                                                                          >  2:81127/1‑138 (MQ=255)
                                           cAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATg                                                                          <  1:81127/138‑1 (MQ=255)
                                              cTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCaa                                                                                                                >  1:309645/1‑97 (MQ=255)
                                              cTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCaa                                                                                                                <  2:309645/97‑1 (MQ=255)
                                              cTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTgg                                                                      >  1:650482/1‑139 (MQ=255)
                                                  cAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGACAAGAGAACGATGCTGGTGAc                                                                  <  2:895675/139‑1 (MQ=255)
                                                               ggCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAgacga                                                     >  2:523085/1‑139 (MQ=255)
                                                               ggCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAgacga                                                     >  2:1485558/1‑139 (MQ=255)
                                                                  gcTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGc                                                  >  1:753662/1‑139 (MQ=255)
                                                                   cTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGAc                                                                  >  2:542235/1‑123 (MQ=255)
                                                                   cTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGAc                                                                  <  1:542235/123‑1 (MQ=255)
                                                                              ttCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATTCTGGTGACCATTGCCAGACGACGc                                                  <  2:1069683/127‑1 (MQ=255)
                                                                              ttCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGc                                                  >  1:1069683/1‑127 (MQ=255)
                                                                                    cGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATa                                <  2:650482/139‑1 (MQ=255)
                                                                                      ccATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTg                                                                       <  2:1182580/98‑1 (MQ=255)
                                                                                      ccATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTg                                                                       >  1:1182580/1‑98 (MQ=255)
                                                                                      ccATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAgc                              >  1:70044/1‑139 (MQ=255)
                                                                                         tCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGcc                           <  2:753662/139‑1 (MQ=255)
                                                                                                         cAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGccgccg           >  1:648435/1‑139 (MQ=255)
                                                                                                                 cGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaagagaa                                                                               >  1:835631/1‑63 (MQ=255)
                                                                                                                 cGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaagagaa                                                                               <  2:835631/63‑1 (MQ=255)
                                                                                                                  gATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAg                                            <  2:1384525/97‑1 (MQ=255)
                                                                                                                  gATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAg                                            >  1:1384525/1‑97 (MQ=255)
                                                                                                                  gATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCga  <  2:1358277/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                        
TGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGA  >  NC_000913/3372122‑3372374

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