Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,707,029 A→G 19.9% V37V (GTA→GTG rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,707,0290AG19.9% 109.5 / 19.1 55V37V (GTA→GTGrsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (23/21);  new base G (7/4);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCC  >  NC_000913/1706933‑1707127
                                                                                                |                                                                                                  
ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTAccc                                                                                                 >  2:629318/1‑100 (MQ=255)
ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTAccc                                                                                                 <  1:964276/100‑1 (MQ=255)
  ctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCt                                                                                               <  1:1694010/100‑1 (MQ=255)
     cATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGgc                                                                                            <  2:747045/100‑1 (MQ=255)
      aTTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGgcg                                                                                           >  1:597777/1‑100 (MQ=255)
        tCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCa                                                                                         <  2:654714/100‑1 (MQ=255)
         cAGAAAAAATAAAATCTGCGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGgcgccg                                                                                        >  2:1206088/1‑98 (MQ=255)
          aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGc                                                                                       >  1:340199/1‑100 (MQ=255)
          aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGc                                                                                       >  1:639752/1‑100 (MQ=255)
           gAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAAAGGTAAACCCCTCCGAAAGGTGCCCCTGCCGCAGCg                                                                                      >  2:1236650/1‑100 (MQ=255)
                aataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTg                                                                                 <  1:851954/100‑1 (MQ=255)
                aataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTg                                                                                 <  2:1031696/100‑1 (MQ=255)
                 ataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGt                                                                                <  1:22001/100‑1 (MQ=255)
                  taaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACACCTGGAGCAGCGTTTTGt                                                                                >  2:841475/1‑99 (MQ=255)
                  caaaatCTGGGATCTCAACGGCGGCATACATCCACCGGAGATGAAAACCCAGTCCAACAGCACACCCATGCGCCAGGTGCCCCTGGCGTAGc                                                                                       >  2:396168/2‑92 (MQ=255)
                       tCTGGGTTTTCAACTGTGGCATCCATCCACCGTAGATGAAAACCCAGTCCAACTGTACACTCCTGCGCCAGGTACCCTTTGCGCAGCGTTTTGTTATTcc                                                                          <  2:360631/100‑1 (MQ=255)
                       tCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTcc                                                                          <  1:1392686/100‑1 (MQ=255)
                         tGGGATTTCAACGGCGGCATCCATCAACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGGCAGATACCCCTGGCGCAGCGTTTTGTTATTCCAc                                                                        >  2:2010634/1‑100 (MQ=255)
                          gggATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACt                                                                       <  1:1012663/100‑1 (MQ=255)
                           ggATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTg                                                                      <  1:1184233/100‑1 (MQ=255)
                              tttCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGaaa                                                                   >  1:1156303/1‑100 (MQ=255)
                                tCAACGGCGGCATTCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACa                                                                 <  1:731735/100‑1 (MQ=255)
                                tCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACa                                                                 >  1:1955580/1‑100 (MQ=255)
                                 cAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACa                                                                 >  1:269/1‑99 (MQ=255)
                                  aaCGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGc                                                               <  1:1585400/100‑1 (MQ=255)
                                    cggcggCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCat                                                             <  2:1586797/100‑1 (MQ=255)
                                         gtatccatccACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTgg                                                        <  1:1920261/98‑1 (MQ=255)
                                         gcATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTgg                                                        <  2:392319/100‑1 (MQ=255)
                                                   aCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAgg                                               <  2:880947/99‑1 (MQ=255)
                                                    ccGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTg                                             >  1:486865/1‑100 (MQ=255)
                                                    ccGGAGATCAAAACCCAGTCCAACGGTACACCCCTGCGCAAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGt                                              <  2:550414/99‑1 (MQ=255)
                                                      ggAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAg                                           <  1:2232612/100‑1 (MQ=255)
                                                       gagaTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGt                                          <  1:1937066/100‑1 (MQ=255)
                                                           tGAAACCCCAGTCCAACGGTAAACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAAAAGCATAGCGGCGCTGAagttgagttg                                        >  2:458305/1‑98 (MQ=255)
                                                           tGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                                      <  2:2364143/100‑1 (MQ=255)
                                                           tGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAACCATCATATTGGCGCTGAAGGTGAGTtgtg                                      >  1:384260/1‑100 (MQ=255)
                                                             aaaaCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCg                                    >  2:591741/1‑100 (MQ=255)
                                                              aaaCCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGt                                   >  1:953750/1‑100 (MQ=255)
                                                               aaCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCg                                    >  2:641333/1‑98 (MQ=255)
                                                                aCCCAGTCCAACGGTACACCCCTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTa                                 >  2:747722/1‑100 (MQ=255)
                                                                    aGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTTTTTTTTTTTTTAAATAGAATATTGGCGCTGAAGGTGAGTTGTGGGTTAgtgt                             >  2:1377391/1‑100 (MQ=255)
                                                                     gTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTc                            <  1:1003793/100‑1 (MQ=255)
                                                                        cAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTCTTATTCCACTGCACCAGCATACTGGCGCCGAAGCTGCCCTGTGCGTTACCGTc                            >  2:644444/1‑97 (MQ=255)
                                                                          aCGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGa                       >  1:2394570/1‑100 (MQ=255)
                                                                               acacCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAg                  >  2:552233/1‑100 (MQ=255)
                                                                                  ccccTGCGCCAGGTGCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAt               <  1:1674493/100‑1 (MQ=255)
                                                                                  ccccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAt               >  2:1692250/1‑100 (MQ=255)
                                                                                   cgcTGCGCCAGGGGTACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAtt              >  2:1863565/3‑100 (MQ=38)
                                                                                   cgcTGCGCCAGGGGCACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAtt              >  1:21256/3‑100 (MQ=255)
                                                                                   cccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGACGGCGATAAAGTAtt              >  1:646022/1‑100 (MQ=255)
                                                                                    ccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTg             >  2:1746844/1‑100 (MQ=255)
                                                                                       gcgcCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTgcgc          >  2:1380023/1‑100 (MQ=255)
                                                                                        cgcCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGcgcg         >  1:2153406/1‑100 (MQ=255)
                                                                                               tACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGcc  <  2:340196/100‑1 (MQ=255)
                                                                                               tACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGcc  <  1:552236/100‑1 (MQ=255)
                                                                                                |                                                                                                  
CTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCC  >  NC_000913/1706933‑1707127

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