Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,019,860 Δ2 bp intergenic (+59/‑24) b3836 → / → b3838 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,019,8600C.97.2% 301.3 / 0.8 71intergenic (+59/‑25)b3836/b3838hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (35/34);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01
*NC_0009134,019,8610C.97.1% 295.2 / 2.0 71intergenic (+60/‑24)b3836/b3838hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (34/34);  minor base G (2/0);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGCCGCCTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGC  >  NC_000913/4019793‑4019924
                                                                   ||                                                                
cAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                         <  1:3251861/76‑1 (MQ=255)
cAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                         >  1:3236161/1‑76 (MQ=255)
 gggTGTATGCCGTGTTTGATAGCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGcc                                                        <  2:584176/75‑1 (MQ=255)
 aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGcc                                                        <  2:555619/76‑1 (MQ=255)
 aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGcc                                                        <  1:196673/76‑1 (MQ=255)
  ggTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCg                                                       >  1:927193/1‑76 (MQ=255)
  ggTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCg                                                       >  1:982388/1‑76 (MQ=255)
   gtgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGc                                                      >  1:1025061/1‑76 (MQ=255)
     gtAATCCTTGTTTTATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                                    <  1:668938/76‑1 (MQ=255)
     gtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCCTCATCGGCCTCG‑‑GCGTTCTGGGGCCCCac                                                    >  2:3733720/1‑75 (MQ=255)
     gtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCac                                                    >  2:578217/1‑75 (MQ=255)
     gtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                                    <  1:1143430/76‑1 (MQ=255)
     gtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                                    >  1:3990131/1‑76 (MQ=255)
      tAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                                   >  2:589231/1‑76 (MQ=255)
      tAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                                   <  1:53469/76‑1 (MQ=255)
        aTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGa                                                 <  1:1682004/76‑1 (MQ=255)
        aTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGa                                                 >  2:505123/1‑76 (MQ=255)
         tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGAc                                                >  2:1032271/1‑76 (MQ=255)
           cGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTg                                              >  2:1112718/1‑76 (MQ=255)
            gtgtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGc                                             >  1:3155575/1‑76 (MQ=255)
            gtgtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGc                                             >  2:772030/1‑76 (MQ=255)
              gtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCt                                           >  1:3273630/1‑76 (MQ=255)
              gtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCt                                           >  2:2762757/1‑76 (MQ=255)
               tttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCtg                                          >  1:3372135/1‑76 (MQ=255)
                ttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCtgt                                         >  2:3797424/1‑76 (MQ=255)
                ttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCtgt                                         <  2:766460/76‑1 (MQ=255)
                ttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGCCCTCGGCGGTTTGGGTCCTCAACTACTGCCtgg                                         >  2:3834836/1‑75 (MQ=255)
                 tGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCtgtg                                        >  1:3181286/1‑76 (MQ=255)
                 tGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCtgtg                                        <  1:1112719/76‑1 (MQ=255)
                  gATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGCCTGCCTGCgg                                       >  2:716030/1‑76 (MQ=255)
                  gATAGCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                       <  2:3602957/76‑1 (MQ=255)
                   atatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGc                                      <  1:1511002/76‑1 (MQ=255)
                   atatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGc                                      >  2:2146428/1‑76 (MQ=255)
                   atatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGCCCGCAACGACTGCCTGTGGc                                      >  2:3833822/1‑76 (MQ=255)
                    tatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCg                                    <  1:4221223/76‑1 (MQ=255)
                     atCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                                   >  2:890775/1‑76 (MQ=255)
                     atCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                                   >  1:3614128/1‑76 (MQ=255)
                     atCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                                   <  1:3237733/76‑1 (MQ=255)
                     atCGGTTTGAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                                   <  2:1157755/76‑1 (MQ=255)
                      tgggTTTTAGCGAACGGCTATTGGTGTTCATCATCGGCCGCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGt                                  <  2:172120/74‑1 (MQ=255)
                       cGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgggg                                 >  1:4818264/1‑74 (MQ=255)
                       cGGTGTTAGCGAACGGCTATGGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTa                                 <  2:4078494/76‑1 (MQ=255)
                        ggTTTTAGCGAACTGCTTTTGGTGTTCATCATCGGCCGCGGCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaa                                <  2:7369/76‑1 (MQ=255)
                        ggTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgggga                                >  1:2132906/1‑73 (MQ=255)
                        ggTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaa                                <  1:2638067/76‑1 (MQ=255)
                        ggTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaa                                <  1:3323659/76‑1 (MQ=255)
                          ttttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGaaaa                              >  1:1776391/1‑76 (MQ=255)
                           tttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaaaa                             <  2:3631461/76‑1 (MQ=255)
                           tttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGaaaaa                             >  1:3570339/1‑76 (MQ=255)
                            ttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGaaaaaa                            >  1:657038/1‑75 (MQ=255)
                             tAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACg                           <  2:3129629/76‑1 (MQ=255)
                             tAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGaaaaaag                           >  1:4116726/1‑74 (MQ=255)
                              aGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACgg                          <  1:679393/76‑1 (MQ=255)
                              aGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGAAAAAAgg                          >  1:3137405/1‑76 (MQ=255)
                                ggAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCGGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTa                        <  2:4325069/75‑1 (MQ=255)
                                 gAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGTAAAAACgcag                        >  2:3056376/1‑73 (MQ=255)
                                 gAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGGAAAAAGGGAg                       >  1:1766899/1‑76 (MQ=255)
                                  aaCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGc                      <  1:2081295/76‑1 (MQ=255)
                                   aCTGCTATTGTTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCg                     <  1:1068401/76‑1 (MQ=255)
                                   aCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCg                     <  1:4470530/76‑1 (MQ=255)
                                       gtATTGGTGTTCATCATCTGCCTCG‑‑TCTGTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCt                 <  2:3570339/75‑1 (MQ=255)
                                       cTATTGTTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCt                 <  1:2044249/76‑1 (MQ=255)
                                        tATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGGAAAAAGGTAGggggggg                >  1:2687271/1‑73 (MQ=255)
                                        tATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGAAAAAAGGTAGCGGGGTg                >  1:2682848/1‑76 (MQ=255)
                                         aTTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGAAAAAATGGAGCGGGGTgg               >  1:469073/1‑76 (MQ=255)
                                            tgtggTCATCGGCGGCCTCGTCG‑‑TTCGGGGGCGGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGAtt            <  2:1408586/75‑1 (MQ=255)
                                             gtgtTCGTCTGCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTc           <  2:4596585/76‑1 (MQ=255)
                                                    tcagcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgtt    <  2:4089667/76‑1 (MQ=255)
                                                     catcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttg   <  2:2094208/76‑1 (MQ=255)
                                                     catcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttg   <  1:818128/76‑1 (MQ=255)
                                                      atcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttgc  <  1:1319699/76‑1 (MQ=255)
                                                                   ||                                                                
CAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGCCGCCTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGC  >  NC_000913/4019793‑4019924

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