Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,473 +G intergenic (+40/‑193) 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,4721.G98.6% 186.2 / NA 73intergenic (+39/‑194)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (33/38);  minor base T (0/1);  total (33/40)
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.98e-01

GCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATG  >  NC_000913/4008410‑4008543
                                                                    |                                                                        
ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCAtt                                                                  <  2:2041517/75‑1 (MQ=255)
ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCAtt                                                                  >  1:1830542/2‑76 (MQ=255)
 cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTa                                                                 >  1:3924921/1‑76 (MQ=255)
 cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTa                                                                 <  1:4411708/76‑1 (MQ=255)
 cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTa                                                                 <  1:3575133/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:2873027/74‑1 (MQ=255)
                     tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGt                                             <  2:1434903/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:3435984/76‑1 (MQ=255)
                          aCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGaa                                        >  2:3618591/1‑76 (MQ=255)
                            gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAg                                      <  1:4318325/76‑1 (MQ=255)
                            gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAg                                      >  1:2512601/1‑76 (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:2481421/76‑1 (MQ=255)
                                       cgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                           >  2:2733047/1‑76 (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:1466152/1‑76 (MQ=255)
                                           gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc                       <  2:3924921/76‑1 (MQ=255)
                                           gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc                       <  1:109247/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)
                                                          gcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa        >  1:1469020/4‑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     >  2:2664109/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)
                                                                    |                                                                        
GCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATG  >  NC_000913/4008410‑4008543

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: