Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,780,371 C→T A252V (GCC→GTC)  ydiQ → putative electron transfer flavoprotein subunit YdiQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,780,3710CT97.3% 128.4 / ‑5.8 37A252V (GCC→GTC) ydiQputative electron transfer flavoprotein subunit YdiQ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/15);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACAGAACGTAAGCACATCATTATCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGCCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAA  >  NC_000913/1780293‑1780445
                                                                              |                                                                          
aaaCAGAACGTAAGCACATCATTATCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGtcc                                                                          >  1:317788/1‑81 (MQ=255)
 aaCAGAACGTAAGCACATCATTATCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCct                                                                         >  1:317791/1‑81 (MQ=255)
      aaCGTAAGCACATCATTATCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCctgaact                                                                    >  1:317775/1‑81 (MQ=255)
       aCGTAAGCACATCATTATCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCctgaact                                                                    >  2:317778/1‑80 (MQ=255)
        cGTAAGCACATCATTATCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCctgaactg                                                                   >  1:317771/1‑80 (MQ=255)
               acaTCATTATCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTAc                                                           <  1:317741/81‑1 (MQ=255)
                 atcatTATCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACgg                                                         <  1:317742/81‑1 (MQ=255)
                       aTCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGagaag                                                    >  2:317794/1‑80 (MQ=255)
                       aTCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGagaag                                                    <  1:317743/80‑1 (MQ=255)
                       aTCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGagaag                                                    <  1:317744/80‑1 (MQ=255)
                        tCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGagaaga                                                   <  1:317745/80‑1 (MQ=255)
                           aTAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAAc                                                >  2:317776/1‑80 (MQ=255)
                            tAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGa                                              >  2:317783/1‑81 (MQ=255)
                            tAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGa                                              <  1:317746/81‑1 (MQ=255)
                             aaCGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGAt                                             >  1:317787/1‑81 (MQ=255)
                              aCGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGAt                                             <  2:317747/80‑1 (MQ=255)
                                  ttCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGt                                         >  1:317777/1‑80 (MQ=255)
                                  ttCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGt                                         >  2:317774/1‑80 (MQ=255)
                                     gCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAAt                                     >  2:317784/1‑81 (MQ=255)
                                             ccATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGc                             >  2:317780/1‑81 (MQ=255)
                                              cATTGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCg                            >  1:317781/1‑81 (MQ=255)
                                                 tGCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCt                         >  2:317789/1‑81 (MQ=255)
                                                  gCCGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTg                        >  2:317786/1‑81 (MQ=255)
                                                    cGAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTggg                      <  1:317748/81‑1 (MQ=255)
                                                     gAGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGt                     >  2:317798/1‑81 (MQ=255)
                                                      aGCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTc                    >  1:317792/1‑81 (MQ=255)
                                                       gCTGGCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCt                   <  1:317749/81‑1 (MQ=255)
                                                          ggCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGTGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCtt                  >  1:317790/1‑79 (MQ=255)
                                                           gCGGAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCttt                 <  1:317750/79‑1 (MQ=255)
                                                             ggAACATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCg             <  1:317751/81‑1 (MQ=255)
                                                               aaCATCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGAt           <  1:317752/81‑1 (MQ=255)
                                                                   tCTGAAGAAAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATc       >  1:317793/1‑81 (MQ=255)
                                                                       aagaaAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCt     <  1:317756/79‑1 (MQ=255)
                                                                       aagaaAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCt     <  1:317755/79‑1 (MQ=255)
                                                                       aagaaAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCt     <  1:317753/79‑1 (MQ=255)
                                                                       aagaaAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTg    >  2:317785/1‑80 (MQ=255)
                                                                       aagaaAGCCCTGAACTGTAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATc       <  2:317754/77‑1 (MQ=255)
                                                                        agaaAGTCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGaa  >  1:317800/1‑81 (MQ=255)
                                                                              |                                                                          
AAACAGAACGTAAGCACATCATTATCGATAACGATTCGCCGGAGGCCATTGCCGAGCTGGCGGAACATCTGAAGAAAGCCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAA  >  NC_000913/1780293‑1780445

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