Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,310,449 G→C 9.4% L364V (CTG→GTG)  PP_1144 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,310,4490GC 9.4% 150.9 / 5.7 51L364V (CTG→GTG) PP_1144membrane protein
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base C (3/2);  total (26/25)
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.40e-01

GGTGGTCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTG  >  NC_002947/1310306‑1310593
                                                                                                                                               |                                                                                                                                                
ggtggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTg                                                                                                                                            <  1:887751/150‑1 (MQ=255)
ggtggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTg                                                                                                                                            <  1:847450/150‑1 (MQ=255)
     tCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGt                                                                                                                                             <  2:946825/144‑1 (MQ=255)
     tCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGt                                                                                                                                             >  1:946825/1‑144 (MQ=255)
       aTGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCgg                                                                                                                                     >  2:231508/1‑150 (MQ=255)
            cATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTTCGGGTAGc                                                                                                                                >  2:364/1‑150 (MQ=255)
            cATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCgg                                                                                                                                     >  2:259143/1‑145 (MQ=255)
            cATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCgg                                                                                                                                     <  1:259143/145‑1 (MQ=255)
            cATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGc                                                                                                                                >  1:1078600/1‑150 (MQ=255)
             aTGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGATGGGCCAGGTGGCGGGCGTCGTTGCGGCTGTTGAGGCTGGTGTCCACAAAGAAGTCGTCGGGACCAAAGCGCGCAACCAGGGCGGGGTTGCGGGTAGCg                                                                                                                               >  2:222002/1‑150 (MQ=255)
               gcgcCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATTCGCTTGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGctgt                                                                                                                             <  1:1342467/150‑4 (MQ=255)
               gcgcCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGc                                                                                                                             <  2:678917/150‑1 (MQ=255)
               gcgcCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGc                                                                                                                             <  1:272140/150‑1 (MQ=255)
                 gccgccCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCACGACCTGGTTGCgcgcttt                                                                                                                                <  2:123305/145‑7 (MQ=255)
                 gccgccCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCACGACCTGGTTGCgcgcttt                                                                                                                                >  1:123305/1‑139 (MQ=255)
                   cgcccgccCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTg                                                                                                                         >  1:633424/1‑150 (MQ=255)
                          ccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGt                                                                                                                  >  2:1003424/1‑150 (MQ=255)
                            aGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCAGGGCCCGGAGGACCGCGTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTg                                                                                                                >  1:602938/1‑150 (MQ=255)
                                    ggggTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCAc                                                                                                        <  1:389329/150‑1 (MQ=255)
                                          tttCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGag                                                                                                  <  1:1394626/150‑1 (MQ=255)
                                              aTCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagaagc                                                                                              >  2:176681/1‑150 (MQ=255)
                                               tCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGGCCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTg                                                                                                                                            >  2:1266243/1‑103 (MQ=255)
                                               tCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCGTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTg                                                                                                                                            <  1:1266243/103‑1 (MQ=255)
                                                ccGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTc                                                                                            <  1:1003424/150‑1 (MQ=255)
                                                                tGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTgg                                                                            <  2:633424/150‑1 (MQ=255)
                                                                      aGGTGGCGGCCTTCGTCGCGGCTGTTCACGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAgg                                                                      <  2:678373/150‑1 (MQ=255)
                                                                                 ttCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGAt                                                           >  2:1164945/1‑150 (MQ=255)
                                                                                             gTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGtt                                                              >  2:983203/1‑135 (MQ=255)
                                                                                             gTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGtt                                                              <  1:983203/135‑1 (MQ=255)
                                                                                              ttCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATgg                                              >  1:270345/1‑150 (MQ=255)
                                                                                                  ggCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGa                                          <  2:1008527/150‑1 (MQ=255)
                                                                                                          aCCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAgg                                  >  1:1418209/1‑150 (MQ=255)
                                                                                                          aCCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAgg                                  <  1:222002/150‑1 (MQ=255)
                                                                                                              cacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGa                                                            >  1:595399/1‑120 (MQ=255)
                                                                                                              cacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGa                                                            <  2:595399/120‑1 (MQ=255)
                                                                                                                 aGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGa                           <  2:868469/150‑1 (MQ=255)
                                                                                                                     cTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGc                       <  1:231508/150‑1 (MQ=255)
                                                                                                                     cTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGc                       >  1:111487/1‑150 (MQ=255)
                                                                                                                        gtcgCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAg                    >  2:1072433/1‑150 (MQ=255)
                                                                                                                        gtcgCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAg                    >  2:687240/1‑150 (MQ=255)
                                                                                                                          cgCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGc                  <  2:270345/150‑1 (MQ=255)
                                                                                                                              accaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGtt              >  1:355497/1‑150 (MQ=255)
                                                                                                                              accaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGtt              <  1:687240/150‑1 (MQ=255)
                                                                                                                              accaAAGCGCGCAACCAGGTCGTGGATGTTGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGtt              >  1:750031/1‑150 (MQ=255)
                                                                                                                                 aaaGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTc           >  2:210484/1‑150 (MQ=255)
                                                                                                                                 aaaGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGATGGATATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTc           <  2:750031/150‑1 (MQ=255)
                                                                                                                                   tacccgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGc         >  2:158182/5‑150 (MQ=255)
                                                                                                                                   aGCGCGCAACCAGGTCGTGGTTGAGGAAAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGc         <  1:364/150‑1 (MQ=255)
                                                                                                                                      ccgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGaaa                                                 >  2:249042/2‑107 (MQ=255)
                                                                                                                                      ccgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGaaa                                                 <  1:249042/106‑1 (MQ=255)
                                                                                                                                          aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTg  >  1:821028/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GGTGGTCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTG  >  NC_002947/1310306‑1310593

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