Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,926,782 G→A *221* (TAG→TAA) 
D686D (GAC→GAT) 		 exoX →
ptrB ←		 exodeoxyribonuclease 10, DNA exonuclease X
protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,926,7820GA98.0% 162.1 / NA 51*221* (TAG→TAA) 
D686D (GAC→GAT) 		exoX
ptrB		exodeoxyribonuclease 10, DNA exonuclease X
protease II
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/26);  minor base C (0/1);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCG  >  NC_000913/1926647‑1926902
                                                                                                                                       |                                                                                                                        
tGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATcc                                                                                                                       >  1:409725/1‑139 (MQ=255)
    aaaTACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAg                                                                                                                   <  1:1609049/139‑1 (MQ=255)
         ccGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTa                                                                                                              >  2:154809/1‑139 (MQ=255)
          cGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAg                                                                                                             <  2:959532/139‑1 (MQ=255)
                    cGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGt                                                                                                   >  2:1214696/1‑139 (MQ=255)
                                 tGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAgc                                                                                      <  2:1696621/139‑1 (MQ=255)
                                    cGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGcgcg                                                                                   <  2:1402098/139‑1 (MQ=255)
                                      aaCGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                 <  2:1657077/139‑1 (MQ=255)
                                      aaCGCGATCCGGGCTATCTGCGCTGATTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTGGAACATTGAAACATCATCTGGAACATAATTACACCTCAGGCGTAGCGGGTAATAATCCCTCCGCCTGCGAGAc                                                                                 <  1:1360501/139‑1 (MQ=255)
                                           gATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaa                                                                            <  2:1660901/139‑1 (MQ=255)
                                                            cTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCt                                                           >  2:1118122/1‑139 (MQ=255)
                                                            cTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCt                                                           <  1:276135/139‑1 (MQ=255)
                                                             tGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTa                                                          <  1:1061236/139‑1 (MQ=255)
                                                                ttATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                 >  1:1162366/1‑129 (MQ=255)
                                                                ttATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                 <  2:1162366/129‑1 (MQ=255)
                                                                         aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGc                                                                                                            <  2:1669125/77‑1 (MQ=255)
                                                                         aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGc                                                                                                            >  1:1669125/1‑77 (MQ=255)
                                                                         aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCCACGCCTTCGTACg                                              >  1:1705609/1‑139 (MQ=255)
                                                                             tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                 <  2:458280/116‑1 (MQ=255)
                                                                             tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                 >  1:458280/1‑116 (MQ=255)
                                                                               gACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTaa                                        <  1:1214696/139‑1 (MQ=255)
                                                                               gACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAACATATTCCATCGCTACGCCTTCGt                                                 <  2:334373/130‑1 (MQ=255)
                                                                               gACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAACATATTCCATCGCTACGCCTTCGt                                                 >  1:334373/1‑130 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                 >  2:945091/1‑109 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                 <  1:945091/109‑1 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                 <  2:1443856/109‑1 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                 >  1:1443856/1‑109 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGt                                                 <  2:1218012/125‑1 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGt                                                 >  1:1218012/1‑125 (MQ=255)
                                                                                     aTGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACg                                                        <  2:173420/117‑1 (MQ=255)
                                                                                     aTGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACg                                                        >  1:173420/1‑117 (MQ=255)
                                                                                         gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGc                                                                                                                 >  1:1822270/1‑56 (MQ=255)
                                                                                         gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGc                                                                                                                 <  2:1822270/56‑1 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGc                                                                                                                     >  1:1602461/1‑47 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGc                                                                                                                     <  2:1602461/47‑1 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCa                                                                                                                    <  2:380049/48‑1 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCa                                                                                                                    >  1:380049/1‑48 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCgg                                                                                                          <  2:1793640/58‑1 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCgg                                                                                                          >  1:1793640/1‑58 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTAc                         >  2:79007/1‑139 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTAc                         >  1:499426/1‑139 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTAc                         >  1:175441/1‑139 (MQ=255)
                                                                                                         aaCACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt              >  2:589156/1‑139 (MQ=255)
                                                                                                         aaCACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt              <  1:1118122/139‑1 (MQ=255)
                                                                                                              tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc         <  1:789561/139‑1 (MQ=255)
                                                                                                              tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc         >  1:334056/1‑139 (MQ=255)
                                                                                                              tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCCACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc         <  2:1705609/139‑1 (MQ=255)
                                                                                                                aaaCATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCa       <  2:499426/139‑1 (MQ=255)
                                                                                                                     ttATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCg  >  2:514991/1‑139 (MQ=255)
                                                                                                                      tATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAgcgc                                                                                    >  2:1547478/1‑56 (MQ=255)
                                                                                                                      tATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAgcgc                                                                                    <  1:1547478/56‑1 (MQ=255)
                                                                                                                                       |                                                                                                                        
TGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCG  >  NC_000913/1926647‑1926902

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