Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,027,863 C→A 14.2% A28D (GCT→GAT)  PP_4432 → Xaa‑Pro aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,027,8630CA14.2% 74.8 / 2.7 28A28D (GCT→GAT) PP_4432Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (10/14);  new base A (2/2);  total (12/16)
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.53e-01

AACTGAATCACCAGGCAACACAATAAAACCATACCCCAATCATGAGGGATCGATCATGTCTGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGTCGTTCTACGTCCACCAGTGCGTACTGCTGGCCATGGAAGGTGAGCCTGTTTGGTATGGCCGCAACCAAG  >  NC_002947/5027726‑5028003
                                                                                                                                         |                                                                                                                                            
aacTGAATCACCAGGCAACACAATAAAACCATACCCCAATCATGAGGGATCGATCATGTCTGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGa                                                                                                                                  >  1:385365/1‑150 (MQ=255)
       tCACCAGGCAACACAATAAAACCATACCCCAATCATGAGGGATCGATCATGTCTGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAgctgctg                                                                                                                           <  2:385365/150‑1 (MQ=255)
              gCAACACAATAAAACCATACCCCAATCATGAGGGATCGATCATGTCTGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGTCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTa                                                                                                                    <  1:268068/150‑1 (MQ=255)
                  cacaATAAAACCATACCCCAATCATGAGGGATCGATCATGTCTGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGa                                                                                                                <  2:277943/150‑1 (MQ=255)
                  cacaATAAAACCATACCCCAATCATGAGGGATCGATCATGTCTGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGa                                                                                                                >  1:77889/1‑150 (MQ=255)
                     aaTAAAACCATACCCCAATCATGAGGGATCGATCATGTCCGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGAccc                                                                                                             >  2:145517/1‑150 (MQ=255)
                     aaTAAAACCATACCCCAATCATGAGGGATCGATAATGTCCGAAACAATAGTGCACGTTCCTTTTACACGGAAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGATAGGAACATTGAGCTGCTGCTGGTTACTGAccc                                                                                                             >  1:9912/1‑150 (MQ=255)
                      aTAAAACCATACCCCAATCATCAGGGATCGATCATGTCTGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCt                                                                                                            <  1:501920/150‑1 (MQ=255)
                                                       aTGTCTGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGt                                                                           <  2:475291/150‑1 (MQ=255)
                                                          tCTGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTgg                                                                        >  2:63027/1‑150 (MQ=255)
                                                           cTGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGt                                                                       <  1:63027/150‑1 (MQ=255)
                                                              aaaTTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGTCGt                                                                    <  1:153278/150‑1 (MQ=255)
                                                                aTTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGTCGTTc                                                                  <  2:269391/150‑1 (MQ=255)
                                                                       tgAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGTCGTTCTACGTcc                                                           >  1:494166/1‑150 (MQ=255)
                                                                             ttCAGTTTACACGGGAGTAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGATAGGAACATTGAGCTGTTTCTGGTTACTGGCCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGTCGTTCTACGTCCACCAGt                                                     <  2:9912/150‑1 (MQ=255)
                                                                                       acGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTgg                                                                        >  1:368327/1‑121 (MQ=255)
                                                                                       acGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTgg                                                                        <  2:368327/121‑1 (MQ=255)
                                                                                            aggagTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGATAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTg                                                                         >  1:507767/1‑115 (MQ=255)
                                                                                            aggagTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGATAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTg                                                                         <  2:507767/115‑1 (MQ=255)
                                                                                                      ccAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTa                                                                                  >  2:165800/1‑96 (MQ=255)
                                                                                                      ccAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTa                                                                                  <  1:165800/96‑1 (MQ=255)
                                                                                                        aGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGTCGTTCTACGTCCACCAGTGCGTACTGCTGGCCATGGAAGGTGAGc                          <  2:77889/150‑1 (MQ=255)
                                                                                                           gTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGAccc                                                                                                             <  1:383207/64‑1 (MQ=255)
                                                                                                           gTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGAccc                                                                                                             >  2:383207/1‑64 (MQ=255)
                                                                                                                       tCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGTCGTTCTACGTCCACCAGTGCGTACTGCTGGCCATGGAAGGTGAGc                          >  2:280172/1‑135 (MQ=255)
                                                                                                                       tCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGTCGTTCTACGTCCACCAGTGCGTACTGCTGGCCATGGAAGGTGAGc                          <  1:280172/135‑1 (MQ=255)
                                                                                                                        ccGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGTCGTTCTACGTCCACCAGTGCGTACTGCTGGCCATGGAAGGTGAGCCTGTTTGGTATGGCCg          >  1:141683/1‑150 (MQ=255)
                                                                                                                                aaaTGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGTCGTTCTACGTCCACCAGTGCGTACTGCTGGCCATGGAAGGTGAGCCTGTTTGGTATGGCCGCAACCAAg  <  2:494166/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
AACTGAATCACCAGGCAACACAATAAAACCATACCCCAATCATGAGGGATCGATCATGTCTGAAATTGTTGTGAACCTTCCGTTTACACGGGAGGAGTACGCCCAGCGTATCGCCAAAGTCCGTGCAGAAATGCAGGCTAGGAACATTGAGCTGCTGCTGGTTACTGACCCTTCGAACATGGCTTGGCTGACCGGCTACGACGGTTGGTCGTTCTACGTCCACCAGTGCGTACTGCTGGCCATGGAAGGTGAGCCTGTTTGGTATGGCCGCAACCAAG  >  NC_002947/5027726‑5028003

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