Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,525 +G intergenic (+92/‑141) b4555 → / → yigM hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,5241.G99.1% 313.1 / ‑3.5 116intergenic (+91/‑142)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (48/67);  total (48/68)
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.90e-01

Only 100 of 121 total aligned reads displayed.

AAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGT  >  NC_000913/4008451‑4008594
                                                                            |                                                                       
acaTGCTGTCGAGGGCGGGGAAAGGCTGCATTATGTGCAGTGCGCGCCAGCGTCTGAAAGACCTTCATCCCGAGCt                                                                          <  1:2584880/74‑1 (MQ=255)
aaaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt                                                                          <  1:3760596/76‑1 (MQ=255)
  aTGCTGGCGATGGCGGGGAAAGGCTGCATGATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg                                                                       <  2:3448801/76‑2 (MQ=255)
   tGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgga                                                                      >  1:730843/1‑74 (MQ=255)
    gCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                                                     <  1:2076657/76‑4 (MQ=255)
       gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                                                                  >  2:3016000/1‑76 (MQ=255)
         cGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                                                                <  2:3796450/76‑1 (MQ=255)
         cGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                                                                <  2:3835462/76‑1 (MQ=255)
         cGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                                                                <  1:2434551/76‑1 (MQ=255)
           aTGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                                                              >  1:3514535/1‑76 (MQ=255)
              gcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                                                           <  2:538722/73‑1 (MQ=255)
               cggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc                                                          <  2:3054200/74‑1 (MQ=255)
                ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                                                         >  1:594061/2‑76 (MQ=255)
                 gggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGCAGTGATTGGTACTa                                                        >  2:1580058/1‑76 (MQ=255)
                 gggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTa                                                        >  1:3900375/1‑76 (MQ=255)
                  ggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAGCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                                       <  2:4094724/76‑1 (MQ=255)
                   gAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                                     <  2:1146351/76‑1 (MQ=255)
                     aaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc                                                    >  2:2895791/4‑76 (MQ=255)
                     aaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc                                                    >  2:737636/4‑76 (MQ=255)
                      aggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                                                   <  1:4548023/74‑1 (MQ=255)
                       ggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg                                                  >  1:1435341/2‑76 (MQ=255)
                       ggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg                                                  <  1:1211114/75‑1 (MQ=255)
                        gCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                                                 <  1:1967565/76‑1 (MQ=255)
                        gCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                                                 <  1:673755/76‑1 (MQ=255)
                          tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacg                                               >  2:4000117/1‑76 (MQ=255)
                          tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacg                                               >  2:4152501/1‑76 (MQ=255)
                           gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGGAGTGATTGGTACTAATGCcgacga                                              >  2:1422163/1‑76 (MQ=255)
                           gCATTAGGGGCAGGGCGCACCAGCGGCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga                                              <  2:150081/76‑1 (MQ=255)
                            cATTATGGGCAGTGCGCGCCAGCGTCTGAAAGACCTGCATCCCGAGCTGTAAGTGATTTGTACTAATGCCgacgac                                             <  1:3405968/76‑1 (MQ=255)
                            cATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGATGTGGAAGTGATTGGTACTAATGCCgacgac                                             <  1:2559158/76‑1 (MQ=255)
                            cATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                                             <  1:2516127/76‑1 (MQ=255)
                            cATTATGGGCAGTGCGCACCAGCGTCTGAAAGAACTTCATCCCGAGCTGGACGGGATTGGTACTAATGCCgatgac                                             >  2:3241215/1‑76 (MQ=255)
                            cATTAGGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                                             <  2:3646296/76‑1 (MQ=255)
                             aTTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGacgacg                                            <  1:2410777/76‑1 (MQ=255)
                              ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc                                           <  1:3542838/76‑1 (MQ=255)
                              tgaTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc                                           <  1:511002/74‑1 (MQ=255)
                               tATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                                          <  2:4465115/76‑1 (MQ=255)
                               tATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                                          <  1:2892318/76‑1 (MQ=255)
                                aTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGCCCACGCgg                                         >  2:4799969/1‑76 (MQ=255)
                                aTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCgg                                         >  2:3589337/1‑76 (MQ=255)
                                aTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCgg                                         >  2:3584729/1‑76 (MQ=255)
                                 tGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt                                        <  1:4103645/76‑1 (MQ=255)
                                 tGGGCAGGGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt                                        <  2:582257/76‑1 (MQ=255)
                                  gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTg                                       <  2:2545943/76‑1 (MQ=255)
                                  gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTg                                       <  1:2213306/76‑1 (MQ=255)
                                   ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGc                                      <  1:898295/76‑1 (MQ=255)
                                    gCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGcc                                     <  2:3238767/76‑1 (MQ=255)
                                         gcgcaccaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAtt                                >  2:2356529/4‑76 (MQ=255)
                                         gcgcaccaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAtt                                <  2:1732972/73‑1 (MQ=255)
                                          cgcaccaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTa                               <  2:3097448/74‑1 (MQ=255)
                                            caccaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGTTATTAATGCAGACGCACCGGTGTCGtattatt                             >  1:3983709/1‑75 (MQ=255)
                                            caccaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc                             <  2:2848143/76‑1 (MQ=255)
                                             accaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt                            >  2:2134340/1‑76 (MQ=255)
                                             accaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt                            >  2:1001007/1‑76 (MQ=255)
                                              ccaGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTg                           <  1:4864459/76‑1 (MQ=255)
                                               caGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAGGCGGTGCCGCATTACCTGc                          <  2:1822726/76‑1 (MQ=255)
                                               caGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc                          <  1:4335048/76‑1 (MQ=255)
                                                aGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCg                         <  1:2019843/76‑1 (MQ=255)
                                                aGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCg                         <  2:1441815/76‑1 (MQ=255)
                                                aGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCg                         >  2:287322/1‑76 (MQ=255)
                                                aGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCg                         <  1:4075634/76‑1 (MQ=255)
                                                 gCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGt                        <  1:4793629/76‑1 (MQ=255)
                                                  cGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTAACTGCGTa                       >  1:3678532/1‑76 (MQ=255)
                                                  cGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGATCCTAATGCCGCCGACGCGGTGCCGCATTATCTGCGTa                       <  1:4512455/76‑1 (MQ=255)
                                                   gTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaa                      <  2:1533823/76‑1 (MQ=255)
                                                   gTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaa                      >  2:417211/1‑76 (MQ=255)
                                                    tCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaaa                     <  1:3764795/76‑1 (MQ=255)
                                                     cTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAc                    <  2:3852194/76‑1 (MQ=255)
                                                     cTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAc                    >  2:1159331/1‑76 (MQ=255)
                                                     cTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAc                    >  1:3484655/1‑76 (MQ=255)
                                                      tGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTg                           >  2:4023065/1‑68 (MQ=255)
                                                      tGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTg                           <  1:4023065/68‑1 (MQ=255)
                                                      tGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAct                   <  1:598939/76‑1 (MQ=255)
                                                       gAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAActc                  <  2:2153935/76‑1 (MQ=255)
                                                       gAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAActc                  >  2:1227101/1‑76 (MQ=255)
                                                       gAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAActc                  >  1:4430839/1‑76 (MQ=255)
                                                         aaGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCCTTATCTGCGTAAActctt                >  2:31234/1‑75 (MQ=255)
                                                         aaGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTa                <  1:71536/76‑1 (MQ=255)
                                                         aaGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTa                <  1:417211/76‑1 (MQ=255)
                                                          aGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga                                              >  1:1176460/1‑45 (MQ=255)
                                                          aGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga                                              <  2:1176459/45‑1 (MQ=255)
                                                          aGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTACACTCTAt               >  2:792859/1‑76 (MQ=255)
                                                          aGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTAt               <  1:1955800/76‑1 (MQ=255)
                                                          aGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGCGCCGCATTATCTGCGTAAACTCTAt               >  2:9189/1‑76 (MQ=255)
                                                           gACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTAtt              >  1:3798502/1‑76 (MQ=255)
                                                             gcTTCGTCCCGAGGTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTa            <  2:88020/75‑1 (MQ=255)
                                                              cTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACTCGGTGCCGCATTATCTGCGTAAACTCTATTTAt           <  2:3181370/76‑1 (MQ=255)
                                                              cTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTAt           <  2:2019076/76‑1 (MQ=255)
                                                               ttCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATc          <  2:3779939/76‑1 (MQ=255)
                                                                tCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCg         <  2:4375132/76‑1 (MQ=255)
                                                                 cATCCCGAGCTGGAAGTGATTGGGACTAATGCCGACGACGAGGTGCCGCATTATCTTCGTAAACGCTATTTATAGt        >  1:4623292/1‑76 (MQ=255)
                                                                    cccGAGCTGGGAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAAt     <  1:1422164/76‑1 (MQ=255)
                                                                      cGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCg   >  1:3100030/1‑76 (MQ=255)
                                                                      cGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCg   >  1:3096441/1‑76 (MQ=255)
                                                                      cGAGCTGGAAGTGATTGGGACGAATGCCGACGACTCTGTGCCGCATTATCTGCGTAAACTCTATTTATCGTaagcg   >  1:3100405/1‑73 (MQ=255)
                                                                      cGAGCTGGAAGAGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATATGCGTAAACTCTATTTATCGTAATCg   >  1:1163193/1‑76 (MQ=255)
                                                                       gAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGt  >  1:4806687/1‑76 (MQ=255)
                                                                       gAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGt  <  1:307316/76‑1 (MQ=255)
                                                                       gAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGt  >  2:2531386/1‑76 (MQ=255)
                                                                       gAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGt  >  2:3430864/1‑76 (MQ=255)
                                                                            |                                                                       
AAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGT  >  NC_000913/4008451‑4008594

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