Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,734,055 G→A 9.1% V18I (GTC→ATC)  PP_3301 → RND transporter membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,734,0550GA 9.1% 121.2 / 3.5 44V18I (GTC→ATC) PP_3301RND transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base G (23/17);  new base A (2/2);  total (25/19)
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.95e-01

AGAAGTCTCGCCCCTTAGCAGGGTGTGACCATTGACAATATTGGTCACTGGTAGAAAAATGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGTCGCTCGCACACAGAGTGACCTAGGCTTCCGA  >  NC_002947/3733909‑3734201
                                                                                                                                                  |                                                                                                                                                  
aGAAGTCTCGCCCCTTAGCAGGGTGTGACCATTGACAATATTGGTCACTGGTAGAAAAATGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTATc                                                                                                                                                  >  8:152035/1‑149 (MQ=255)
aGAAGTATCGCCCCTTAGCAGGGTGTGACCATTGACAATATTGGTCACTGGTAGAAAAATGGGATTCCTATATCGCCTTTTGAGAAAGGGATCACATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGACCTTATc                                                                                                                                                  >  8:156098/1‑149 (MQ=255)
  aaGTCTCGCCCCTTAGCAGGGTGTGACCATTGACAATATTGGTCACTGGTAGAAAAATGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGc                                                                                                                                                >  5:220226/1‑149 (MQ=255)
     tctcGCCCCTTAGCAGGGTGTGACCATTGACAATATTGGTCACTGGTAGAAAAATGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTg                                                                                                                                             >  5:3987/1‑149 (MQ=255)
       tcGCCCCTTAGCAGGGTGTGACCATTGACAATATTGGTCACTGGTAGAAAAATGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCg                                                                                                                                           >  3:87064/1‑149 (MQ=255)
           cccTTAGCAGGGTTTGACCATTGACAATATTGGTCACTGGTAGAAAAATGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGa                                                                                                                                       <  1:131150/149‑1 (MQ=255)
                     ggTGTGACCATTGACAATATTGGTCACTGGTAGAAAAATGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCt                                                                                                                             <  6:220226/149‑1 (MQ=255)
                               ttGACAATATTGGTCACTGGTAGAAAAATGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGAcc                                                                                                                   <  2:215943/149‑1 (MQ=255)
                                  aCAATATTGGTCACTGGTAGAAAAATGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCga                                                                                                                <  6:3987/149‑1 (MQ=255)
                                                        aaaTGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGCCCCGAGAGACACACCGCCGGTCGGGa                                                                                           >  5:244722/1‑148 (MQ=255)
                                                           tGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTc                                                                                                                                                  <  6:25010/90‑1 (MQ=255)
                                                           tGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTc                                                                                                                                                  >  5:25010/1‑90 (MQ=255)
                                                                   cTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGacacac                                                                                                       >  3:204845/1‑125 (MQ=255)
                                                                   cTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGacacac                                                                                                       <  4:204845/125‑1 (MQ=255)
                                                                   cTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCg                                                                                <  4:87064/148‑1 (MQ=255)
                                                                       atCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTa                                                                             >  3:204504/1‑147 (MQ=255)
                                                                        tCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAg                                                                          >  6:63037/1‑149 (MQ=255)
                                                                         cGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGt                                                                                             <  7:180733/129‑1 (MQ=255)
                                                                         cGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGt                                                                                             >  8:180733/1‑129 (MQ=255)
                                                                                   gaAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGc                                                               >  1:178342/1‑149 (MQ=255)
                                                                                             ccATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTAGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCAcgt                                                     >  4:149202/1‑149 (MQ=255)
                                                                                                       gCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCg                                           >  7:25600/1‑149 (MQ=255)
                                                                                                           ccTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGag                                        >  7:127778/1‑148 (MQ=255)
                                                                                                            cTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTATCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGCGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGagta                                      <  7:156098/149‑1 (MQ=255)
                                                                                                            cTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTATCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGagta                                      <  7:152035/149‑1 (MQ=255)
                                                                                                                cgccATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTgg                                  >  6:135886/1‑149 (MQ=255)
                                                                                                                 gccATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAg                                                                          <  5:89817/108‑1 (MQ=255)
                                                                                                                 gccATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAg                                                                          >  6:89817/1‑108 (MQ=255)
                                                                                                                   cATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGtcg                               >  4:195330/1‑149 (MQ=255)
                                                                                                                     tAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGtcgct                             <  2:146398/149‑1 (MQ=255)
                                                                                                                       gTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGtcgctc                            <  8:127778/148‑1 (MQ=255)
                                                                                                                        tCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACAGCAGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGtcgctcg                           >  5:89769/1‑148 (MQ=255)
                                                                                                                          ggTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGTCGCTCGca                         <  5:135886/148‑1 (MQ=255)
                                                                                                                                 tACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGTCGCTCGCACACAGAGt                 >  1:154255/1‑149 (MQ=255)
                                                                                                                                    tCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGagta                                      >  6:11213/1‑125 (MQ=255)
                                                                                                                                    tCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGagta                                      <  5:11213/125‑1 (MQ=255)
                                                                                                                                    tCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGtcgctc                            >  6:186338/1‑135 (MQ=255)
                                                                                                                                    tCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGtcgctc                            <  5:186338/135‑1 (MQ=255)
                                                                                                                                      ccACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGTCGCTCGCACACAGAGTGAcc             >  5:127807/1‑148 (MQ=255)
                                                                                                                                        aCTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGTCGCTCGCACACAGAGTGACCTAg          <  6:89769/149‑1 (MQ=255)
                                                                                                                                             cccTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCg                                           >  2:53305/1‑111 (MQ=255)
                                                                                                                                             cccTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCg                                           <  1:53305/111‑1 (MQ=255)
                                                                                                                                              ccTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTgg                                  <  3:71557/119‑1 (MQ=255)
                                                                                                                                              ccTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTgg                                  >  4:71557/1‑119 (MQ=255)
                                                                                                                                                  gTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGTCGCTCGCACACAGAGTGACCTAGGCTTCCGa  <  2:178342/147‑1 (MQ=255)
                                                                                                                                                  gTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGTCGCTCGCACACAGAGTGACCTAGGCTTCCGa  <  4:204504/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGAAGTCTCGCCCCTTAGCAGGGTGTGACCATTGACAATATTGGTCACTGGTAGAAAAATGGGATTCCTATATCGCCTTTTGAGAAAGGGATCCCATGCGCCTGCGTCCTGCCGCCATAGTCGGTTGCTTACTCCCACTCGCCCTTGTCGCTTGCGGGGACTCTTCACCTGGAAAAGACCCGAGAGACACACCGCCGGTCGTGAAAGCGGTCTCGGTACAGGACGCGTCAGCCCTTGCACGTTCGTTTTCCGGAGTAGTGGTCGCTCGCACACAGAGTGACCTAGGCTTCCGA  >  NC_002947/3733909‑3734201

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: ‑