Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,109,062 A→T 7.7% I104N (ATT→AAT)  PP_1881 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,109,0620AT 7.7% 150.3 / 3.4 52I104N (ATT→AAT) PP_1881hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base T (2/2);  total (26/26)
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.87e-01

GTGAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGCGTGC  >  NC_002947/2108916‑2109209
                                                                                                                                                  |                                                                                                                                                   
gtgAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATgg                                                                                                                                                  <  1:1220218/150‑1 (MQ=255)
  gAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGcc                                                                                                                                                <  1:490659/150‑1 (MQ=255)
      cATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCATTGGCCCTgt                                                                                                                                            <  1:636225/150‑2 (MQ=255)
        tatTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGAGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAg                                                                                                                                          >  2:369489/1‑150 (MQ=255)
                gCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTc                                                                                                                                  >  2:256875/1‑150 (MQ=255)
                             aCCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAg                                                                                                                     <  2:152856/150‑1 (MQ=255)
                             aCCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAg                                                                                                                     >  1:248328/1‑150 (MQ=255)
                                 gTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGa                                                                                                                 >  1:861691/1‑150 (MQ=255)
                                        gCGAATCCACAGTTTGCCATTGCGTGGCATGGACGGCTGGTGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAg                                                                                                          <  1:1247833/150‑1 (MQ=255)
                                             tCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTa                                                                                                     >  1:658616/1‑150 (MQ=255)
                                                acaGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCATTGGCCCTgt                                                                                                                                            <  1:624746/108‑2 (MQ=255)
                                                acaGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCATTGGCCCTgt                                                                                                                                            >  2:624746/1‑107 (MQ=255)
                                                       gCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGagta                                                                                           <  2:658616/150‑1 (MQ=255)
                                                                 ggAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGa                                                                                 <  1:256875/150‑1 (MQ=255)
                                                                  gAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGAt                                                                                <  2:1068402/150‑1 (MQ=255)
                                                                     tGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATctg                                                                             <  1:369489/150‑3 (MQ=255)
                                                                             tGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGt                                                                     >  1:599939/1‑150 (MQ=255)
                                                                                 gAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCa                                                                 >  1:26268/1‑150 (MQ=255)
                                                                                 gAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCACTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCa                                                                 <  2:567331/150‑1 (MQ=255)
                                                                                       ccAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCgg                                                                                                                                >  1:237204/1‑81 (MQ=255)
                                                                                       ccAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCgg                                                                                                                                <  2:237204/81‑1 (MQ=255)
                                                                                         aGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGAc                                                         >  1:508653/1‑150 (MQ=255)
                                                                                            gCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCAt                                                      <  2:26268/150‑1 (MQ=255)
                                                                                             cTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATg                                                     >  1:899333/1‑150 (MQ=255)
                                                                                              tAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGc                                                    >  1:1233587/1‑150 (MQ=255)
                                                                                               aaCGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAg                                                                >  1:1153986/1‑137 (MQ=255)
                                                                                               aaCGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAg                                                                <  2:1153986/137‑1 (MQ=255)
                                                                                                   gAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAg                                                                <  1:1151723/133‑1 (MQ=255)
                                                                                                   gAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAg                                                                >  2:1151723/1‑133 (MQ=255)
                                                                                                   gAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCggg                                                 >  1:798116/1‑148 (MQ=255)
                                                                                                   gAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCggg                                                 <  2:798116/148‑1 (MQ=255)
                                                                                                   gAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGa                                               >  1:71961/1‑150 (MQ=255)
                                                                                                   gAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCAACAGGGCCAATGGCCGTGAAGGTGTTTTCGGCGACCGTTAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGa                                               <  1:726501/150‑1 (MQ=255)
                                                                                                    aaaGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGAt                                              <  2:861691/150‑1 (MQ=255)
                                                                                                          tcCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGAc                                        <  1:106518/150‑1 (MQ=255)
                                                                                                                cACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGCCAGGTCAGACCATGCGGGGATCTCGACGGTCAg                                  >  1:232235/1‑150 (MQ=255)
                                                                                                                 aCGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAgg                                 <  2:71961/150‑1 (MQ=255)
                                                                                                                   gAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTc                               >  1:1069524/1‑150 (MQ=255)
                                                                                                                          gggAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTc                        <  2:899333/150‑1 (MQ=255)
                                                                                                                                 aCTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCggg                 >  1:115868/1‑149 (MQ=255)
                                                                                                                                   ttttCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGGTCACGACGTTCAGGTCTGTTTTCATGCGGGCt               >  2:391300/1‑150 (MQ=255)
                                                                                                                                      tCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCAGACGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTtg            <  1:391300/150‑1 (MQ=255)
                                                                                                                                       cttcAGGGCCATTGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTtgt           >  2:1230149/4‑150 (MQ=255)
                                                                                                                                             ggCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGt                            <  1:84443/127‑1 (MQ=255)
                                                                                                                                             ggCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGt                            >  2:84443/1‑127 (MQ=255)
                                                                                                                                              gCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGCGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGCGt    <  2:248328/150‑1 (MQ=255)
                                                                                                                                              gCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCATTGGGGCAGTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGCGt    >  2:1138980/1‑150 (MQ=255)
                                                                                                                                              gCCAATGGACCTAAAGGTGTTTTCGGCGCCAGTAAAGGTGAACTGCAGTTGTAACGTCGACGAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGCGt    <  2:132633/150‑1 (MQ=255)
                                                                                                                                                cAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCggg                                                 >  2:307361/1‑103 (MQ=255)
                                                                                                                                                cAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCggg                                                 <  1:307361/103‑1 (MQ=255)
                                                                                                                                                cAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGCGTGc  >  1:697213/1‑150 (MQ=255)
                                                                                                                                                cAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGAACAAGGCGCGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGCGTGc  >  2:358366/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTGAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGCGTGC  >  NC_002947/2108916‑2109209

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