Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,970,186 G→A 7.1% R267C (CGT→TGT)  tap ← methyl‑accepting protein IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,970,1860GA 7.1% 174.6 / 3.1 56R267C (CGT→TGT) tapmethyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base A (2/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAA  >  NC_000913/1970042‑1970326
                                                                                                                                                |                                                                                                                                            
tttGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGca                                                                                                                                            >  1:197861/1‑147 (MQ=255)
          cGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTAcc                                                                                                                                 >  1:287875/1‑148 (MQ=255)
                  ttAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGcac                                                                                                                                           >  4:308215/1‑130 (MQ=255)
                  ttAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGcac                                                                                                                                           <  3:308215/130‑1 (MQ=255)
                        gCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAg                                                                                                                    <  5:240921/147‑1 (MQ=255)
                          tCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCct                                                                                                                 >  1:40857/1‑148 (MQ=255)
                              tACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCctgctg                                                                                                             <  8:175169/148‑1 (MQ=255)
                               aCTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCctgctg                                                                                                             >  7:60929/1‑147 (MQ=255)
                                   gcggcCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCAt                                                                                                          <  2:197861/146‑1 (MQ=255)
                                          tttGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAg                                                                                                 <  2:287875/148‑1 (MQ=255)
                                            tgtgCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGa                                                                                                <  2:40857/147‑1 (MQ=255)
                                                ccAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTgg                                                                                            <  3:301301/147‑1 (MQ=255)
                                                 cAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGc                                                                                           >  4:278148/1‑147 (MQ=255)
                                                    cGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCaaa                                                                                        >  1:227785/1‑147 (MQ=255)
                                                     gATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCaaaaa                                                                                      <  3:278148/148‑1 (MQ=255)
                                                         ccgccTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGc                                                                                  <  8:60929/148‑1 (MQ=255)
                                                          cgccTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTc                                                                                                   >  8:431824/1‑130 (MQ=255)
                                                          cgccTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTc                                                                                                   <  7:431824/130‑1 (MQ=255)
                                                             ctgctgTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGc                                                                                                                        >  2:249911/1‑106 (MQ=255)
                                                             ctgctgTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGc                                                                                                                        <  1:249911/106‑1 (MQ=255)
                                                               gctgTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGAt                                                                            >  3:23908/1‑148 (MQ=255)
                                                               gctgTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACAGTACCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGa                                                                             >  8:125318/1‑147 (MQ=255)
                                                                ctgTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGa                                                                             >  8:47486/1‑146 (MQ=255)
                                                                ctgTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGa                                                                             <  7:47486/146‑1 (MQ=255)
                                                                  gTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGAtctc                                                                         <  4:23908/148‑1 (MQ=255)
                                                                    tCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCAGTTACCGTCCCACGCAAAg                                                                                                                    <  3:16548/103‑1 (MQ=255)
                                                                    tCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAg                                                                                                                    >  4:16548/1‑103 (MQ=255)
                                                                     agttccGACTTGAGAGCTCGTTATTGCATGCGACAATCTCCGCGATACCAATGTGCTTCTCCTGGCTTCCCTTACGCACATCCCTTACCGTCCAACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCAt                                                                       <  1:337785/142‑1 (MQ=255)
                                                                               tGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTa                                                                                                                                   >  5:42240/1‑77 (MQ=255)
                                                                               tGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTa                                                                                                                                   <  6:42240/77‑1 (MQ=255)
                                                                                tagagaTCGTTATTGCCTGCGAAAATCTCCGCGATACCAATGTGCATCTCCTGGGTTCCCTTACGCACATCACTTACCGTCCCACGCAATGCCTGCTGCATGGTCTTCAGACTGGGAAAAATGGCGGTGATCTCATTACGACCat                                                              <  5:31430/144‑1 (MQ=255)
                                                                                gagagaTCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCatat                                                            <  2:227785/147‑1 (MQ=255)
                                                                                      tCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCaa                                                     >  3:394384/1‑148 (MQ=255)
                                                                                              gCCTGCGACAATCTCCGCTATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCCATCGGAcg                                              <  8:404208/147‑1 (MQ=255)
                                                                                                  gCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTAc                                                                                                                                  <  1:212258/59‑1 (MQ=255)
                                                                                                  gCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTAc                                                                                                                                  >  2:212258/1‑59 (MQ=255)
                                                                                                   cGACAATCTCCGCGATACCAATGTGCATCTCCTGGCGTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGcc                                          <  6:283955/146‑1 (MQ=255)
                                                                                                     aCAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCATCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATAAACCGCAATCGGACGCCCCAGAt                                      >  6:260368/1‑148 (MQ=255)
                                                                                                          ctcCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTAccc                                 >  4:395703/1‑148 (MQ=255)
                                                                                                              gcgATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCgcag                             >  4:394923/1‑148 (MQ=255)
                                                                                                                gATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCgcagca                           <  3:394923/148‑1 (MQ=255)
                                                                                                                gATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCgcagca                           <  7:125318/148‑1 (MQ=255)
                                                                                                                   aCCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCtgctgc                                                                                                            <  1:214912/64‑1 (MQ=255)
                                                                                                                   aCCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCtgctgc                                                                                                            >  2:214912/1‑64 (MQ=255)
                                                                                                                            cATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCaaaa               >  6:144804/1‑148 (MQ=255)
                                                                                                                            cATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACAGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCaaaa               >  7:109587/1‑148 (MQ=255)
                                                                                                                            cATCTCCTGGCTTCCCTTACACACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCAGTCaaaa               >  6:262473/1‑148 (MQ=255)
                                                                                                                            cATCTCCTGGCTTCCCTTACACACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTAAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATAACGCAATCGGACGCGCCAGATTACACGAAGAAATGCTGTCaaaa               >  6:268604/1‑148 (MQ=255)
                                                                                                                              tctcCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATg             >  5:305468/1‑148 (MQ=255)
                                                                                                                                 ccTGGCTTCCCTTACGCACATCACTTACCGTCCCCCGCAAAGCCTGCTGCATGGTCGTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACACGCCAGATTACCCGCAGCAATGCTGTCAAAATGGc           >  4:170918/1‑147 (MQ=255)
                                                                                                                                  cTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTa         >  2:379590/1‑148 (MQ=255)
                                                                                                                                         cccTTACACACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATa   <  5:262473/147‑1 (MQ=255)
                                                                                                                                         cccTTACACACATAAATTACCGTCCCACGCAAAGCATGCTGCATTGTCTTCAGACTGGCAAAAATGGAGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATa   <  5:268604/147‑1 (MQ=255)
                                                                                                                                          ccTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATaa  <  4:394384/147‑1 (MQ=255)
                                                                                                                                             tACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATa   >  4:33021/1‑143 (MQ=255)
                                                                                                                                             tACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATa   <  3:33021/143‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
TTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTTCAGACTGGCAAAAATGGCGGTGATCTCATTACGACCATATACCGCAATCGGACGCGCCAGATTACCCGCAGCAATGCTGTCAAAATGGCTACCGATAA  >  NC_000913/1970042‑1970326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑