Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,260,338 T→A 9.2% Q161L (CAG→CTG)  PP_2857 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,260,3380TA 9.2% 148.1 / 4.3 54Q161L (CAG→CTG) PP_2857hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base A (2/3);  total (26/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 = 8.96e-01

TGGTGATTCTCCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGAT  >  NC_002947/3260191‑3260465
                                                                                                                                                   |                                                                                                                                
tGGTGATTCTCCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTg                                                                                                                                 >  6:191545/1‑149 (MQ=255)
tGGTGAATCTCCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTg                                                                                                                                 <  2:48928/149‑1 (MQ=255)
  gTGATTCTCCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTgg                                                                                                                                <  8:260465/148‑1 (MQ=255)
       tctcCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACta                                                                                                                          >  2:104979/1‑148 (MQ=255)
           cGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCAc                                                                                                                            >  5:299113/1‑143 (MQ=255)
           cGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCAc                                                                                                                            <  6:299113/143‑1 (MQ=255)
            gTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGt                                                                                                                     <  5:191545/149‑1 (MQ=255)
               cGCTATTAGGTTAAGTGTCGTTCATAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGa                                                                                                                   <  8:277164/148‑1 (MQ=255)
                    ttAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGACCAgg                                                                                                              <  3:146076/148‑1 (MQ=255)
                    ttAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGACCAgg                                                                                                              <  1:104979/148‑1 (MQ=255)
                       ggTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGAt                                                                                                          <  1:166715/149‑1 (MQ=255)
                        gTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGa                                                                                                           >  7:12506/1‑147 (MQ=255)
                        gTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGa                                                                                                           >  2:251513/1‑147 (MQ=255)
                        gTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATa                                                                                                         <  7:217677/149‑1 (MQ=255)
                        gTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATAACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGa                                                                                                           >  2:247718/1‑147 (MQ=255)
                         ttAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCAGTCGGTTTCCCCCCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGGCACTGGTAGTGACCTGGCGATAg                                                                                                        >  7:61763/1‑149 (MQ=255)
                            aGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCct                                                                                                     >  3:308585/1‑149 (MQ=255)
                             gtgtCGTTCTTAACTGCAGGGTCAATCCGCGACCCAGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCctg                                                                                                    >  1:207517/1‑149 (MQ=255)
                                tCGTTCTTAACTGCAGGGTCAATCCGCGACCCAGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCctgctg                                                                                                 <  2:207517/149‑1 (MQ=255)
                                  gTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTgg                                                                                                >  8:79983/1‑148 (MQ=255)
                                   ttCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGaa                                                                                              >  2:200388/1‑149 (MQ=255)
                                   ttCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGaa                                                                                              >  6:97304/1‑149 (MQ=255)
                                   ttCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACACCGGCACTCTCGTTCCACTGGATCACACCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGAGATAGCCTGCTGGaa                                                                                              >  6:229213/1‑149 (MQ=255)
                                      ttAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGt                                                                                           >  3:137899/1‑149 (MQ=255)
                                        aaCTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGt                                                                                           <  4:308585/147‑1 (MQ=255)
                                        aaCTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGt                                                                                           <  7:79983/147‑1 (MQ=255)
                                                ggTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGAc                                                                                                                  <  6:294683/116‑1 (MQ=255)
                                                ggTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGAc                                                                                                                  >  5:294683/1‑116 (MQ=255)
                                                      tCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgtt                                                                           >  2:233370/1‑149 (MQ=255)
                                                         gcgACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttgtt                                                                        <  8:12506/149‑1 (MQ=255)
                                                              ccGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGAt                                                                     >  6:201520/1‑147 (MQ=255)
                                                                          gtggtTGACCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCt                                                       <  8:79812/146‑1 (MQ=255)
                                                                           gtgtTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTg                                                      >  8:182898/1‑149 (MQ=255)
                                                                            tgtTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGt                                                     <  5:201520/149‑1 (MQ=255)
                                                                                ggCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑Gc                                                 <  1:200388/148‑1 (MQ=255)
                                                                                ggCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCa                                                >  4:185424/1‑149 (MQ=255)
                                                                                 gCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttg                                                                          <  8:86280/123‑1 (MQ=255)
                                                                                 gCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttg                                                                          >  7:86280/1‑123 (MQ=255)
                                                                                 gCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCAt                                               <  3:309289/149‑1 (MQ=255)
                                                                                        ggTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGAc                                         <  7:182898/148‑1 (MQ=255)
                                                                                         gTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACg                                        >  4:121857/1‑148 (MQ=255)
                                                                                                         ctctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTg                       <  4:137899/149‑1 (MQ=255)
                                                                                                          tctGGTTCAACTGGATCACCCCTCAGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTGc                      >  7:196727/1‑149 (MQ=255)
                                                                                                            tGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTGc                      <  1:247718/147‑1 (MQ=255)
                                                                                                            tGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTGc                      <  1:251513/147‑1 (MQ=255)
                                                                                                                   aCTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGACAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCa             >  5:104957/1‑149 (MQ=255)
                                                                                                                         tCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATg       >  7:247610/1‑149 (MQ=255)
                                                                                                                             cccTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGa   <  3:121857/149‑1 (MQ=255)
                                                                                                                              ccTGGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGAGCTGGGTATAGCCTGCTGTAATGTCAGTTGCTCTTGTGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTGGGGTTTTGCACCAATGCTGAt  <  8:196727/149‑1 (MQ=255)
                                                                                                                              ccTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTCGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGa   <  1:233370/149‑1 (MQ=255)
                                                                                                                                          gTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAg                                                                  >  3:32437/1‑74 (MQ=255)
                                                                                                                                          gTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAg                                                                  <  4:32437/74‑1 (MQ=255)
                                                                                                                                                gCCAGGTCACTACCAGTGACCAGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTg    >  6:274443/1‑129 (MQ=255)
                                                                                                                                                gCCAGGTCACTACCAGTGACCAGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTg    <  5:274443/129‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                
TGGTGATTCTCCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTT‑GCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGAT  >  NC_002947/3260191‑3260465

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