Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,036,373 C→T *578* (TAG→TAA recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,036,3730CT98.2% 178.9 / NA 57*578* (TAG→TAArecJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/29);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGCTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTAAAAG  >  NC_000913/3036244‑3036506
                                                                                                                                 |                                                                                                                                     
gTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGc                                                                                                                              <  2:1527139/139‑1 (MQ=255)
   cTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAg                                                                                                                           >  2:402346/1‑139 (MQ=255)
    tGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGcc                                                                                                                             >  2:1228533/1‑136 (MQ=255)
    tGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGcc                                                                                                                             <  1:1228533/136‑1 (MQ=255)
      aaTGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAGATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGata                                                                                                                        <  1:1312233/139‑1 (MQ=255)
      aaTGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGata                                                                                                                        <  1:136868/139‑1 (MQ=255)
          cGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGt                                                                                                                    >  2:1353839/1‑139 (MQ=255)
                                        gTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAAt                                                                                                                                  >  1:614874/1‑95 (MQ=255)
                                        gTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAAt                                                                                                                                  <  2:614874/95‑1 (MQ=255)
                                          cTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGa                                                                                    <  2:1589300/139‑1 (MQ=255)
                                          cTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGa                                                                                    <  1:895642/139‑1 (MQ=255)
                                                tCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCg                                                                              <  2:470545/139‑1 (MQ=255)
                                                    ttATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGtt                                                                            <  1:751231/137‑1 (MQ=255)
                                                    ttATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGtt                                                                            >  2:751231/1‑137 (MQ=255)
                                                          aGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                    <  1:1353839/139‑1 (MQ=255)
                                                          aGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                    <  2:570411/139‑1 (MQ=255)
                                                            cTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGa                                                                                    <  2:682524/121‑1 (MQ=255)
                                                            cTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGa                                                                                    >  1:682524/1‑121 (MQ=255)
                                                               gTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTta                                                               >  2:1494261/1‑139 (MQ=255)
                                                               gTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTta                                                               >  1:1270069/1‑139 (MQ=255)
                                                                tCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                                  >  1:1168315/1‑119 (MQ=255)
                                                                tCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGaaa                                                                                  <  2:1168315/119‑1 (MQ=255)
                                                                    aaTGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGc                                                          <  2:1270069/139‑1 (MQ=255)
                                                                    aaTGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGc                                                          <  1:402346/139‑1 (MQ=255)
                                                                          gTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGTTAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTg                                                    >  1:374491/1‑139 (MQ=255)
                                                                                    cGGGATCGTACCCAATCCACGCTCTTTTTCATAGCGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCCGTTGCCGCGAAACTCGTTGAGATCGAGCTTATAAGCCAGTTGGACTTCgcgc                                          <  2:1086901/139‑1 (MQ=255)
                                                                                     gggATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCg                                                                     <  1:769212/111‑1 (MQ=255)
                                                                                     gggATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCg                                                                     >  2:769212/1‑111 (MQ=255)
                                                                                                   tcgaCGCTCTTTTTTATAGAGCAGATGACGCTCAATTGGCCAGATATTGACGATGATAATTTGCTGGCTGCGGTTACCGCGAAACTCGTTGATATGGAGCGCATAAGCCAGTTGCAATTCGCGGTCGCCGTTATCCGGc                           >  1:1133043/4‑139 (MQ=255)
                                                                                                   tCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGc                           >  1:1142708/1‑139 (MQ=255)
                                                                                                   tCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGc                           >  2:749007/1‑139 (MQ=255)
                                                                                                   tCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGc                           >  1:508442/1‑139 (MQ=255)
                                                                                                        gCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCagag                      >  1:535918/1‑139 (MQ=255)
                                                                                                          tctTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAggg                    <  1:1097732/139‑1 (MQ=255)
                                                                                                             tttttATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGcc                                     <  1:706089/119‑1 (MQ=255)
                                                                                                             tttttATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGcc                                     >  2:706089/1‑119 (MQ=255)
                                                                                                                  ataGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCg            <  1:749007/139‑1 (MQ=255)
                                                                                                                    agagAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGAc          >  2:1388223/1‑139 (MQ=255)
                                                                                                                      agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                    >  2:1385198/1‑79 (MQ=255)
                                                                                                                      agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                    <  1:1385198/79‑1 (MQ=255)
                                                                                                                      agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTcgcg                                           <  2:864224/104‑1 (MQ=255)
                                                                                                                      agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTcgcg                                           >  1:864224/1‑104 (MQ=255)
                                                                                                                      agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCgcgc                                          >  2:1600859/1‑105 (MQ=255)
                                                                                                                      agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCgcgc                                          <  1:1600859/105‑1 (MQ=255)
                                                                                                                      agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACAt        >  2:1458036/1‑139 (MQ=255)
                                                                                                                       gaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACAtt       <  1:1342455/139‑1 (MQ=255)
                                                                                                                        aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGcc                                     >  2:388433/1‑108 (MQ=255)
                                                                                                                        aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGcc                                     <  1:388433/108‑1 (MQ=255)
                                                                                                                        aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCg            >  2:689029/1‑133 (MQ=255)
                                                                                                                        aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCg            <  1:689029/133‑1 (MQ=255)
                                                                                                                        aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTa      >  1:1462346/1‑139 (MQ=255)
                                                                                                                        aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTa      >  2:715450/1‑139 (MQ=255)
                                                                                                                        aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTa      >  2:1268971/1‑139 (MQ=255)
                                                                                                                          gaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTaa     >  1:1458566/1‑138 (MQ=255)
                                                                                                                          gaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTaa     <  2:1458566/138‑1 (MQ=255)
                                                                                                                            tGACGTTAAATTGGCCAGATATTGTCGATGTTAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTAAAAg  <  2:374491/139‑1 (MQ=255)
                                                                                                                            tGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCAAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTAAAAg  <  1:1388223/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                     
GTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGCTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTAAAAG  >  NC_000913/3036244‑3036506

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