Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,469 +G intergenic (+36/‑197) 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,4661.G98.4% 162.2 / ‑5.0 64intergenic (+33/‑200)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base G (27/35);  minor base T (0/1);  total (28/36)
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.95e-01

GCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTG  >  NC_000913/4008410‑4008534
                                                            |                                                                       
ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCAtt                                                         <  2:2041517/75‑1 (MQ=255)
ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCAtt                                                         >  1:1830542/2‑76 (MQ=255)
 cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTa                                                        <  1:3575133/76‑1 (MQ=255)
 cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTa                                                        >  1:3924921/1‑76 (MQ=255)
 cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTa                                                        <  1:4411708/76‑1 (MQ=255)
  tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTAt                                                       >  1:1965006/1‑76 (MQ=255)
   gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATg                                                      >  1:2470019/1‑76 (MQ=255)
    aaGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATgg                                                     <  2:2711712/76‑1 (MQ=255)
    aaGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATgg                                                     <  2:1089164/76‑1 (MQ=255)
          tgtattgCGTTGGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTg                                               <  2:2035636/73‑1 (MQ=255)
                    ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCcgcg                                     >  2:3397167/4‑76 (MQ=255)
                     tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCgc                                    >  2:1147751/3‑75 (MQ=255)
                     tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGt                                    <  2:1434903/74‑1 (MQ=255)
                     tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGt                                    <  2:2873027/74‑1 (MQ=255)
                       gtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt                                  <  1:2379453/76‑1 (MQ=255)
                         gACGGGATGAACGACGCCGAAATGCTGTCGACGGCGGCGAAAGGCTGCATTATAGGCAGCGCACACCAGCGTCtca                                >  1:3894669/1‑74 (MQ=255)
                          aCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTgca                               >  1:3810150/1‑74 (MQ=255)
                          aCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGaa                               >  2:3618591/1‑76 (MQ=255)
                          aCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGaa                               <  2:3435984/76‑1 (MQ=255)
                            gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAg                             >  1:2512601/1‑76 (MQ=255)
                            gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAg                             <  1:4318325/76‑1 (MQ=255)
                            gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGCGCGCACCACCGTCTGaaac                             >  2:898294/1‑75 (MQ=255)
                             gggtGAACGACGCGGAAATGCTGTGGGGGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGa                            <  2:4739062/73‑1 (MQ=255)
                             ggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGa                            <  2:1680721/76‑1 (MQ=255)
                              gATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAc                           <  1:4349596/76‑1 (MQ=255)
                               aTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCCGTGCGCACCAGCGTCTGCAACAcc                          >  2:1443325/1‑76 (MQ=255)
                               aTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAcc                          >  2:3123186/1‑76 (MQ=255)
                               aTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAcc                          >  2:4304880/1‑76 (MQ=255)
                                tGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCt                         <  1:3361916/76‑1 (MQ=255)
                                tGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCt                         <  2:3339859/76‑1 (MQ=255)
                                 gAACGACGCCGAAATGCTGTCGATGACGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCtt                        <  1:2092071/76‑1 (MQ=255)
                                  tgcgacgCCGTGATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCATTTGGCGCCAGCTTCTGAAAGACCTTc                       <  2:3167314/74‑1 (MQ=255)
                                     gacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                    >  2:554602/1‑76 (MQ=255)
                                     gacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGGGCGCACCAGCGTATGAAAGACCTTCCTc                    >  2:1316775/1‑76 (MQ=255)
                                     gacgCCGAAATGCGGGCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                    <  2:3667113/76‑1 (MQ=255)
                                      acgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                   >  2:2990825/1‑76 (MQ=255)
                                      acgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                   <  2:3045707/76‑1 (MQ=255)
                                       cgtcGAAATGCGGCCGTTGGGTGGGAAATGCTGTCTTATGGGCCGTGCGCACCAGCGTCTGACCGACCTTCATccc                  <  2:1932311/73‑1 (MQ=255)
                                       cgCCGAGATGCTGGCGAGGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACGTTCATccc                  <  2:972160/76‑1 (MQ=255)
                                       cgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                  >  2:2733047/1‑76 (MQ=255)
                                       cgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                  <  2:2481421/76‑1 (MQ=255)
                                        gCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCCGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                 >  1:3054200/1‑76 (MQ=255)
                                        gCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                 >  1:1076628/1‑76 (MQ=255)
                                         ccGAGATGCTGGCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                <  2:265884/76‑1 (MQ=255)
                                         ccGAATTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGTCCTTCCTCCCGa                <  1:1443252/76‑1 (MQ=255)
                                         ccGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                <  1:3326156/76‑1 (MQ=255)
                                         ccGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                >  2:2528308/1‑76 (MQ=255)
                                          cGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg               >  1:4465115/1‑76 (MQ=255)
                                          cGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg               >  1:4074132/1‑76 (MQ=255)
                                           gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc              <  1:109247/76‑1 (MQ=255)
                                           gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc              >  1:1466152/1‑76 (MQ=255)
                                           gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc              <  2:3924921/76‑1 (MQ=255)
                                           gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc              <  2:2818113/76‑1 (MQ=255)
                                           gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGCGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc              >  2:3773430/1‑76 (MQ=255)
                                            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    <  1:2434551/76‑1 (MQ=255)
                                                     cGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt    <  2:3796450/76‑1 (MQ=255)
                                                     cGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt    <  2:3835462/76‑1 (MQ=255)
                                                       aTGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg  >  1:3514535/1‑76 (MQ=255)
                                                            |                                                                       
GCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTG  >  NC_000913/4008410‑4008534

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