Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,498,482 T→C 7.0% I508V (ATC→GTC)  PP_3988 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,498,4820TC 7.0% 166.0 / 3.0 57I508V (ATC→GTC) PP_3988hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base C (2/2);  total (30/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 = 8.78e-01

CTCTTGCTAACAGCACCAGACAAATATTTATGCTTATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTGGGGACAATTTTCTTC  >  NC_002947/4498334‑4498624
                                                                                                                                                    |                                                                                                                                              
ctctTTCTAACAGCAACAGACAAATATTTATGCTTATTCATCAACTGCAGCAATTCACTGTGCAGCCAGTTTTATATGATCTGATAATGCGGATGCGAAACATTAAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGAt                                                                                                                                                <  5:278691/149‑1 (MQ=255)
   ttGCTAACAGCACCAGACAAATATTTATGCTTATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATcaa                                                                                                                                             >  5:240730/1‑148 (MQ=255)
   ttGCTAACAGCACCAGACAAATATTTATGCTTATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATca                                                                                                                                              >  5:123603/1‑148 (MQ=255)
   gtGCTAACAGCACCAGACAAATATTTATGCTTATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAAAGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATc                                                                                                                                               <  6:240730/146‑1 (MQ=255)
             cACCAGACAAATATTTATGCTTATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATctct                                                                                                                                   <  2:185079/149‑1 (MQ=255)
             cACCAGACAAATATTTATGCTTATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATctct                                                                                                                                   <  7:178554/149‑1 (MQ=255)
               ccAGACAAATATTTATGCTTATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTg                                                                                                                                  >  5:123238/1‑148 (MQ=255)
                     aaaTATTTATGCTTATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGAGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGAATCTAAACCTTTGATCACAAAAATCTCTGCAGaaa                                                                                                                            >  5:342922/1‑148 (MQ=255)
                      aaTATTTATGCTTATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATc                                                                                                                          >  2:125215/1‑149 (MQ=255)
                          tttATGCTTATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGcc                                                                                                                      <  6:123238/149‑1 (MQ=255)
                                 ttATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAAcc                                                                                                                >  2:84639/1‑148 (MQ=255)
                                   aTTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAg                                                                                                             >  7:122664/1‑149 (MQ=255)
                                     tcatcaACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGt                                                                                                            >  7:126972/1‑148 (MQ=255)
                                      catcaACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCGCTATCAATATTTAGCGCAGCATCTACACCTTTGATGACAAAATGCTCTGCAGAATTCTGCGGCACCCTAGTac                                                                                                          >  8:287779/1‑148 (MQ=255)
                                        tcaACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGc                                                                                                        >  8:324902/1‑149 (MQ=255)
                                            cTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCt                                                                                                       >  6:263576/1‑146 (MQ=255)
                                                 gCAATGCGCTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAAcc                                                                                               <  7:317011/149‑1 (MQ=255)
                                                  cAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCt                                                                                              <  7:44182/149‑1 (MQ=255)
                                                        aCTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGACCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACt                                                                                        >  7:90995/1‑149 (MQ=255)
                                                        aCTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGACCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACt                                                                                        >  7:96218/1‑149 (MQ=255)
                                                                   aGTTTTTTATGATCTGATAATGCGGATGCGAAACATAGAAAGAATCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAActc                                                                             <  7:202381/149‑1 (MQ=255)
                                                                     ttttttATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTg                                                                                                                                  <  3:39829/94‑1 (MQ=255)
                                                                     ttttttATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTg                                                                                                                                  >  4:39829/1‑94 (MQ=255)
                                                                                 tGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGt                                                               >  5:229126/1‑149 (MQ=255)
                                                                                   aTAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGACCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGtt                                                              <  8:96218/148‑1 (MQ=255)
                                                                                   aTAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGACCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGtt                                                              <  8:90995/148‑1 (MQ=255)
                                                                                     aaTGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTg                                                                                                                                  >  1:331660/1‑78 (MQ=255)
                                                                                     aaTGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTg                                                                                                                                  <  2:331660/78‑1 (MQ=255)
                                                                                       tGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCt                                                         <  1:125215/149‑1 (MQ=255)
                                                                                         cGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTg                                                        <  6:123603/148‑1 (MQ=255)
                                                                                             tGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGgcgc                                                    >  1:123488/1‑148 (MQ=255)
                                                                                                gAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCAtt                                                <  7:287779/149‑1 (MQ=255)
                                                                                                 aaaCATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCAGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCAtt                                                >  5:275039/1‑148 (MQ=255)
                                                                                                     aTTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTc                                           <  5:263576/149‑1 (MQ=255)
                                                                                                     aTTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTc                                           >  3:189859/1‑149 (MQ=255)
                                                                                                     aTTGAAAGATTCCATATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTc                                           <  1:84639/149‑1 (MQ=255)
                                                                                                       tGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCACTGCAGAAATCTGCCGAAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCt                                          >  1:99564/1‑148 (MQ=255)
                                                                                                         aaaGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGa                                        <  7:324902/148‑1 (MQ=255)
                                                                                                         aaaGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGa                                        >  7:202691/1‑148 (MQ=255)
                                                                                                                 ccTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACg                                >  3:291686/1‑148 (MQ=255)
                                                                                                                  cTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGcc                              <  8:122664/149‑1 (MQ=255)
                                                                                                                  cTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGc                               >  3:159586/1‑148 (MQ=255)
                                                                                                                    aTCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCAt                            >  3:95266/1‑149 (MQ=255)
                                                                                                                     tCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCAtt                           <  2:123488/149‑1 (MQ=255)
                                                                                                                       aaTATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAACCAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATtgt                         >  6:314576/1‑149 (MQ=255)
                                                                                                                       aaTATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATtg                          <  4:291686/148‑1 (MQ=255)
                                                                                                                            ttAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCGCAGACTCGTGCTACTTCAAAAAACTCTCATCAACCAAGGTTCCGCTGGCGCCATTGATTGGGATAGGCACGCCATTGTGCTcc                    >  6:73287/1‑149 (MQ=255)
                                                                                                                               gcgcAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTg                 <  4:159586/149‑1 (MQ=255)
                                                                                                                               gcgcAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTg                 <  4:189859/149‑1 (MQ=255)
                                                                                                                                cgcAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTgg                <  5:237433/149‑1 (MQ=255)
                                                                                                                                 gcagcaTCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTggg               >  5:37517/1‑149 (MQ=255)
                                                                                                                                 gcagcaTCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTggg               <  8:202691/149‑1 (MQ=255)
                                                                                                                                       tCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTGGGGACAAt         >  8:350281/1‑149 (MQ=255)
                                                                                                                                           aaCCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTGGGGACAAtt        <  6:275039/146‑1 (MQ=255)
                                                                                                                                            ctcTTTGATCACAAAACTCTCTGCAGAACTCTGCCGCAACCTAGAAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAACCAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTGGGGACAATTttct    <  2:298680/147‑1 (MQ=255)
                                                                                                                                               tttGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTGGGGACAATTttct    >  1:298680/1‑146 (MQ=255)
                                                                                                                                                   aTCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTGGGGACAATTttcttc  <  7:590/144‑1 (MQ=255)
                                                                                                                                                   aTCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTGGGGACAATTttcttc  >  8:590/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CTCTTGCTAACAGCACCAGACAAATATTTATGCTTATTCATCAACTGCCGCAATGCACTGTGCAGCCAGTTTTTTATGATCTGATAATGCGGATGCGAAACATTGAAAGATTCCCTATCAATATTTAGCGCAGCATCTAAACCTTTGATCACAAAAATCTCTGCAGAAATCTGCCGCAACCTAGTAAGCTCAGAAACCTGATACTTCAAAAAACTCTCATCAAACAACGTTCCGCTGGCGCCATTGATTCTGATAAGCACGCCATTGTGCTCCTTGGGGACAATTTTCTTC  >  NC_002947/4498334‑4498624

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