Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,747,389 Δ1 bp coding (572/591 nt) b4540 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,747,3880C.99.2% 537.9 / ‑4.2 119coding (571/591 nt)b4540hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (63/55);  total (64/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

TTGCGCACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATCCTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGA  >  NC_000913/2747313‑2747459
                                                                           |                                                                       
ttGCGCACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATc                                                                         >  1:5141024/1‑76 (MQ=255)
tcgcgcACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATc                                                                         <  1:2793264/74‑1 (MQ=255)
  gcgcACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATctg                                                                       >  2:4345996/1‑74 (MQ=255)
  gcgcACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATctg                                                                       <  2:199290/76‑3 (MQ=255)
  gcgcACTATCGGGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATcta                                                                       <  2:5537172/76‑3 (MQ=255)
     cACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGaa                                                                   <  2:4910599/76‑1 (MQ=255)
      aCTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGaaa                                                                  <  1:1032215/76‑1 (MQ=255)
      aCTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGaaa                                                                  <  2:538472/76‑1 (MQ=255)
      aCTATCGGGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGaaa                                                                  <  2:5805769/76‑1 (MQ=255)
       cTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGaaca                                                                 >  2:2949207/1‑74 (MQ=255)
        caTCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGaaaaa                                                                <  2:5636093/75‑1 (MQ=255)
          tCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAAt                                                              <  2:4425732/76‑1 (MQ=255)
           tgggCACGAATGACGCGCGCAACTGTCCGTTGGCAATCGGGATCTGCTGCTGTCGGTGCTCGCT‑CTGGAAACAATg                                                             <  2:5785855/73‑1 (MQ=255)
            gTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGa                                                            >  1:1500478/1‑76 (MQ=255)
            gTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGa                                                            <  1:38176/76‑1 (MQ=255)
            gTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGa                                                            <  1:4983375/76‑1 (MQ=255)
            gTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGa                                                            <  2:1209231/76‑1 (MQ=255)
              gCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGGT‑CTGGGAAAAATGAcc                                                          >  1:4028981/1‑76 (MQ=255)
              gCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGAcc                                                          <  1:1730421/76‑1 (MQ=255)
              gCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGAcc                                                          >  1:895100/1‑76 (MQ=255)
              gCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGAcc                                                          >  1:1506381/1‑76 (MQ=255)
               cACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCg                                                         <  1:1422221/76‑1 (MQ=255)
                aCGACTCCCGCTCACAAATTTCCCGTCGCAATCAAGATCTGCTGCTGTCGGTGCTGGAT‑CTGGAAAACACGAcctt                                                        <  1:2877227/76‑3 (MQ=255)
                aCGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGt                                                        >  2:4577753/1‑76 (MQ=255)
                aCGAATCACGCTCACAAATTTCCCCTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGt                                                        >  2:4581919/1‑76 (MQ=255)
                 cGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGtt                                                       >  2:4459321/1‑76 (MQ=255)
                 cGAATCACGATCACAAATTTCCCGACGCAATCAGGATATGCTGCTGTCGGCGCTCGGT‑CTGGGAAAAATGACCGtt                                                       >  1:1709349/1‑76 (MQ=255)
                  gAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTtg                                                      >  2:5175808/1‑76 (MQ=255)
                   aaTCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTtga                                                     >  1:5693285/1‑76 (MQ=255)
                    aTCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTtgat                                                    <  1:2520182/76‑1 (MQ=255)
                     gcACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTtgatg                                                   <  2:503675/75‑1 (MQ=255)
                      cACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTtgatga                                                  >  1:3181830/1‑76 (MQ=255)
                       aCGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAc                                                 >  1:89743/1‑76 (MQ=255)
                         gCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAcat                                               >  1:3982869/1‑76 (MQ=255)
                          cTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAcatc                                              >  2:3106775/1‑76 (MQ=255)
                          cTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAcatc                                              >  1:5789756/1‑76 (MQ=255)
                           tCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGATCGTTGATGAcatca                                             >  1:3179809/1‑76 (MQ=255)
                           tCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAcatca                                             >  2:1386174/1‑76 (MQ=255)
                            cacaAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAcatcat                                            >  1:166720/1‑76 (MQ=255)
                             acaAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATg                                           >  2:453310/1‑76 (MQ=255)
                               aaaTTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGt                                         >  2:2824132/1‑76 (MQ=255)
                                aaTTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTg                                        <  1:5558158/76‑1 (MQ=255)
                                 aTTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGgcgc                                       >  1:5019513/1‑76 (MQ=255)
                                  tttCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGcc                                      <  1:2740598/76‑1 (MQ=255)
                                   ttCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCcg                                     <  2:2544133/76‑1 (MQ=255)
                                    tCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGcgccgc                                    >  2:5173980/1‑76 (MQ=255)
                                    tCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCcgc                                    >  1:273980/1‑76 (MQ=255)
                                    tCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCcgc                                    >  2:3306138/1‑76 (MQ=255)
                                     cccGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCcgcg                                   >  1:1179166/1‑76 (MQ=255)
                                     cccGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCcgcg                                   <  2:5021546/76‑1 (MQ=255)
                                      ccGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCgcgc                                  >  2:1475449/1‑76 (MQ=255)
                                      ccGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCgcgc                                  >  2:2522300/1‑76 (MQ=255)
                                       cGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCa                                 <  2:5554216/76‑1 (MQ=255)
                                        gTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAg                                >  1:1012800/1‑76 (MQ=255)
                                        gTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAg                                >  1:1962871/1‑76 (MQ=255)
                                        ctCGCCATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAg                                <  2:113847/75‑1 (MQ=255)
                                         tCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGt                               <  2:4418967/76‑1 (MQ=255)
                                          cGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTg                              >  2:5282977/1‑76 (MQ=255)
                                          cGCAATCAGGATATGCTGCTGTCGGGGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGAAgcg                              >  1:4521008/1‑74 (MQ=255)
                                           gCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGa                             >  1:5508602/1‑76 (MQ=255)
                                           gCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGa                             <  1:5261456/76‑1 (MQ=255)
                                           gCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGa                             >  2:545966/1‑76 (MQ=255)
                                           gCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGa                             <  2:3634695/76‑1 (MQ=255)
                                             aaTCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGaaa                           <  1:2522300/76‑1 (MQ=255)
                                              aTCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATCGTGTCGCGCAGTGAAAt                          >  1:5129587/1‑76 (MQ=255)
                                                cAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTa                        <  1:5198066/76‑1 (MQ=255)
                                                 aGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTAt                       >  1:4322065/1‑76 (MQ=255)
                                                 aGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTAt                       <  1:2923031/76‑1 (MQ=255)
                                                 aGGATATGCTGCTGTCGGAGCTCGAT‑CTGGAAAAAAGGACCGTTGATGACCTCATGGGGCCGCGCAGTGAACttcc                       >  2:1623240/1‑74 (MQ=255)
                                                  ggATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATc                      <  1:3605122/76‑1 (MQ=255)
                                                  ggATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATc                      <  1:357920/76‑1 (MQ=255)
                                                   gATATGCTGCTCTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCg                     <  1:3820326/76‑1 (MQ=255)
                                                    atatGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGATGCCGCGCAGTGAAATTATCgg                    <  1:4570766/76‑1 (MQ=255)
                                                     tatGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGt                   <  1:3555815/76‑1 (MQ=255)
                                                      atGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTa                  <  1:3695848/76‑1 (MQ=255)
                                                       tgctgcTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTAt                 <  2:241350/76‑1 (MQ=255)
                                                       tgctgcTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTAt                 >  2:1719182/1‑76 (MQ=255)
                                                        gctgctGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTAtt                <  1:2918778/76‑1 (MQ=255)
                                                        gctgctGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTAtt                >  1:2032568/1‑76 (MQ=255)
                                                        gctgctGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTAtt                >  2:1587508/1‑76 (MQ=255)
                                                         ctgctgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTAtcg               <  1:3992541/76‑3 (MQ=255)
                                                         ctgctgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTg               >  1:3353605/1‑76 (MQ=255)
                                                          tgctgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGa              >  1:3237235/1‑76 (MQ=255)
                                                           gctgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGat             >  1:3261599/1‑76 (MQ=255)
                                                             tgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGatat           >  1:161961/1‑76 (MQ=255)
                                                              gTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATc          >  1:2384216/1‑76 (MQ=255)
                                                               tCGGTGCTCGAT‑CTGGAATAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCa         <  1:810446/76‑1 (MQ=255)
                                                               tCGGTGCTCGAT‑CTGGAAAAAATTACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCa         <  1:1466676/76‑1 (MQ=255)
                                                               tCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCa         >  1:3415145/1‑76 (MQ=255)
                                                                cGGTGCTCGAT‑CTGGAAAAAATGCCCGTTGATGACATCATGGTGCCGCGCATTGAAATTATCGGTATTGATATCaa        <  1:1893263/76‑1 (MQ=255)
                                                                cGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCaa        >  2:1034561/1‑76 (MQ=255)
                                                                 ggTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAAc       >  1:1552993/1‑76 (MQ=255)
                                                                  gTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACg      <  1:2838934/76‑1 (MQ=255)
                                                                  gTGCTCGA‑CCTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCCGAGAAGTTATTGGTATTGATATAAAcc      >  2:3679106/1‑75 (MQ=255)
                                                                  gTGCGCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAAca      <  2:5708710/76‑2 (MQ=255)
                                                                    gCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACgct    >  2:5082198/1‑74 (MQ=255)
                                                                    gCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGAt    <  1:2520277/76‑1 (MQ=255)
                                                                     cTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATg   >  1:702347/1‑76 (MQ=255)
                                                                     cTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATg   <  2:3749853/76‑1 (MQ=255)
                                                                     cTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATg   <  1:4403871/76‑1 (MQ=255)
                                                                           |                                                                       
TTGCGCACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATCCTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGA  >  NC_000913/2747313‑2747459

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